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Transonic Flow Past Wedge Sections 


H. W. LIEPMANN* anp A. E. BRYSON, JR.t 
California Institute of Technology 


SWMMARY 


The paper gives the results of an interferometric wind-tunnel 
gudy of transonic flow past wedge sections. Both 
fow with a local supersonic zone and supersonic flow with a de- 
tached shock wave have been investigated. 
fibution on these wedges and the drag coefficients as functions 
ofthe Mach Number have been obtained. 
pared with general theoretical considerations and recent theo 
retical work of Guderley and especially of Vincenti and of Cole. 
Based on theory and experiment, the behavior of the drag curve 
fom M = 0 through \J = 
petely determined. 


subsonic 
The pressure dis- 


The results are com- 


1 to supersonic flow can be com- 


SYMBOLS 


free-stream Mach Number 
local Mach Number 
free-stream velocity 
components of the local velocity 
critical velocity 
perturbation velocity component 
Ma, Ms! refers to specific characteristic values of the 
ma, ms ' Mach Number 
6 = wedge angle (semiapex angle) 
C wedge chord 
pressure coefficient 
drag coefficient based upon C, @ (i.e., pressure 
drag coefficient of the half wedge) 
ratio of specific heats 
Cre(y + 1)'78 eo ', reduced pressure coefficient 
Co(y + 1)'/3/6°/*, reduced drag coefficient 
(M? — 1)/[(y + 1)6]’”, transonic similarity 
parameter 
distance along the axis of the wedge measured 
from the apex 


(I) INTRODUCTION 


T SWE PAPER GIVES THE RESULTS of a detailed experi- 

mental investigation of two-dimensional trans- 
onic flow past wedge sections—i.e., profiles consisting of 
e————___ 


’ 


Presented under the title “Experiments in Transonic Flow’ 
at the Aerodynamics Session, Annual Summer Meeting, I.A.S., 
Los Angeles, July 12-14, 1950. Received August 9, 1950. 

* Professor of Aeronautics, Guggenheim Aeronautical Labora- 
tory. 


t Graduate Student. 


a simple wedge at zero angle of attack followed by a 
straight section. The significance of the data and the 
comparison with theory are discussed. 

Investigations of transonic-flow problems are inher- 
ently not easy. In theory, one has to deal with a prob- 
lem that cannot be handled with linearized equations 
without losing its essential features. In experiment, 
one is faced with the well-known difficulties of investi- 
gating the flow field near Mach Number 1—i.e., flow 
with extremely large lateral extent of the perturbation 
field of a body. In the last few years, progress in 
transonic-flow theory has been steady. The investi- 
gations of von Karman, Busemann, Guderley, and 
others have led to a general understanding of the main 
features of transonic flow. Similarly, the technique 
of wind-tunnel experiments, especially the develop- 
ment of variable Mach Number supersonic tunnels and 
of interferometric measuring techniques, etc., has 
advanced considerably in the last few years. 

During the course of experimental research of vis- 
cous phenomena and of boundary-layer shock-wave in- 
teraction carried out at GALCIT, under N.A.C.A. 
sponsorship, an interferometer was built. Before ap- 
plying the interferometric technique to a problem as in- 
volved as the boundary-layer shock-wave interaction 
problem, it was decided to try to investigate in detail 
an essentially nonviscous transonic-flow field of large 
dimensions. In view of the status of transonic theory 
and the lack of experimental results in this range, it 
was clear that a study of the simplest transonic-flow 
problem was of great theoretical and practical interest. 
Hence, an investigation of transonic flow past wedge 
sections was undertaken. The choice of the section, 
angles, etc., was based upon the following considera- 
tions: 

(1) Two-dimensional flow is preferable, since a 
theoretical study of two-dimensional flow is easier than 
axially symmetric flow and since, furthermore, the 
evaluation of interferograms is more direct and faster. 
Obviously, setting up of two-dimensional flow past a 
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body is not easy at all, but this difficulty has to be ac- 
cepted. 

(2) The extremely large lateral extent of the per- 
turbation field of a body in the transonic region makes 
it necessary to reduce size—e.g., the chord of the model 

relative to the wind-tunnel dimensions as much as 
possible. This, in turn, tends to reduce the Reynolds 
Number of the model and thus increase the importance 
of boundary-layer and, in general, viscous phenomena. 
Since this investigation was mainly concerned with a 
nonviscous-flow problem, it was obviously necessary 
to reduce the viscous effects as much as possible. The 
very nature of the flow does not permit the increase of 
size of the model with respect to a given wind tunnel. 
However, one can restrict the importance of viscous 
phenomena to a large extent by choosing the geometry 
of a model such that the boundary layer is in a region of 
favorable pressure gradient. Separation problems are 
thus avoided, and the importance of boundary-layer 
effects is reduced. . Based on this consideration, it was 
decided to investigate a wedge followed by a straight 
section rather than a double-wedge section. This 
choice of the model does not greatly restrict the investi- 
gation because the important part of a double-wedge 
section in nonviscous transonic flow for W/ > | is actu- 
ally the forebody. The flow over the afterbody is 
supersonic throughout and, hence, can be treated with 
fairly standard methods. 

Clearly, viscous phenomena cannot be completely 
eliminated in this fashion. For example, there is still 
a slight change of the effective shape of the wedge by the 
existence of the boundary layer, as well as a ‘rounding 
off’’ of the shoulder corner, etc. However, these ef- 
fects are comparatively slight compared with large 
effects such as separation. 

(3) Both theoretical advantages 
need focus the interest upon thin sections. 
angles chosen were, thus, 4.5°, 7.5°, and 10° semi- 
angles. It was not considered feasible to measure 
wedges with angles less than 4.5° semiangle due to 
structural difficulties. On the other hand, for small 
apex angles the transonic similarity law should allow the 
extension of data without additional measurements. 

Because of the precautions outlined in (2) and the 
use of interferometric measuring methods, it proved 
possible to measure the complete pressure distribution 
as a function of Mach Number on wedges of length 
not greater than '/, in. and to obtain reliable pressure 
drag data in spite of the fact that the Reynolds Num- 
ber of the sections was only about 6 X 10%. 

Curiously enough, there does not seem to exist much 
previous experimental work on wedge sections in the 
transonic régime. Interferometric investigations of 
transonic flow of this type, carried out at Braunschweig 
after the war and reported by Pack,' are not too ex- 
tensive and reliable, since the wind-tunnel flow was 
nonuniform, since the tunnel seems to have been choked 
during at least some of the measurements, and since 


practical 
The wedge 


and 
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only one supersonic Mach Number was investigate 
No other systematic experimental investigation of th: 
problem seems to have been published so far, 

Theoretical results consisted, at the beginning of th 
present research, essentially of von Karman’s gener] 
similarity considerations,’ of Busemann’s general gj. 
cussion of flows with detached shock waves,’ anq , 
papers by Guderley and by Guderley and Yoshihap 
respectively——one is a general discussion of transopj 
flow patterns,’ the other is a detailed computation ¢j 
flow past thin sections at a free-stream Mach Numbe 
of unity.° During the course of the present experi 
ments, it was learned that Vincenti,® at the Ames Aer 
nautical Laboratory, had carried out a detailed cop, 
putation of flow past double-wedge sections with 
detached shock wave. The cited theoretical papers 
have the following in common: that they deal with th 
problem in transonic approximation but are rigoroy 
within this frame. That is, no further simplifying as 
sumptions are introduced, and in Guderley’s anj 
especially Vincenti’s work the Tricomi equation jp 
the hodograph is actually solved under the prescribed 
boundary conditions. Hence, these papers differ con. 
siderably from others not cited here, where other sim- 
plifying assumptions have been made, such as an a- 
sumed shape of the shock, etc. 

Vincenti’s work has so far not been published 
However, because of the cooperation of the N.A.C.A, 
it was possible to exchange information between the 
Ames Aeronautical Laboratory and GALCIT, and, 
hence, a partial comparison of the present results with 
Vincenti’s theoretical work is included here. For 
further details with respect to the theory and compar: 
son with experiment, reference is made to forthcoming 
N.A.C.A. publications. 

In the last few weeks, J. D. Cole, of GALCIT,’ sue. 
ceeded in finding a solution of the transonic equations 
for a free-stream Mach Number less than unity, from 
which the drag of the wedge can be obtained. Henee, 
it proved possible to compare the measured drag co- 
efficients with theory both in the high subsonic region 
(Cole) and in the low supersonic region (Vincenti 
The cooperation of Vincenti and Cole in making it pos 
sible to use the unpublished results is gratefully ac- 
knowledged. 

The authors would also like to acknowledge the co- 
operation of the other members of the Transonic Re 
search Group at GALCIT. Especially the help of $ 
Dhawan and G. Solomon was essential. 


(II) APPARATUS AND METHODS 


(1) Wind Tunnel 


The measurements were carried out in the GALCIT 
4- by 10-in. Transonic Wind Tunnel. For a descrip- 
tion of the tunnel, and especially of the flexible nozzle 
employed, see reference 8. 
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2) Interferometer 

An interferometer of somewhat unusual construction 
was used in this investigation. The interferometer is 
described in detail in reference 9. 


3) Models 

Accurate lapped steel wedges were used. The gen- 
eral geometry of these wedge sections is shown in Fig. 1. 
Two pressure orifices symmetrically located on the sides 
of the wedge were used both to set the angle of attack 
«o zero and to identify one of the interferometer fringes. 
4n orifice near the shoulder was used as a second check 
on the measured Mach Number and fringe identifica- 


tion. 


4) Two-Dimensional Flow 


To obtain two-dimensional conditions in the wind 
tunnel as nearly as possible, the span of the wedge sec- 
tion has to be carefully chosen. If the span is too long 
and extends well into the boundary layer on the wall, 
strong shock-wave boundary-layer interaction on the 
sidewalls results. If, on the other hand, the span of 
the wedge is chosen too short, one is faced with three- 
dimensional tip effects. Between these two limits 
there exists an optimum ratio of span to tunnel width. 
The optimum was found by gradually shortening the 
span of a given wedge section until the shock wave had 
reached its best definition as judged by both schlieren 
and interferometer observation. For the present tun- 
nel, this turned out to be a ratio of model span to tunnel 
width of 0.875. The Ames Aeronautical Laboratory 
has made shadow photographs of tips of wedge sections 
with detached shock waves near the boundary layer on 
aflat wall. These shadow photographs were taken at 
right angles to the flow direction and the span so that 
the plane of the sections and the wall boundary, as 
well as a section through the shock, became visible. 
These photographs, which thanks to the cooperation 
of the N.A.C.A. were made available to the authors, 
confirm the fact that there exists an optimum tip dis- 
tance from the wall such that very nearly two-dimen- 
sional flow is obtained over the section. More details 
and further checks will be reported later in N.A.C.A. 


publications. T 


5) Wind-Tunnel-Wall Interference 


For the measurements of flows with detached shock 
waves, the interference of the top and bottom walls is 
zero, since the subsonic zone was kept well away from 
the walls. In the measurements in subsonic flow re- 
ported here, some interference must exist. However, 
the ratio of the height of the tunnel to the chord of the 


t It should be remembered that a profile that does not touch 
the sidewalls is not equivalent to a three-dimensional body with 
aN aspect ratio given by its span chord ratio but is more nearly 
equivalent to an infinitely long profile with periodic gaps equal 
to twice the distance of the edges from the sidewalls. 
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airfoil is so large that the effect is minor and consists 
mainly in a slightly altered free-stream Mach Number. 
This latter effect was taken into account. 


(III) CHARACTERISTIC FEATURES OF TRANSONIC FLOW 
PAST WEDGE SECTIONS 


This is not the place to discuss the transonic problem 
in general; however, it appears useful and necessary 
to draw attention to some specific and important fea- 
tures of transonic flow past single-wedge sections. 
Some of these characteristics have been discussed be- 
fore; some have, to the authors’ knowledge, not been 


noticed or explicitly stated. 


(1) Characteristic Values of the Free-Stream Mach 

Number 

For a wedge section of half apex angle @ in uniform 
parallel flow, there exist the following characteristic 
Mach Numbers: 

(a) M,, the free-stream Mach Number at which 
the shock wave attaches. 

(b) Ms, the Mach Number at which the flow behind 
the attached shock wave is exactly sonic. 

(c) Mo, the Mach Number corresponding to the so- 
called Crocco point—i.e., the Mach Number at which 
a characteristic change in the curvature of the attached 
shock wave occurs. 

Mo is not too important for the present investiga- 
tions. Reference is made to Guderley’s paper‘ for a 
more detailed discussion of J». 

(d) The critical Mach Number ./-—1.e., the free- 
stream Mach Number for which the local Mach Num- 
ber first reaches unity—1is zero for a wedge due to the 
sharp corner at the shoulder. 

Hence, neglecting viscosity and boundary-layer 
effects, the transonic region for a wedge section extends 
from M = 0to M = Ms. 
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(2) Characteristic Values of the Local Mach Number 


The local Mach Number at the apex of the wedge and 
the Mach Number immediately before and behind the 
shoulder corner have characteristic values in transonic 
flow. 

(a) The Mach Number at the apex of the wedge 
m,, is always zero, since a stagnation point exists there. 

(b) The Mach Number at the shoulder just before 
the turn, ms,, is always equal to unity. 

(c) It is impossible to reach the free-stream Mach 
Number on the flat part of the wedge afterbody by an 
isentropic Prandtl-Meyer expansion around the corner. 
The flow around the corner overexpands first and re- 
compresses through one or more shock waves. 

Statement (a) is trivial; (b) and (c) follow directly 
from simple considerations: For supersonic flow with a 
detached shock wave, the sonic point must be at the 
shoulder of the wedge; if the sonic point lies between 
the apex and the shoulder, then the flow behind the 
sonic point would be supersonic and the shoulder could 
not influence the forebody and the shock wave. Since 
for a given wedge angle the distance of the shoulder 
from the apex—i.e., the chord C—is the only charac- 
teristic length of the problem, it must determine the 
distance of the detached shock wave, etc. Hence, the 
shoulder must influence the forebody, and the sonic 
point cannot lie ahead of the shoulder. Furthermore, 
subsonic compressible flow around a sharp corner is 
impossible. This statement is plausible, since, here, 
incompressible flow would reach infinite velocities and 
compressible flow at least supersonic velocities. Math- 
ematically, this statement is shown to be true by ob- 
serving that, if subsonic flow around a sharp corner 
were possible, the corner point would map into a line 
in the hogodraph plane. This requires that the 
Jacobian of the transformation vanish, which is im- 
possible for subsonic flow (see, e.g., reference 10). In 
supersonic flow, this latter fact is possible and leads to 
the simple Prandtl-Meyer solution, in which the turn- 
ing angle from a given Mach Number determines the 
flow uniquely. Statement (b) may be proved for sub- 
sonic flow by a comparatively simple argument. 

It is also easy to show that the sonic line leaves the 
wedge shoulder normal to the direction of the flow—i.e., 
normal to the flanks of the wedge—and that the sonic 
line begins with infinite curvature. 

Concerning (c), one may remark that, if completely 
isentropic recompression were possible, for transonic 
flow with a subsonic free-stream Mach Number the 
drag would be zero. However, this does not seem rea- 
sonable. 

One can argue in more detail as follows: According 
to (b) the Mach Number ahead of the turn is equal to 
unity for the whole range of free-stream Mach Num- 
bers from 0 up to Ms. The total angle of turn is al- 
ways 6. However, Prandtl-Meyer flow starting from 
sonic velocity will lead to a Mach Number ms,, which 
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is uniquely determined by @. 
ms, # M; for example, in subsonic flow obviously ms, > 
M. Consequently, the boundary conditions cay onh 
be met by a combination of expansion and shock waves 
It appears theoretically difficult to find the exact for, 
of the shock waves. One may prove, though, that at 
least one wave has to originate obliquely from the coy. 
ner and that an entropy integral extended over th, 
shock or system of shocks must be equal to the tot, 
pressure drag over the forebody.'” 


4 . 
Clearly, in general 


Experimental observations of shock formations pea 
the corner are strongly influenced by boundary-laye 
effects and are not easily evaluated. 


(3) Transonic Approximation and Transonic Similarity 


In the so-called transonic approximation, one cop. 
siders the difference of any Mach Number in the field 
from unity as a small parameter. This implies that 
the wedge angle @ is also small, and one keeps in this 
approximation, in general, only the largest terms jy 
dand in M? — lorm? — 1. 

Based upon this approximation the following results 


are easily obtained: 

(3/4) [(y + 18] (1 

ms,2 — 1 = (3/2) [(y + 16] (2 
Ms? — 1 = 2” [(y + 10] (3 


M,?-—-1= 


Furthermore, the following relations obtained in the 
same approximation are useful for later applications 
If MJ denotes the Mach Number ahead of a normal 
shock wave, m, the Mach Number behind, then 


M?—1=1-—m,? (4 


Furthermore, let 1/ be the stream Mach Number and 
m the local Mach Number at a certain point, and as- 
sume again that 


and 


m*—1 =p << ! 


then, we have, for the pressure coefficient at this point, 


—— — = Cp= : (u — v) (9 
(1/2)p.U.? y¥+ 1 
The transonic similarity rules are usually obtained from 
the differential equation of motion and the boundary 
conditions.2. Eqs. (1) to (5) are consistent with simi- 
larity rules. 

Eqs. (1), (2), (3),and (4) suggest the use of (7*.— 1)+ 
[((y + 1)0]” as independent variable and (m? — 1) + 


[(y + 1)@]* as dependent variable and, hence, sug: 


gest the existence of a functional relation of the form 
m? — | _73 M*—1 }) 
(y+ el ~ Uy + als 


From Eq. (5) for Cp, it follows that 
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Cey + 1)” j Wwz-1 ) dCp 2 dv 
2), =g \ "1 f (6) = 1— 

0 (vy + 1)8] dp y¥+1 dp 


since Cp, the drag coefficient, is simply equal to the 
rage pressure coefficient times the angle #, we have 


Croly + 1)” aw M?—1 \ 
6° U(y + 1e]?s 


Eqs. (6) and (7) are the well-known similarity laws of 


ave 


(7) 


transonic flow. 

In the following, Cp(y + 1)”'/0 will be denoted 
by Cp and (M? — 1)/[(y + 1 )0]’* by & It is also pos- 
sible to determine the form of the characteristic refer- 
ence length / of the shock wave—e.g., the length from 
the centerline to the sonic point on the shock wave. 

The drag coefficient Cp is related to the entropy 
change across a shock and, thus, to the Mach Number 
ahead of the shock wave through a relation of the form 


Cy ~ [CMP — 1)87(y + 1)?] C/C) (8) 
valid for transonic flow (e.g., reference 12). 
Thus, by comparing Eq. (8) with Eq. (7), we find 
(M? — 1)? 1 ” j ; ~~ } 
~™ i 7 . 2 > 
Grif c Gv+D Uy + 1Ieys 
(/C) ~ [(y + 1)0] (1/E*)h (*) 
or 
l ) a=] 4 
(y+ 1ée]’* =k yg = RE) (Y) 
‘ | Ul(y + els : 


4, Slope of the Drag Curve at \/ = 1 


A simple and intuitive consideration enables one to 
determine the slope of the Cp versus .V/ relation at 1/ = 
1, This slope will be found useful in joining Guder- 
ley’s Cp value at J = 1 with Vincenti’s drag curve 
for J > 1 and for a comparison with experiment. The 
reasoning is not restricted to this particular case but 
should apply to any profile for the part ahead of the 
sonic point. 

Eq. (3) expresses the well-known fact that the Mach 
Number downstream of a weak, normal shock wave is 
a much below unity as the Mach Number upstream 
That is, there exists a certain sym- 
For ./ very near unity, the de- 


is above unity. 
metry about AJ = 1. 
tached shock wave is far away from the wedge and is 
nearly normal; the Mach Number just downstream of 
the shock is slightly subsonic. For the pressure dis- 
tnbution on the wedge, it is thus irrelevant whether 
this subsonic Mach Number at large distances is due 
to the presence of a shock wave or due to the fact that 
the velocity at infinity is slightly below sonic. Conse- 
quently, the pressure or Mach Number distribution on 
the wedge should be independent of Mach Number in 
the neighborhood of MJ = 1. To find the slope of the 
drag curve, one has only to differentiate Eq. (4) with 
fespect to u = Af? — 1; thus, 


The above consideration implies that vy = m* — 1 is 


independent of » as MW — 1—+.e., if u — O—hence, at 


M = 1, we find 

(dCp/dy), = 9 = 2/(¥ + 1) 
Transforming ‘this into the proper transonic similarity 
form, we find, for the drag coefficient, 


_ (=*) a 2 = 
a dy un =0 - Y +1 


(10) 


Co(y + 1)” 
d > 
l 6”? 
y+ 1 M? — 1 
¢ 2 
l(y + 1)6]” »=0 


and, consequently, 
dCp/dt = 2 (11) 


Actually, the measurements at high subsonic and low 
supersonic velocities indicate that the Mach Number 
distribution over the forebody is independent of the free- 
stream Mach Number not only as 1J — 1 but even for 
finite small differences from AJ = 1. The same result 
was also recently found by Bleakney and Griffith" using 
the Princeton shock tube. The recent theory of Cole 
furnishes essentially the same result analytically and 
leads to the correct slope [Eq. (11)! as @— 0.t 


(IV) GENERAL REMARKS CONCERNING COMPARISON 
OF EXPERIMENTAL DATA WITH APPROXIMATE THEORY 


A short discussion is included here about a_phil- 
osophy of comparing experimental data with approxi- 
mate theories. Recently, extensive discussion has 
arisen in connection with the application of small per- 
turbation theory of axial symmetric flow to experi- 
mental results. Especially, Van Dyke" has called at- 
tention to certain difficulties in choosing the proper 
parameters in an approximate theory. The following 
discussion is related to Van Dyke’s arguments but 
differs slightly in the point of view. The transonic 
theory referred to here involves two different kinds of 
One is a physical approximation— 


approximations. 
The second is a 


the assumption of nonviscous flow. 
mathematical approximation that involves the use of 
i.e., the difference 
i.e., the 


the highest order terms in 1/* — 1 
of the Mach Number from unity 
deflection angle. 

It is the second type of approximation which will be 
These approximations make certain 


and in @ 


discussed here. 
functions of the flow field equivalent for the theory— 
i.e., within the approximation certain functions of the 


+ In Pack’s report,! reference is made to an unpubiished paper 
of Macoll and Codd, and it is indicated that similar considerations 
However, the authors have so far been 
G. Guderley 


have been made there. 
unable to obtain a copy of Macoll and Codd’s work. 
has independently made similar considerations recently. 
derley’s result differs from the one given here by a different defi- 


Gu- 


nition for Cp. 
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Fic. 3. Infinite fringe interferogram (detached shock). 
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TRANSONIC FLOW 
For example, (P — P..) + 
(1/2)p«U*, —2(u'/a*), —2(u’/U) are equivalent 
within the approximation; hence, either of the last 
two expressions may be taken to define Cp. Simi- 
larly, M? — 1 and 2(.M/ — 1) are equivalent, and, hence, 
the sinilarity parameter has the equivalent forms 
(M2 — 1)/0" or (MVM — 1)/0 *, The experimental 
data of this study, and indeed practically all useful 


velocity, etc., are identical. 


data, are obtained at values of 6 and MW for which 
the differences between equivalent expressions in the 
theory are not always negligible, and the question arises 
of how theory and experiment are to be evaluated. 

(1) One may consider that the approximate theory 
js correct in the limit as 1J ~ 1,@— 0. This point of 
view is essentially implied if one uses mathematical 
order symbols as has been done extensively by the 
British school in linearized supersonic-flow theory. 
For the experimenter, this means that he should ob- 
tain data for a set of conditions for several 86 and M and 
then establish the limiting law by extrapolating to 
#-0and W— 1. This latter result must then agree 
with the approximate theory. 

(2) One may also look at the results of the theory as 
valid for small but finite values of MW? — 1 and @. 
Then the theory does not provide a definite result (as 
for example a drag curve) but rather a region, more or. 
less narrow, of “‘equivalent’’ drag coefficients. It is 
not possible to distinguish between these within the 
framework of the theory. For example, such a region 
would be provided if Cp were alternatively defined by 
—2(u'’/a*), —2(u'/U), or (P — P.)/(1/2)p.U”. 
The experiments should then fall within the limits of 
a band and can be used to decide which one of the 
equivalent expressions is the best in the sense that it 
extends the validity of the approximate theory over the 
widest possible range of the experimental parameters. 
In some cases a decision of this type can also be fur- 
nished by comparison with other theoretical results if 
the regions of validity overlap. Here, it can be seen, 
for example, that 1/* — 1 is a more reasonable choice 


than 2(.7 — 1), since with 1/* — 1 the transonic theory * 


goes over smoothly into Ackeret’s thin airfoil theory. 
This has already been shown by von Karman.’ 

The above considerations led to the presentation 
chosen for the results of this investigation. It was 
found difficult to adopt the first point of view—i.e., to 
extrapolate the data to the limit 9 0. This difficulty 
is essentially due to the nonlinear approach to this 
limit and to the experimental scatter. 

The second point of view was accordingly adopted in 
presenting Fig. 13. Here, the theoretical indeter- 
minacy is indicated by showing the range of Cp’s be- 
tween the limits defined by —2(u’/a*) and —2(u’/U). 
It is clearly seen how the precision of the theory in- 
creases as 9 — 0. The same convention for the pres- 
entation of supersonic theoretical results is adopted in 
Fig.11. The same point of view can be adopted in pre- 


PAST 





WEDGE SECTIONS 751 
senting the slope at M = 1.¢ dCp/dy can, of course, be 
obtained also using the exact form of the pressure co- 


efficient. In this case one has, instead of Eq. (10), 


(=) 24, Ew 
du n=0 Y 77 1 


: Mew. ‘i 
1+ (m? — » | 
y+ | 


In first approximation with m? — 1 = » << 1, Eq. 
(10) results. Again, the approximate theory cannot 
distinguish between Eqs. (10), (11), or, for example, 


(<) ; 2 ( v ) 
= i+ 
du 0 ; i oe YT 1 


u 
Similar results obviously hold for Cp. Thus, the slopes 
of Cp at M = 1 are presented as fans in Fig. 13, the 
larger slope corresponding to Eq. (10) and the smaller 
slope to Eq. (13). The difference is seen to be very 
small indeed for small 6. 


(vy - ad 


f (12) 


(13) 


(V) RESULTS OF THE MEASUREMENTS 


(1) General Character of the Flow Field 


Figs. 2 to 4 show typical interferograms of the flow 
phenomena using the so-called “‘infinite fringe’ picture. 
The fringes shown correspond to lines of constant 
density and, thus, nearly to lines of constant Mach 
Number. The difference between successive fringes is 
approximately 0.05 Mach Number. For the actual 
evaluation it is more convenient to use the superposi- 
tion method,’ using a finite fringe picture such as Fig. 5 
and a finite fringe picture of the undisturbed flow. The 
result is a picture as shown in Fig. 6. The lines of 
constant density appear here again. From _ photo- 
graphs such as these, it is not difficult to evaluate the 
pressure field, determine the sonic lines, etc. In all 
these evaluations the increase in entropy through the 
shock wave was taken into account. 

Fig. 7 shows typical tracings of the whole flow field 
about a wedge at different values of the free-stream 
Mach Number, obtained from superimposed finite 
fringe exposures. 


(2) Pressure Distributions 


The pressure distributions on one wedge section as 
functions of the free-stream Mach Number are shown 
in Figs. 8a and 8b. The slow and continuous variation 
from a typical subsonic to a typical supersonic dis- 
tribution is clearly seen. For high subsonic free- 
stream Mach Numbers and for supersonic Mach Num- 
bers with the bow wave detached, the sonic velocity is 
reached a little ahead of the corner. This is due to the 
boundary layer, which effectively rounds off the corner. 
After ‘‘attachment’’ of the shock wave, there exists a 
small subsonic zone around the very apex of the wedge. 
Clearly, the presence of the boundary layer distorts the 


+ This was pointed out to the authors by W. G. Vincenti 





7 


EO) 


O42 





JOURNAL 


OF 


THE 


Hi} 
/ fyi} 


hi} 


4j 


AERONAUTICAL 


Hif|ii (i 
Hy 


} 


ij 


ij 
pit 


ii) 





47 





AAA 
BAL } 


ALE 


~ 
~ 
~ 
~ 
— 
~~ 
~ 


TH 





Fyc. 7. 


Fic. 6. Typical superimposed finite fringe interferogram. 
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1 = 1.278; 1/2 angle of wedge = 10°. 
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sured pressure distribution. This is evident from 
attached shock plots of Fig. 8a and also from a 
comparison of one of the measured pressure distribu- 
tions with one computed by Vincenti, Fig. 9. The 
order of magnitude of the change in pressure distribu- 
tion and the qualitative behavior can be obtained from 
naive considerations of boundary-layer displacement. 
The layer is laminar, and, since the Reynolds Number 
based upon the chord is approximately 6 X 10*, one can 
write for the displacement thickness 6*/C at X/C 


mee 
the 


approximately 
6 lyr 
: - LX 


fo ~ 0.007 ¥7, 
C 
6 103 ( 
tic xX 


The additional pressure change thus corresponds 
roughly to the change produced by a paraboloidal body 
with about 0.7 per cent maximum thickness at X = C. 
This turns out to be of the right order of magnitude 
as the observed deviations, and, hence, it seems rea- 
sonably certain that most of the deviations from the 
theoretical shape are due to boundary-layer influence. 
Since the additional pressure distribution is nearly 
antisymmetrical fore and aft, the drag is hardly in- 
fluenced. 

The influence of the expansion around the shoulder 
edge is felt upstream due to the boundary layer. Hence, 
sonic velocity is reached a little ahead of the corner. 
Experiments on shock-wave boundary-layer interaction 
have indicated that the distance of upstream influence 
of a compression in the laminar boundary layer is of the 
order of 20 to 40 times the boundary-layer thickness. 
For an expansion one expects less, since the boundary 
layer becomes thinner because of the gradient. In 
the present experiments, 


e 


0 ” 
ee ~ 2X 1 
& V6 X 10¢ 


Hence, to find the sonic point 10 per cent ahead of the 
corner amounts to an influence through the boundary 
layer of roughly 56. Fig. 10 shows the movement of the 
sonic point with J/ as compared with the inviscid 
theory for a wedge with @ = 7.5°. 

Fig. 11 shows the local Mach Number m at the 
mid-chord point of a wedge with @ = 7.5° as function 
of the free-stream Mach Number /. The theoretical 
results of Vincenti and Guderley are shown, as well 
as the slope at 1/ = 1, according to the general con- 
siderations outlined earlier. The local Mach Number 
appears to have a constant value over the range from 
about J = 0.8 to 1.2, Fig. 12 demonstrates the same 
effect by showing selected local Mach Number distri- 
butions for various M. 


t The pressure distribution shown corresponds to -the-mest 
unfavorable case, since the relative thickness of the boundary 
layer to the wedge is obviously largest for the smallest wedge 
angle. 
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(3) Drag Coefficient 


Fig. 13 presents the variation of the drag coefficient 
versus free-stream Mach Number for the three wedge 


angles. The 
and Vincenti are 


theoretical 
shown. 


values 
The 


results of 


Cole, Guderley, 


Cole and 


Vincenti are presented as vertical lines rather than as 


points. 


For M > 1, the upper limit corresponds to 


evaluating Cy on the basis of taking Cp = —2(u’/a*); 
the lower limit of the lines corresponds to taking Cp = 


—2(u'/U); the reverse is true for M < 1. 
ley’s value at \/ = 1, these two limits coincide. 


For Guder- 
Plot- 


ting these results in this fashion shows clearly that the 
experiments and theory agree within the limits expected 
and that the theory becomes rapidly more precise for 
The comparatively large scatter 


smaller wedge angles. 
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of the experimental drag measurements at subsopip 
speeds is due to the fact that the drag from the Pressure 
distribution appears as the small difference betwee 
fairly large numbers. 

Finally, Fig. 14 shows the drag coefficients plotted in 
the fashion of von Karman’s similarity rules. [t i 
clearly seen that the measurements for the thre 
wedges now lie much closer than before (Fig. 13) py 
that systematic deviations still exist. This is obviously 
due to the fairly large values of @ in these measurements 
In cases such as these, the similarity rules can be ey. 
pected to apply only approximately. 


(VI) CONCLUDING REMARKS 


The present measurements of pressure distribution 
and drag coefficient of simple wedge sections in trap. 
sonic flow confirm, within the limits of accuracy of both 
theory and experiment, the theoretical results of yon 
Karman, Guderley, and especially the recent work of 
Vincenti and Cole. 

The curve of pressure drag coefficient versus Mach 
Number for a simple wedge section rises rapidly at 
high subsonic Mach Numbers and reaches a maximum 
at a slightly supersonic Mach Number. Beyond this 
maximum the drag coefficient decreases as expected. 

The pressure or Mach Number distribution on the 
wedge has a stationary value at J/ = 1, and, hence, 
the slope of the drag coefficient-Mach Number curve 
can easily be determined at / = 1. In transonic simi- 
larity parameters, dCp/dé = 2 at M = 1. It is found 
experimentally that the Mach Number distribution 
remains, in fact, practically constant through M = 1. 
The local Mach Number at a specific point on the 
wedge increases with .\/ in subsonic flow, then remains 
constant through 1/ = 1, and increases again in the 
supersonic régime. The range in which the local Mach 
Number remains constant is for thin wedges given 
approximately by 


—1/2 < (M? — 1)/[(y + 10] < 1/2 
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First U. S. National Congress of Applied Mechanics 


Plans have now been formulated to hold the First U.S. National Congress of Applied Mechanics on 
The Illinois Institute of Technology will bz host to the Congress. 


Members of the Institute, one of the sponsoring societies, are cordially invited to submit papers, 


which should constitute original research in applied mechanics. 
words or the equivalent in equations, tables, and diagrams.) Abstracts must be submitted to the Chair- 


man of the Editorial Committee before April 14, 1951; complete papzrs must be submitted prior to June 


Additional information about the Congress may be obtained by writing to N. M. Newmark, 
Secretary of the Congress, University of Illinois, Urbana, III. 


(Their length must not exceed 5,000 


All members are invited to attend. 
































_ A Study of the Response of an Airplane 
Landing Gear Using the Differential Analyzer 


WALTER C. HURTYt 
Unwersity of California at Los Angeles 


SUMMARY 


The response of a conventional airplane landing gear is studied 
by deriving and solving the differential equations of motion for 
the system. The resistance of the air-oil shock strut depends 
upon an orifice damping coefficient and an exponent, y, which 
determines the air pressure as a function of volume. A study 
of experimental ‘‘drop-test”’ results on similar systems is made 
to determine these constants. The nonlinear equations of 
motion are solved by means of the differential analyzer for (1) 
the landing condition as simulated by a drop test and (2) the 
taxiing condition in which a periodic ground contour function is 
assumed. The differential analyzer solutions for the landing 
condition are compared with drop-test results on the same system. 


SYMBOLS 


components of displacement, ve- 
posi- 


, te 8 = horizontal 
locity, and acceleration, respectively; 
tive if directed in direction of motion of the 
airplane 

Ss 3 = vertical components of displacement, velocity, 

and acceleration, respectively; positive if 
directed up 

= angle included between the shock-strut center- 

line and the vertical 

@ = angle included between a tangent to the ground 
contour at the point of contact and the hori- 


zontal 

M = one-half of total airplane mass less m 

m = mass of “unsprung system’ including wheel, 
tire, axle, and shock-strut piston 

Ps = axial force applied to shock strut 

G(y2) = normal ground force applied to tire 

So’, fo" = derivatives of the ground contour function taken 


with respect to x¢ 
Lb = one-half of total aerodynamic lift on airplane 
acceleration of gravity 


£ = 

X = wave length of ground sinusoid 

a = amplitude of ground sinusoid 

w = angular frequency of excitation 

Pa = air pressure in shock-strut cylinder 

Ap, = pressure drop across orifice 

Py = fluid pressure in shock strut below orifice 

d., dy, dp = diameters, respectively, of shock air chamber, 
fluid chamber at orifice, and piston 

a,, ay, Gp = cross-section areas corresponding to diameters 
d,, dys, and dp 

"1 = displacement of shock-strut piston measured 


from fully extended position 


Received January 3, 1950. 
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thank Prof. L. M. K. Boelter, Dean of the College of Engineering, 
University of California at Los Angeles, for making available the 
use of the differential analyzer for the solution of the equations. 

{ Assistant Professor, College of Engineering. 


Vo = tire deflection 

F(y1) = force due to air pressure; function of y, 

R(m) = damping function due to flow through orifice 

Va = air volume in strut cylinder; function of y 

V, = air volume when piston is fully extended 

y = constant for polytropic compression and ex. 
pansion of air 

Nir yr! = initial values of y, and y. at beginning of motion 

Q = quantity of fluid flow through orifice 

, = orifice area 

Ca = orifice discharge coefficient 

Cy = orifice pressure coefficient 

ky = shock-strut constant = a,/V, 

h = height of free drop in drop test 


Subscripts A, B, and C refer, respectively, to the airplane center 
of gravity, the center of the wheel, and the point of contact of 
wheel with ground. Subscript so refers to snubber orifices. 


INTRODUCTION 


pling DESIGN OF AIRPLANE STRUCTURES, it is becom- 
ing increasingly important to take into considera- 
tion transient states of stress, as well as steady-state 
stresses resulting from periodic forces applied to the 
airplane. An example of the emphasis placed upon the 
determination of such dynamic stresses may be seen in 
the many efforts that have been made in recent years 
to determine stresses in the wings of an airplane re- 
sulting from the landing impact.’ * The growing ur- 
gency of this problem has been brought about by two 
trends in airplane design: first, the use of thin wings 
demanded by aerodynamic considerations connected 
with high-speed flight; second, the disposition of items 
of considerable mass along the wing span outboard of 
the landing gear attachment points. Both factors 
reduce the natural frequencies of the wing in both 
bending and torsion and, hence, increase the likelihood 
of exciting large amplitude vibrations by means of 
ground forces. 


The calculation of bending and torsional vibrations 
of an elastic wing can be made by well-established 
methods of analysis provided the forces applied by the 
landing gear are known as functions of time. It is the 
purpose of this paper to outline a method of deter- 
mining such force functions in specific cases rather than 
to investigate the wing vibrations. Strictly speaking, 
the problem cannot thus be separated into two parts, 
since the motions of the shock-absorbing system and 
those of the vibrating wing are coupled. However, it 
has been shown? that the landing-gear forces can be 
determined with but a small error on the conservative 
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Fic. 1. 


Diagram of the system. 


side by considering the airplane as a rigid body, and this 
procedure is followed in the present investigation. 

The work reported here was carried out as part of a 
general analysis made in connection with the design of 
aparticular airplane. In scope, the problem consisted 
of the determination, by analytical means, of the 
forces imposed on the wing by the landing-gear strut in 
landing and taxiing the airplane. Since the bending 
modes of vibration were of particular interest, only the 
vertical forces transmitted to the wing were found. 
The equations of motion derived for the two degree of 
freedom system are nonlinear and require for their 
solution the use of numerical analysis techniques or an 
analog computer. In this case, use was made of the 
mechanical differential analyzer located on the campus 
of the University of California at Los Angeles and cper- 
ated by the College of Engineering. This machine is 
extremely versatile and well adapted to the solution of 
nonlinear differential equations. 


EQUATIONS OF MOTION 


The system, shown in Fig. 1, consists of a mass JJ, 
representing the airplane; a conventional air-oil shock 


/ * 
j(cos ¢ + 2’ sin g) + jeV1 + ae’? = Faae’ + wo" a 


2. 9 *s . a oh ot ~ 
Niece sin g + Po V1 + g0'? = a2’ + zot( ss Vi sin ¢ 


AN AIRPLANE 


LANDING GEAR 


strut; and a mass, m, representing the wheel, tire, axle, 
and shock-strut piston. Equations are derived for the 
motion of the system as the airplane rolls over an uneven 
ground surface at a constant velocity. From these 
equations we may write those pertinent to the landing 
case without difficulty. Energy dissipated by the tire 
is neglected, inasmuch as it will be small compared with 
that dissipated by the shock strut. Referring to Fig. 1, 
we may obtain the following velocity relationships 
useful in the derivation: 


X4 =A%c +wnsing + sina 
$4 = 2¢ — 1 cos gy — Wo cos a 


p . ts (1) 
Xp = Xe + Jo sin a 


33 = 2c — 32 COS a 


Differentiating with respect to time, the following ac- 
celeration equations are written: 


t, = % + ji sin g + je sin a + joa cosa 
34 = 2% — VW CcOsS y — Ve cos a+ Voc sina (2) 
Xp = X¢ + Je sin a + Jo COS @ - 


— j2 cos a + Sod sin a 


= ~ 
“B oc 


The equations of equilibrium of the vertical components 
of the forces at points A and B are 


—Mz, — Mg+L+ F,cos¢ = 0 (3) 


—mzZ, — mg + G(ye) cosa — F,cosg = 0 (4) 


We substitute for 3, and 2, their equivalents given by 
equations and express the equilibrium equations in 
terms of the vertical components of the accelerations 
at the ground contact point C. This, in turn, may be 
calculated if we know the velocity of the airplane 
and the ground contour function. 


Sc = Xc tana (5) 
Sco = Xc tan a + Xca sec? a (6) 
tan a = 2%’, sec?a = 1 + 2?) (7) 
” * 4 ( 

& = feo" (1+ 2e']e 


Finally, expressing %¢ in terms of x, from Eqs, (1) 
we can write the two equations of equilibrium as fol- 


lows: 











g oo it (9) 
_ ——== — cos 
mv 1 + 2"? m ’ 


















































758 JOURNAL OF THE AERONAUTICAL SCIENCES—DECEMBER, 1950 
These equations are perfectly general in that they apply gDetlection- inches ee 
for any ground contour function we may choose as 8 
long as the tire is tangent to the ground at not more Inflation press. 205psi 
; : « 225 psi— — 
than one point. Although no fundamental difficulty , Saneas-. ee 
: ‘ a . be 7 
is encountered in attempting to solve the equations by Se 
the differential analyzer, they may be simplified for we” 
most practical cases without appreciable loss of ac 2 = 6 
curacy. In the present problem, the strut angle ¢ is ee meen 
small for the taxiing attitude of the airplane. Also, 
the ground slope 2,’ is usually small compared with © 5 
unity. Thus, Eqs. (8) and (9) may be simplified to ls 
: ; ieiiieoe 
the following form: 225 psi— 2 
8 4 245 psi- _ ae a 
© il sand 
Ay + Ve = Eee’ + S720" + g[l _ (Z; Mg)| - ae 7 —<— 4 
(F,/M) (10) —" - 
2 == 
jo = Rate’ + X40" + g — {[Glye) — Fl/m} (11) Am Impact a3 
a— a 
Considering the airplane velocity to be constant, the , 
equations are finally arranged in the following dimen- °° " ta ” bin ” sa 
sionless form he TR 00 Fa : 
siomiess Io : . . : : : hee 
Fic. 2. Static and impact load-deflection characteristics of tir 
1 Gye) F, ‘a ”) L - a 
= - _— - _ (12) 
g mg mg M Mg L— cas 7 | 
- a . 7 \ a." 
Ye X4 “20 G ( Vo ) F 8 , rx f ) 
= — : l (13) | \ a 
g g mg mg : 
t 4 
It can easily be verified that Eqs. (12) and (13) apply 
to the landing condition on a level surface if we make = = 
zc" zero and apply the proper initial conditions. = 7 CYLINDER 
It was of particular interest in this investigation to | — de PISTON 
study the response of the system to excitation pro- Ory “i 
duced by a periodic ground function. Therefore, the 4 
F . Fane oe : ‘ : AIR iS NA 
simplest periodic function, the sinusoidal wave, was f ; ry The a 
f | : 
chosen. - i 4 while 
4 4 ’ 
SNUBBER ORIFICES 4 j Wem 
Zc = (a/2) [1 — cos(2rx¢/X)] (14) HN E | : 
SNUBBER VALVE i ia “y ANTI 
° ° ° ° — Zn id t | WW 
This may be written as a function of time, and the Hi la aN ||| [ 
. . ‘ PhP ot | HT . 
dimensionless acceleration term may be expressed as | We To de 
WA may | 
L720" z= (aw? 2¢) cos wt (15) ORIFICE | 
HTT 
It can-be shown that the quantity aw*/2g is propor- METERING FIN tH Hl | ies 
tional to the ground reaction that would be developed WN | | Vin 
. ae Wa a 
if the system were rigid. WT TA | 
- WA exten 
HTT A 
THE SHOCK-ABSORBING SYSTEM Hal | e 
| + 
x . . - Wet tl || } . 
I'he shock-absorbing system consists of a conven- tl 1 HH Subst 
tional air-oil shock strut, such as used on most military { | tH | 
and commercial airplanes, and a low-volume high load- 
‘ ; ‘ > ; ‘ ee Fic. 3. Shock-strut design. 
capacity tire designed for airplanes having limited ; ' —e The ¢ 
space for landing-gear retraction. The static and im- and t! 
pact load-deflection characteristics are shown in Fig. During the compression stroke, the fluid pressure be- tion. 
2. Fig. 3 shows the shock-strut design with the work- Jow the orifice is equal to the pressure p, and the pres- tempe 
ing parts essential to this discussion. The main orifice sure drop across the orifice Ap,. tion t 
is metered as shown so that its area is a function of sion f 
. “4: ‘ - ; . Pr = Pa + Ap, (16 
piston position. Syubber orifices are provided for re- : us to 
bound control and are effective only during the ex- The compressive force on the piston is found by mult tion ¢ 
pansion stroke. plying these pressures by corresponding piston areas. the st 
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Fic. 4. Tire load-deflection function, G(y-). 


fF, = p,(@, — ay) + pty — P(A, — Gy) 


= paly + Ap, ay (17) 


The air pressure p, is a function of piston position only 
while Ap, is a function of both position and velocity. 
We may write 


F, - F(y,) + R(y1) 51? (18) 


To determine F(y,), it is assumed that the air pressure 
may be found from 


PaVa’ = constant (19) 


where y is a suitably chosen constant. We may write 
/, in terms of the air volume when the piston is fully 


extended. 


V.= V, (20) 


— Apy\1 
Substituting in Eq. (19), the force F(3;) is obtained. 


F(y,) = constant X [a,/(V, — ap¥1)7] (21) 


The constant depends upon the position of the piston 
and the force applied to it at the beginning of the mo- 
Assuming that the air is initially at ambient 
temperature, we may consider the initia] piston posi- 
tion to have been reached by an isothermal compres- 
sion from the fully extended position. This permits 
us to express F(y;) in terms of the initial piston posi- 
tion and the force F, which holds the piston against 


the stops when extended. 


tion, 
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F(y:) = (1 — Ryn)?" [F,/( — Ky) (22) 


To determine the orifice function R(y,), we use the 
equation for turbulent flow through orifices or other 
restrictions. 


Q = CyAoV 2gh! (23)8 
This equation may be used to give the pressure drop 
across the orifice in terms of orifice area and piston 
velocity. 


Ap, = (3.9 XK 10~*)/cqa?) (a;/Ao) i? (24) 


From this, R(y:) can be determined. For the com- 
pression stroke it is given by the equation 


) 


(3.9 X 10-5/cq?) (a/A,2) (25) 


R(y, )= 


During the expansion stroke, the flow conditions are 
more complex because of the behavior of the return 
flow through the main orifice and flow through the 
snubber orifices. As the piston moves outward, oil in 
the chamber above the orifice is forced by air pressure 
to flow downward through the main orifice. The 
pressure conditions are 


Pr = Pa — AP 
Ps = Pa + Ap. 


Under these conditions the force on the piston is given 


by 


F, = pala, — as) + pps — ps(a, — Ay) 


= pf, — Apay — APp.(G, — G,) 26) 
and the expression for R(y;) is 
3.9 X 10 af 3.9 X 10-5 (a, — a,)* 
R(y1) = = - = r 7 
Ca” A 0” Caso” A so. 
(27) 


If the outward velocity of the piston becomes great 
enough, a critical velocity would be reached at which 
the resistance to flow through the main orifice would 
become equal to the pressure causing it. This con- 
dition would be characterized by the equation 

Ap, = Pa (28) 
At this velocity and beyond, the pressure py would 
drop to zero, and Eq. (26) would become 


F, = pPe(ap — ay) — Ap,,(a, — a) (29) 
In this event we choose to continue to use Eq. (18), 
with suitable modification of R, which now becomes a 
function of velocity as well as displacement. We 
obtain it by letting Ap, = p, in Eq. (26) and substi- 
tuting p, from Eq. (22). 


pa ~% (1 — kp)" F, _ 3.9 X 10- (a, — a,)* 
a, (1 — kp)? jn? bag? A,,’ 
(30) 
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Because Eqs. (12), (13), and (18), which make up 
the system of equations to be solved, are nonlinear, 
their solution is sought in numerical form. Therefore, 
it is necessary to introduce the numerical data to be 
used in the specific example to be considered. 

The tire function G(y2) is obtained from the tire 
manufacturer's data shown in Fig. 2. The impact 
characteristics are used in the analysis. However, the 
static characteristics are used to establish initial con- 
ditions; hence, a rigorous treatment would require the 
use of a different load curve for each initial tire deflec- 
tion considered. These, in turn, depend upon the 
weight of the airplane and the aerodynamic lift. Actu- 
ally, the differences involved are small and are ac- 
counted for by using two curves—one for the landing 
condition and one for the taxiing conditions. These 
curves are plotted in Fig. 4. 

‘At the present time, it appears that the constants 
y and cq can only be obtained by direct experiment. 
As a part of this general investigation, a number of 
landing-gear drop-test oscillograph records were an- 
alyzed in an attempt todetermine them. These records 
were obtained from drop tests on shock struts similar 
in size and design to the one considered in this prob- 
lem. Although many such drop tests have been made, 
not many include the direct measurement of air pres- 
throughout the compression and expansion 
Eight records in which air-pressure records 


sures 
strokes. 
were made were available for study, and the values of 
yy and cy, were based on results of the analysis of these 
records, Pressure-volume points from one of the 
records are plotted on logarithmic coordinates in Fig. 
5. It is seen that these points group fairly closely 
about the straight line representing y = 1.1, although 
in some of the records the scatter is somewhat greater 
than in the one shown. In all cases, the compression 
stroke begins with a y close to the adiabatic constant 
for air, decreasing in value as the compression stroke 
progresses and then dropping still further at the end of 
the stroke and during the expansion stroke. The loop 
thus formed probably represents the heat loss through 
the walls of the strut during the process. From the 
records studied, it is concluded that y is not a con- 
stant but that the heat transfer characteristics should 
be included in a rigorous study of such a system. How- 
ever, from an energy standpoint, the error involved 
in taking y as a constant is not great, and this simpli- 
fication of the problem seems to be justified. For the 
solution of the present problem, a value of 1.1 was 
chosen. 

The determination of cy from the records is more 
difficult, and the values found show considerable 
scatter when plotted. It was found convenient to use 
an orifice pressure coefficient, c,, defined by 


Ap, = ¢)V,2 (31) 
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test record. 
The relationship between the two coefficients is 
C = 3.9 X 10-*/c,4? (32 


The velocity of flow through the orifice is dependent 
upon piston velocity and net orifice area, both of which 
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Fic. 7. Orifice damping function R(y,) 
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Fic. 8. Metering pin and orifice 


may be determined from the piston displacement trace 
and metering-pin dimensions. 
V. = (dy A Vy (33) 


0 0 


The pressure drop may be found if we know the damp- 
ing force. In analyzing the records, this force was 
assumed to be the difference between the total com- 
pressive force on the piston and that due to the air pres- 
sure already determined. This neglects the effects of 
friction, which should be taken into account at this 
pont. Attempts to evaluate the frictional force were 


TABLE | 
Numerical Data 


Total stroke of shock strut 9.00 in. 
Static to fully compressed 1.66 in. 
Air piston diameter, d, 5.5 in. 
Air volumes: Compressed 19.67 — cu.in 
Static 59.02 cu.in 
Extended 233.49 cu.in 
Fluid piston diameter, d, 4.75 in 
Main orifice diameter 0.6875 in 


Rebound check: Inside diameter of cylinder, d, 75 in. 
Number of snubber orifices 


Diameter of snubber orifices }.0625 in. 


Specific gravity of fluid ). 834 

Airplane gross weight 16,460 Ibs. 
Weight on main wheels in statie position 40,580 Ibs. 
Weight of unsprung mass, per wheel 420 Ibs. 


not successful because the test data were not suffi- 
ciently complete for this purpose. However, near the 
end oi the compression stroke where the damping force 
approaches zero, friction was found to be unimportant, 
as indicated by the fact that the measured force due to 
air pressure was nearly equal, in all cases, to the total 
force. Computed pressure drops were plotted against 
orifice velocity. A typical plot is shown in Fig. 6. The 
points are widely scattered but fall between c, = 2 X 
10> and c, = 5 X 10. An average value of 3.5 X 
10~° is used in the majority of the differential analyzer 
solutions, although for comparison partial solutions of 
the drop-test condition were made using values ranging 
between 2.5 X 10°-° and 4.6 X 10~°. For reverse flow 
through the main orifice, cy is taken as 0.52, which is 
the value generally used for re-entrant tubes. The 
snubber orifices are small drilled holes with square 
edges, for which the coefficient is 0.61. The remain 
ing data required to determine /F(y,) and R(y) de- 
pend on the dimensions of the shock strut and are listed 
in Table 1. Curves of R(y;) as used on the differential 
analyzer input table are shown in Fig. 7 and represent 
the function as derived for the metering pin shown in 
Fig. 8, as well as for pins having slightly smaller diam 
eters at the large end. The latter curves were included 
to study the effect of metering-pin taper. 


DIFFERENTIAL ANALYZER SOLUTIONS 


Inasmuch as several excellent papers*~ are avail 
able which describe machines of the type used in the 
solution of the present problem, no further description 
is included here. A schematic diagram of the machine 
circuit is shown in Fig. 9. Equipment used included 
eleven integrators; two input tables for introducing the 
functions G(y2) and R(y,); and three output tables on 
which were plotted y,, yo, 31, 32, and F,/ Mg as func 
tions of time. The function F(y,;) was generated in 
the machine by using two integrators—one with in- 
verted input and output connections for generating the 
logarithm of (1 — k,y1). 

In order to begin the solutions on the differential 
analyzer, the initial conditions to be imposed must be 
known. These conditions determine the initial dis 
placements of the integrators and the input and output 
tables. In the landing conditions, the motion consists 
of two phases: the first beginning at the instant of con- 
tact between the tire and the ground and consisting of 
tire deflection only. The shock-strut piston does not 
move relative to the cylinder during this phase because 
the initial force F, is greater than the external force 
applied to the strut. The solution in this case may be 
found by making minor modifications to the machine 
circuit such that ¥, = # = », = O during this period. 
The initial conditions, then, for the landing condi 
tions are those that exist at the instant of ground con- 
tact. 


h=n=0, »=0, te = V2h 
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Fic. 9. 


In the taxiing cases, the initial conditions were 
determined by considering the airplane to roll into the 
sinusoidal ground, surface from a smooth horizontal 
plane. The system was considered to be in static 
equilibrium initially, and the deflections y, and Ve 
were found as follows. The tire deflection was taken 
from the static load-deflection curve of Fig. 4 for the 
condition 


G(y2) = (M+ m)g — L 


The initial piston deflection was determined by con- 
sidering the air to have undergone an isothermal com- 
pression from the extended position. This permits 
expressing y; in terms of the equilibrium force on the 


piston which, in turn, is equal to Wg — L. Thus, we 
write 
F,/ Mg 
"= es 
k, 1 — (L/Mg) 
The remaining arbitrary initial values are j, = i. = 0. 


These four initial values, together with those for 4, and 
je which may be determined from Eqs. (12) and (13), 
are sufficient to set all the initial displacements in the 
machine. The solutions are obtained by setting the 
machine in motion and by allowing it to continue until 
the desired information is obtained. In the drop-test 
solutions, curves were plotted through the compression 
and subsequent expansion strokes to the point at which 
the tire left the ground on rebound. Beyond this 


point, the equations are no longer valid. In the taxi- 


Schematic diagram of differential analyzer circuit 


ing solutions, the objective was to obtain steady-state 
oscillations of Therefore, the solutions 
were continued until the transients were damped out 
and the motion during one oscillation was reproduced 


or nearly reproduced, in the next. 


the system. 


RESULTS 


The drop-test solutions were carried out in anticipa 
tion of an experimental drop-test program on the same 
system. Accordingly, the data used in the solutions 
were chosen to agree with the anticipated experimental 
conditions. It is unfortunate for the purpose of com- 
parison that the experimental tests, which were con- 
ducted nearly a year after the numerical solutions were 
obtained, were carried out under somewhat different 
conditions so that exact comparisons cannot be made 
However, some conclusions can be drawn which indi- 
cate that the analytical approach using the equations 
derived in this paper can be used with good results. 

Differential analyzer solutions were obtained fot 
the metering pin shown in Fig. 8, as well as for the van 
ations in this pin indicated on the R(y;) curves of Fig. 
7. Curves for a representative solution are shown 1 
Fig. 10. In addition, partial solutions were made for 
a slightly different pin, not shown, which was subse- 
quently used in the experimental tests and to which we 
shall refer as the ‘‘test pin.’’ Solutions for this pin 
were not completed because they indicated that, for 


the gross weight considered and for the 27-in. free 
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drop, the tire would completely bottom, producing a 
peak load in excess of that permitted in the design speci 
feation. However, the ‘“‘test pin’’ was used in the 
experimental tests with a lower gross weight, which 
permitted it to be used without causing the tire to 
bottom. Our comparison, then, must be made for 
diferent gross weights and only over that interval of 
time for which the differential solution was carried out. 
Curves A and B of Fig. 11 show, respectively, the ex 
perimental and analytical results for the ‘‘test pin’’ and 
show close agreement up to the point at which the tire 
bottoms. The weight difference does not nullify the 
comparison, inasmuch as experiments show that weight 
is a secondary factor in determining the load-factor 
curve provided given static tire and piston deflections 
are maintained. Curve C in the same figure is the 
analytical solution for the metering pin shown in Fig. 
§. This pin is 0.1 in. karger in diameter at the end, 
which causes the load to increase more rapidly at the 
beginning of the stroke and gives the higher efliciency 
necessary to absorb the energy without bottoming the 
system. This curve is included to show the corre 
spondence with the experimental curve during the ex 
pansion stroke. Reasonable agreement may be ex 
pected here, inasmuch as the expansion stroke is con 
trolled largely by the snubber orifices, which were the 
same for both curves. In Fig. 12, the same three 
curves are shown plotted against piston position. 

The primary objective in carrying out the solutions 
for the taxiing case was to determine the maximum 
response amplitudes corresponding to a state of reso- 
nance. Calculations carried out for a linearized system 
in which damping was neglected indicated that the 
lower natural frequency was near w = 10 rad. per sec.; 
therefore, a frequency range of 5 S wS 20 was investi 
gated. It was found that a value of 0.25 for aw?/2¢ 
represented the maximum excitation that could be 
permitted without bottoming the system or causing the 
tire to rebound. The response curves for this value of 
aw*/2g plotted over the above frequency range are 
shown in Fig. 13. A maximum response at w» = 10 
rad. per sec. is clearly shown. A representative solu 
tion as plotted by the differential analyzer is shown in 
Fig. 14. The amplitude of the ground sinusoid corre- 
sponding with the maximum excitation at resonance 


may be calculated. 


2¢ faw" 2(3S86) m : 
a= — = (0.25) = 1.93 in. 
2¢ 100 


Ww 


Itis of some interest to compare the natural frequencies 
of the system as obtained from the differential analyzer 
solutions with those calculated using an idealized sys 
tem in which no damping is present and linear tire and 
shock-strut spring characteristics are assumed. In this 
linear system, the actual tire load-deflection curve is re 
placed by a straight line tangent to the curve at the 
point of equilibrium. The shock-strut function F(y:) 
is replaced by a linear function whose slope is obtained 
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by substituting the equilibrium value of y, in the first 
derivative of Eq. (22). The differential equations ob- 
tained for such a system are linear, and their solution 
yields, for this example, angular frequencies of 10 and 
157 rad. per sec. The lower frequency may be verified 
by measuring the period of oscillation from Fig. 14. 
In this mode, the motions of the piston and tire are 
essentially in phase. The higher frequency oscillation 
can be detected during the first 0.2 sec. of motion, par- 
ticularly on the velocity curves. The measured period 
is approximately 0.04 sec., which agrees well with the 
‘calculated frequency. It that in this 
the tire deflection and piston deflection are 180° out of 


is seen mode 


phase. 
CONCLUSIONS 


By means of the equations and methods given in this 
paper, the response of, and loads on, a conventional 
airplane landing gear in landing and taxiing may be 
found. 

Comparison of the differential analyzer solutions 
for the landing condition with experimental drop-test 
results show good correlation for the example treated. 
More extensive and exact comparisons should be made 
for further verification of the method and numerical 
data used. In this connection, the use of more com- 
plete instrumentation in drop-test programs to aid in 
the evaluation of strut air pressures, temperatures, and 
frictional forces would be extremely valuable. On the 
basis of present available data, values of 1.1 and 3.5 X 
10~*, respectively, were chosen for y and c,. 

The natural frequencies of the landing-gear system 
can be calculated with good accuracy by considering 
an idealized system with no damping and linear spring 
characteristics. 
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Analytical and Experimental Studies on 
Dynamic Loads in Airplane Structures 
During Landing 


T. H. H. PIAN* anno H. I. FLOMENHOFT Tt 
Massachusetts Inetitute of Technology 


SUMMARY 


Two assumptions are made in the analytical method for deter 
mining the transient stresses in airplane structures during land 
ing. The first assumption, that the landing impact forces are 
independent of the response of the structure, is justified by a 
simple analytical study. The assumption, that the 
lamping and aerodynamic foxces can be neglected, is verified by 
comparing the results of wind-on and wind-off drop tests of a 
The total transient stress is con 


second 


model wing in a wind tunnel. 
sidered to be the sum of the stresses that would result if the sys 
tem were restrained against vibrating and the additional stresses 
that result from the vibratory motion. This method hes been 
made more suitable fer practical application by the use of a tran- 
sient analyzer. The drop tests also provide an experimental 
check on the analytical method. 


INTRODUCTION 


[" DETERMINING THE TRANSIENT STRESSES in the air- 
plane structure during landing, it is desirable, for 
convenience in practical application, to make the sim- 
plifying assumptions that the landing impact force is 
independent of the response of the structure and that 
the damping and aerodynamic forces may be neglected. 
One purpose of the present paper is to summarize the 
analytical and experimental studies that justify the 
above two assumptions. 

A general approach of the analyses of transient 
stresses in elastic structures has been outlined in refer- 
ence 1. When the two assumptions described in the 
previous paragraph have been incorporated, the prob- 
lem of analyzing the dynamic landing stresses reduces 
to the determination of the response of an undamped 
elastic structure subjected to an external load of pre- 
xribed time history. It has been shown in reference |, 
that the response of an undamped structure can be rep- 
resented by a superposition of the responses of the 
iatural modes of vibration and that the response of 
each mode can be considered as that of an equivalent 
single degree of freedom system, excited by a prescribed 
generalized force. 

Two methods of approach in determining the tran- 
sient stresses have been given in reference 1. In the 


irst method, the total stresses are obtained by the 
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superposition of the stresses corresponding to the dis 
placements in the normal modes of vibration. In the 
second method of approach, separate calculations are 
made for the stresses that would result if the system 
were restrained against vibrating and for the additional 
stresses resulting from the vibratory motion. 


In applying the first method of approach, Biot and 
Bisplinghoff? have proposed the following method for 
estimating the dynamic stresses in an airplane structure 
during landing—adding, without regard to phase dif- 
ferences, the stresses corresponding to the maximum 
displacements in a few of the lower modes. The maxi- 
mum amplitude of stress in each mode is estimated 
from an envelope of ‘“‘dynamic response factors,”’ deter- 
mined from a statistical survey, which indicates the 
response to any force-time history that may be ex- 


pected in the landing. 


Drop tests of an airplane model at the National 
Bureau of Standards* were conducted to verify experi- 
mentally the method of approach of reference 2. The 
results indicate that the use of the envelope of dynamic 
response factors may be too conservative for design 
purposes. The experimental results also show that the 
neglect of phase relations and higher modes will cause 
negligible error in the determination of bending mo- 
However, this may not be 
For example, in the evalua- 


ment near the wing root. 
justified in a general case. 
tion of the inertia forces at the wing tip during landing, 
both of these factors must be included in order to ob- 


tain an accurate solution. 


The method proposed in the present paper for the 
transient landing stress analysis has been based on the 
second method of approach, in which the effects of 
static application of external forces and those due to 
vibratory inertia forces have been separated. Although 
the same method of approach has been followed by 
others,‘~® it has been made more suitable for practical 
application by the use of a mechanical-analogy-type 
analyzer, which has been described in a separate paper.’ 
3y means of this analyzer, the static stresses and the 
vibratory stresses due to the first three modes can be 
superposed, with proper regard taken of their phase 
relations. The accuracy of this method has been 
checked by drop tests of a model wing that was designed 
to simulate an airplane with overhanging engines. 
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Fic. 1. Simple two-mass system representing an elastic wing 
with landing gear. 
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This paper presents the summary of a research pro- 
gram consisting of analytical and experimental studies 
of dynamic loads on airplanes during landing. This 
work was part of a project, sponsored by the Bureau 
of Aeronautics, Department of the Navy, on the in- 
vestigation of stresses in aircraft structures under 


dynamic loading. 


(I) CoupLING EFFECT OF THE FLEXIBILITY OF THE 
STRUCTURE ON LANDING IMPACT FORCE 


It has been pointed out that the problem of deter- 
mining dynamic stresses in an airplane structure dur- 
ing landing can be simplified by assuming the landing 
impact force to be independent of the airplane struc- 
tural flexibility. The validity of this assumption is 
examined by the analytical study that follows. 

To simplify the problem, the following two assump- 
tions are introduced: (1) The landing gear struts are 
massless springs:and can be treated as a single spring 
of an equivalent stiffness A; (2) during impact, the 
structure assumes the same deflected position as its 
first symmetrical mode of vibration. 

The deformation of the structure is defined by the 
normal function ¢ of the first symmetrical mode, which 
is so normalized that the generalized mass of the sys- 
tem is equal to the total mass, .\/, of the airplane. The 
mode coefficient corresponding to the point where the 
landing gear is attached is ¢,. The motion of the air- 
plane is completely defined by specifying the displace- 
ment of the vibration mode, qi, and the position of the 
center of gravity of the structure, g. Thus, under 
these assumptions, the airplane structure can be con- 
sidered as a two degree of freedom system represented 
by a simple two-mass system, as shown in Fig. 1. The 
masses my and m, are connected by a spring of stiff- 
ness k, and a spring of stiffness A is attached to me. 

The motion of the equivalent system is defined by 
qo’, the position of the center of gravity of the system, 
and by Yogi’ and Yiqi’, the displacements of the masses 
my and m, respectively, from the c.g. axis. It can be 
seen that ¥y = —(mo/m1)o. 

The conditions necessary for the two-mass system to 
be equivalent to the actual airplane structure are as 
follows: 

(1) The total masses in both systems are equal 


i.e., 


My +m, = M 


(2) The generalized masses corresponding to 4h 


vibration modes in both systems are equal —i.e,, 
Wo?mo + hi2m, = M | 
or 
Yo?mo [1 + (729/1,)] = wi 


(3) The mode coefficient corresponding to the mass 
my for the two-mass system should be the same as that 
corresponding to the point on the actual structyp 
where the landing gear is attached—i.e., 


Z _ dy) 


(4) The natural frequency of the two-mass systey 
must be the same as the frequency of the first symmet 
rical mode of airplane vibration. The natural fr 
quency of the free two-mass system can be expressed as 


wo = V (mo + m)k/mom 
Thus, it can be written that 
V (my + m)Rk/mom, = a 


where w is the natural frequency of the first symmet 
rical mode of the actual structure. 
‘ Combining Eqs. (1), (2), and (3), we obtain 


m,/mMo = $77 


From Eqs. (1), (4), and (5), we can see that, if J/, 
¢,, and w, of the actual structure are given, 7, m, and 
k can be determined. 

In order to determine the effect of the flexibility oi 
the structure on the impact load, we consider two sys 
tems as shown in Figs. 2a, and 2b. Fig. 2a represents 
a rigid body; Fig. 2b, a two-mass elastic system. Both 
of these systems are assumed to be descending at con 
stant velocity, v, when the “landing spring’ touches 
the ground. It is desired to find the time variation oj 
the ground reaction F(t) in each case. By using % t 





Lm) 
mM. ro | 
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F(t) 


(a) (b) 


Fic. 2. Simplified systems representing rigid and elastic airplanes 
at landing 
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AIRPLANE 


LOADS IN 


indicate the displacement of the mass mp, it can be 


ven that, at the initial condition when the spring begins 


to touch the ground, x0(0) = O and X)(0) = 7. 
(a) Rigid Body 
The differential equation of motion is 


(mo + m)Xo + Kx = O (6) 


Its solution 1s 


Xo = WW wz) sin wl (7) 


where a, = © A/(mo + my) equals the natural fre- 
quency of the landing-gear strut. 
The time-history of the ground reaction can be 


written as 


‘ 


The solutions of these equations are 


o(1 + w) ‘| 3° 


STRUCTU 


| - | B? 
— | SIN Qew,t — = 
Qe (1 + pw) a> 
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F(t) = Kxo 

or (S) 
F(t) = (vK/wz) eid 


(b) Two-Mass Elastic System 


The motion of the mass mm is represented by the 
relative displacement x, with respect to the mass mp. 
The differential equations of motion are 


(mo + my) Xo + mx + Kxo = O 
m (Xo + X41) + ku = 0. 


The initial conditions are: at/ = 0, 


xo = 0, Xo = v 
x =0, 2; 


 . { 
— 1 sin cua 
a) 


xo = 
wy, (ai? — a”) (La + pw) a” 
(9) 
v(1 + pw) pe Bn 
a= - SIN Qew,! — Sin awl 
wy ( a” — ao) \ae ay 
where 
. = my,/ Mo 
free-free frequency of the two-mass system — ay k (mo + mm) K 
p = : : ‘ = 
frequency of the landing gear strut Wr mom my + my 
a = Vil +e+ 6)/2)4+ V[a + u t+ 62)/2]? — B 
wo = Vil +44 8)/2]) - Via +44 B)/2])? — B 


The time-history of the ground reaction for the case of el: 


; vk ( l+u ‘I 8° 
Fil) = - 2) = 
@, \aiy — a2” (1 + p)as- 


We wish to compare the maximum value of impact 
force for the case of the elastic body with that for the 
case of the rigid body. By examining Eqs. (8) and (10), 
we conclude that the ratio between the maximum im- 
pact force for the elastic body and that for the rigid 
body is equal to the maximum value of the expression 
within the parentheses of Eq. (10). The impact force 
ratios for various values of » and 6 have been calcu- 
lated and plotted as in Fig. 3. 


The result can also be represented in terms of the 


following two parameters: 


T; 
(1) Tr = duration of impact, natural period of the 
: considered mode of airplane vibration 
(2) @, = the mode coefficient in the considered 


mode at the point where the landing 


gear is attached. (The mode has 
been so normalized that the gener- 


EOS B? 
= | SIN aw, — = 
as (1 + mw) ae 


istic body is, then, 


(10) 


b . \ 
— | sin arr! ¢ 
ay 


alized mass is the same as the mass 

of the airplane) 
The ratio of 7,;/7Ty is equal to ‘'/28, and the mode 
coefficient ¢, is equal to Vu ‘Fig. 3 also indicates the 
variation of the ratio of impact force for a flexible 
airplane to that for a rigid airplane for different ratios 
of T,/Ty and different values of ¢,. 

This figure shows that, for 0.25 < 77/Ty < 2.5 and 
0 < ¢, < 1, the impact force for the flexible structure 
is 0.71 to one times that for the rigid structure. For 
current design of large landplanes with 7)/7Ty ~ | 
and 0 < ¢, < 0.5, the impact force for the flexible struc- 
ture would show a reduction of less than 10 per cent as 
compared with that for a rigid structure. 

Both experimental and analytical investigations of 
the effect of wing flexibility on the force-time relation 
in the shock strut during landing impact have been 
conducted at the National Bureau of Standards.* In 
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forces. During impact, the landing forces, assyne; It 
to be predetermined, are applied directly to the stry trams 
ture. It is desirable to determine the stresses Corre obtai 
sponding to the transient motion and_ to Superpog ing fi 
those stresses upon those corresponding to the initiy syste 
steady condition. ordin 
| We consider a general case in which the airplane js rigid 
g 60/- a | | subjected to a number of impact loads /\(f), F(t), . Fil) 
a | | F(t). These may be, for example, the normal apq rigid 
ud | drag loads at the landing gears at right- and left-hanq rollin 
g 40 i 2. = | | = sides of the airplane. The airplane is assumed to hay portt 
: | | the rigid body motions of the vertical, pitching, ang J 
= | rolling degrees of freedom. The rigid body rotations accel 
= 20 L | are assumed to be small. trast 
| The deformation of the structure is represented by : ° 
a a : | — the superposition of the displacements in the normal it 
@) 5 10 (5 20 25 4 lI . _ 
Ty modes of vibration of the unrestrained airplane. The 
0 4 5 transient displacement of each of the vibration modes 
FREQUENCY RATIO, say the rth mode, is indicated by g(t). The corte wher 
N.B.S. Analytical Results N.B.S Test Results sponding mode shape is described by a deflection func - 
x p70 & ?, = 445 tion ¢, which is so normalized that the generalized part 
+ ?,=-316 = ¢, =519 mass, ./,, is equal to the total mass, \/, of the air 
oS D=447 pei plane. The generalized mass is defined as 
© $=577 a $, =590 geen: 
oO p= | M, = J |e]? dm a whe 
Fic. 3. Effect of frequency and mass distribution on magnitude — where the integration is carried over the entire air 
of impact force 
plane. 
these investigations a simple two degree of freedom If the damping and aerodynamic forces are neglected, 
dynamic model, havirig a landing gear with both elastic the normal modes are uncoupled. Consequently, the \ 
and damping characteristics, were used. Since the equation of motion for the rth vibration mode is + 
analysis of the problem becomes much more com- Vv Meo. ) ‘ 7 
plicated, only a few cases have been considered. The et ae = Oi “= 70 
result of the test program indicated that for a duration where w, is the rth circular frequency of natural oscilla al 
of impact greater than one and one-half times the — tion and Q,(t) is the generalized force corresponding t i 
natural period of the wing, the force-time relation in the — the rth generalized coordinate. In the present case, “ie 
shock structure was substantially the same as though . es 
the flexible structure had been replaced by a rigid Ot) = SS dp? Fit) 13 if 
body having the same net weight. The peak force ; i” 
for 1.5 < 7T,/Ty < 2.5 showed a reduction of up to 10 where @,;‘” is the normalized mode displacement, in the st 
per cent due to flexibility. The analysis conducted same direction as the external load F,, corresponding V 
at the National Bureau of Standards showed that, for _ to the point of the structure at which F; is applied. I; 
0.231 < T7/Ty < 2.47 and for I > $, > 0, the impact Since the problem is linear, the responses in the vibra- ‘ 
force for the flexible wing was from 0.775 to one times tion modes of the airplane and, hence, the transient 
that lor the rigid wing. These results are also shown tresses in the structure can be obtained by considering 5 
in Fig. 3. _ each external load separately. Corresponding to the 
In conclusion, it may be stated that the effect Ol toad F,(t), the differential equation of motion in the 
wing flexibility is to cause such a slight reduction in the 44, mode is 
shock-strut force that it may be neglected in the de- 
sign of conventional airplanes. Mg, + Me,*q, = bri Filt) a a 
a ee The displacement q: and the acceleration g, can then . 
ledaieieent Cnnciieees heaes wis Sinaumaeen adem determined by applying the method given in reference ' 
AERODYNAMIC ForRCES ARE NEGLECTED or," ea ; 4 
q, = / F(t) sin w(t — t)dr (lo 
In the present treatment we assume that the airplane Me, Jo \ 
is making a landing approach with a constant sinking " 4 
speed, so that prior to landing it is in a state of equilib- -= dri" / “~ a q Fir)dr (16 
rium under the gravity forces and the aerodynamic MW So dr 
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IN AIRPLANE 


LOADS 


It is desired to determine the stress due to the 
transient motion by separating the part that would be 
obtained if the structure were restrained against vibrat- 
ing from the stress due to the vibrating motion of the 
system. The first part is the stress that the engineer 
ordinarily computes when the structure is assumed 
rigid. This stress is a result of the external force 
f(t) and the reaction forces that are caused by the 
rigid body accelerations in the vertical, pitching, and 
rolling degrees of freedom; hence, this stress is pro 
portional to F(t). The second part of the transient 
stress is due to the inertia forces corresponding to the 
the vibration The 
transient corresponding to the acceleration 


in the rth mode is proportional to q,(¢). 


accelerations in various modes. 
stress 
For conven- 


ence in computation, we express F(t) as 


= (Fi) max. Silt) (17) 


“ 


F(t) 


where f(t) is a nondimensional function with a maxt- 
mum amplitude of unity. We see that the stress at a 


particular point P can be written in the following form: 


op = >, Bafidt) + > BiKvlt) (18) 
where 
; " d 
K,(t) = cos a,(é — 7) filr)dr (19) 
0 dr 


As an illustration, we consider an airplane in the 
two-wheel symmetrical landing condition (Fig. 4). 
[he problem is one of determining the time-histories of 
the bending moments, torsional moments, and shears 
along the wing structure due to a vertical landing reac- 
tion '/>F(t) at each landing gear. We assume that the 
eastic axis of the wing is a straight line perpendicular 
to the plane of symmetry of the airplane and that the 
wing section is rigid in the chordwise direction. We 
employ the following notation and sign convention: 


V4 = total mass of the airplane 

I; = pitching moment of inertia of the entire airplane 

m(x) = mass of wing per unit wing span 

I(x = moment of inertia about the elastic axis per unit 
wing span 

S(x = unbalanced moment about the elastic axis per unit 
wing span, positive when c.g. is behind the elastic 
axis 

h(x) = normal function describing the mode shape of bend- 
ing of the elastic axis on the rth symmetrical 
mode 

a(x) = normal function describing the mode shape of twist- 
ing about the elastic axis in the rth symmetrical 
mode 

MV, = generalized mass in the rth mode, defined by 

6/2 
M, = (mh? + Ia? + 2Sha)dx 
—b/2 
Mode function is normalized so that M, = M 
or’ = mode coefficient of the rth symmetrical mode corre- 


sponding to the point in the wing structure where 
the landing gear is located 
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es 
EA. 
lic. 4. Airplane under symmetrical vertical landing load 

e = moment arm of the vertical reaction about the c.g. 
of the airplane, positive when the load is behind 
the c.g. 

a = moment arm of the vertical reaction about the elastic 
axis of the wing, positive when the load is behind 
the elastic axis 

W = location of the landing gear from the centerline of 
the airplane 

S = rigid body vertical displacement, positive upward 

0 = rigid body pitching displacement, positive nose down 

V = shear distribution, positive downward 

B.M. = bending moment distribution, positive when bottom 
surface is in compression 

T.M. = torsional moment. distribution, stalling moment 


positive 
The rigid body acceleration of the airplane consists 
of a vertical translation, 2, and a pitching motion, @. 
These can be expressed in terms of the landing impact 


force F(t): 


2 = F(t)/M (20) 
6 = eF(t)/I (21) 

Again, for simplicity, we write 
F(t) N max. gf (t) (22) 


and 


bo 


, m d . 
K,(t) = Cos w,(t — 7) — f(r)dr (23) 
0 dr 


where #,,,,. is the landing load factor used in the con- 
ventional stress analysis. 
The time-histories of the stresses in the wing struc- 
ture can be written in the form of 
oO = 


Bof(t) + > B,K,(b) (24) 
l 


The coefficients By and B, can be expressed as follows: 
(1) Shear Distribution, V: 


ws eS 
Bo = nox m + F es, <> & 
asia eS 
Bo = Nmexf m+ F dx —-M\, x< x, 


b/2 
B, = Nmax2Or f [mh + Sa) dx 








770 JOURNAL OF THE 





Photograph of the drop test model wing. 


Fic. 5. 


(2) Bending Moment Distribution, B.M.: 


. : eMS 
Bo = Nyax 2 m + dxdx, xX > Xz 
J b/2 6/2 by 
? ; eMS 
Bo = taxZ / m + dxdx — 
b/2 6/2 P 


l 
= M (xr — »|, BS Ky 


be my 
Bo = #96," ] / [mh” + Sa] dxdx (26) 
J 6/2 JS b/2 


(3) Torsional Moment Distribution, T.M.: 


*b/2 /eMI 
Bo = N max.& I i 


"2 feMI ~ l 
By = wor? i + S}dx — = Ma\, «<x 


s)a, 2 > 


*b/2 
B, = Nmax.2r” | [Ta + Sh] dx (27) 
x 


It can be seen that the problem of determining the 
transient stresses in the airplane structure during land- 
ing is reduced to the solution of equations of the same 
form as that of Eq. (24). The solution may be ac- 
complished by numerical or analytical integration. 
However, both 
impractical for the designer. 
type analyzer, which may be used to determine the time- 


processes are tedious and are often 
A mechanical-analogy- 


histories of the stresses in an undamped elastic struc- 
ture, has been developed at the Massachusetts Insti- 
tute of Technology.’ The two main functions of the 
analyzer are: (1) to evaluate the integrals Aj,(f), 
K,(t), ..., ete.; (2) to superpose the various terms of 
Eq. (24) with proper regard taken of their phase rela- 
tions. The first function is accomplished by the use 
of torsional pendulum simulators; the second function, 
by the proper design of the electronic pickup and re- 
cording circuit. 

EXPERIMENTAL VERIFICATIONS OF METHODS 
OF ANALYSIS 


(IIT) 


Experimental verifications of Biot and 
hoff’s method of determining transient stresses in the 


Bispling- 


AERONAUTICAL 
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airplane structures have been made in the Nationa 
Bureau of Standards.* The experimental investiga. 
tion described in this paper was extended to (1) th 
verification of the method of analysis proposed in the 
previous part and (2) the determination of the effects oj 
the aerodynamic forces on the transient stresses in the 
airplane structure during landing. 


Description of the Model 


A model was constructed which is dynamically similar 
to a large twin-engined airplane. It consists essentially 
of an elastic wing with two overhanging engines, g 
rigid fuselage, and a single-strut landing gear mounted 
at the centerline. 

The wing is rectangular in plan form, with structural 
strength provided by two magnesium spars connected 
by Duralumin torque tubes. The airfoil shape is main- 
tained by wood covering, made in sections so that it 
would not contribute stiffness to the primary structure, 
Large concentrated masses attached to the wing struc- 

ure comprised the nacelle section. 

The shock strut consists of a cylinder containing a 
piston and piston rod, at the bottom of which dual model 
airplane wheels are attached. Movement of the piston 
deflects a coil spring and forces oil through an orifice. 
The shape of the impact force-time curve was varied by 
changing the grade of oil and the stiffness of spring and 
by using different combinations of rubber pads on which 
the wheels landed at impact. 

The wing and the shock strut were attached to a 
parallelogram linkage that was connected to a steel 
channel fixed to the floor and could thus undergo only 
translation during the drop. Fig. 5 illustrates the 
general arrangement of the model wing. 


Instrumentation 


Instrumentation for the test was required for meas 
uring the bending and torsional moments at a wing 
root station, 1!/, in. from the center of the model, 
the bending moment outboard of the engine section 
at a station 8*/, in. from the center, the acceleration at 
a point near the wing tip, the acceleration at the fuse- 
lage, and the impact force in the shock strut. 

All quantities except acceleration were measured by 
means of wire-resistance strain gages. Each pickup 
for bending or torsional moment consisted of four 
gages, one at the top and bottom of each wing spat 
For measuring bending moment, the gages on the top 
and on the bottom were connected in series to form a 
pair. For measuring torsional moments, the gage on 
the top of one spar was connected in series with the 
gage on the bottom of the other spar to form a pait. 
Each of these pairs of strain gages was connected to a 
Consolidated bridge balance channel. 

The load in the shock strut was determined by meas- 
uring the strain in a fixed-ended steel beam on which the 
wheels landed. A pair of strain gages were attached 
to the top and bottom surface near the root of the beam 
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connected to a Consolidated bridge 


The signals from the bridge circuits 


and were also 
balance channel. 
were aniplified and recorded by an oscillograph. 
Accelerations were measured with a vacuum tube 
“transducer’’ accelerometer developed at M.I.T. The 
signals from the accelerometers were also fed into the 


recording oscillograph. 


Description of Test 


An important requirement in dropping the model is 
that the wing be in a strain-free condition before release 
in order to minimize the vibrations excited by the sud- 
den removal of the supporting forces. Consequently, 
the model was supported by seven strings, connected 
to a dural beam by means of adjustable screws. The 
model was released through a single wire that held the 
dural bar. 

Drop tests were made in a wind tunnel for air speeds 
of 0, 20, 40, and 60 m.p.h. At each speed, drops were 


STRUCTURES 
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made starting at a height of | in., and increasing in 
'/,-in, steps until a load factor of approximately 5 was 
attained. 


Test Results 


Oscillograph records showing the results of the drop 
tests from different heights and at different air speeds 
obtained. Typical records are shown in Figs. 
6and7. It was found that the impact force-time curves 
corresponding to different values of maximum impact 
However, the curves corre- 


were 


were not exactly similar. 
sponding to the same values of maximum impact but 
to different air speeds were similar. This latter fact 
facilitated the investigation of the effects of aerody- 
namic forces on the transient stresses in the structure 
under the application of the same external impact 
force. 

The test results are presented in the form of plots of 
Maximum moment and maximum acceleration against 
maximum strut force at the four different air speeds, 
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shown in Figs. 8 and 9. These results indicate that 
the aerodynamic forces tend to raise the maximum 
bending moments in the wing structure when subjected 
to the same impact force. This is due mainly to the 
decrease in effective wing angle of attack during the 
landing impact. The decrease in angle of attack causes 
a decrease in lift force, which, in turn, corresponds to an 


increase in bending moment. 


The maximum percentage increases in bending mo- 
ments due to 40- and 60-m.p.h. wing speed are shown 
in Table 1. These two wind-tunnel tests are equivalent 
to full-scale landings of 107 and 160 m.p.h., respec- 
tively. In general, the test results indicate that, for 
the usual landing speeds of less than 100 m.p.h., the 
effect of aerodynamic forces on the transient bending 
and torsional stresses is small. It can also be con- 
cluded, from Fig. 9, that the aerodynamic foree has no 
noticeable effect on the accelerations at the wing root 
and the wing tips. 


TABLE 1 
Maximum Percentage Increase in Bending Moments Dye ; 
Air Speed 


Item 40M.P.H. 60MPH 
Root Bending Moment 9.9% 13.9% 
Out-board Bending Moment 11.3% 18.6% 


Comparison of the records for still-air drop tests ang 
for drop tests with wind on (Figs. 6 and 7, respec. 
tively) reveals the important fact that, although the 
aerodynamic forces have little effect on the maximyy 
transient stresses that occur at the first peak, they tend 
to damp out the oscillations rapidly. 


Analysis 


The following calculations on the model wing were 
made in order to test the practicability and the accy- 
racy of the methods outlined in Section (II) for deter. 
mining the transient stresses in airplane structures dur. 
ing landing. 

The computations were carried out for the case of 
2.5-in. drop at 0 m.p.h. air speed, giving a maximum 
impact force of 45.5 lbs. In order to obtain a better form 
for the comparison of experimental and analytical re- 
sults, the results are presented in terms of nondimen- 
sional quantities. The calculated maximum moments 
and maximum accelerations corresponding to the con- 
dition that the structure is restrained against vibra- 
tory motion were chosen as the reference quantities, 
The time-histories of bending and torsional moments 
and accelerations were expressed in terms of the quan- 
tities (B.M.)/(B.M.)rep, (T.M.)/(T.M. rer, and (Ace. 
(Acc.) rep, respectively. 

For simplicity in analy~'s, the wing is replaced by a 
dynamic model as shown in Fig. 10. The wing 1s 
divided into five sections. The mass of each of these 


sections was considered to be concentrated in the span- 


wise direction at its center of gravity except for the 
center-section mass, which is included with that of the 


fuselage. In the chordwise direction, each section 
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TABLE 2 


Comparison of Results of Tests and Analysis 


774 JOURNAL OF THE AERONAUTICAL 
— Positive Maximum —— 
Item Measured Computed 
(B.M.)r/(B.M.)z,,, 1.36 1.23 
(B.M.)o/(B.M.)o,,, 2.50 2.30 
(T.M.)/(T.M.) ret. 1.98 1.605 
(Acc.)tip/(Acc.) ret. 3.40 3.80 
(Acc.) tus. /(AcC.) ret. 0.908 0.892 


The 
masses in each pair are made equal in the three outboard 
sections and are located at equal distances from the 
chordwise c.g. location. In the engine section the 
masses are distributed so that the location of the front 
mass is coincided with that of the center of gravity of 
the lead weight that was used for simulating the over- 
hanging engine. The front mass of the engine section 
is divided into two elastically connected masses, m, and 
m,, because the support of the lead weight is flexible. 
The first three coupled mode shapes and frequencies 
for the free-free wing with the fuselage restrained 


is replaced by a pair of concentrated masses. 


against pitching motion are computed by a process of 
matrix The flexibility influence coeffi- 
cients of the dynamic model, required in this compu- 


iteration.® 


tation, were experimentally determined. 
The transient stress in the airplane structure is given 


by 


o(t) = Bof(t) + BiKil(t) + BoK2(t) + B3K3(t) (28) 


where f(¢) is the forcing function in nondimensional 
form and A,(t) is the time integral defined previously. 
The constants Bo, B,, B., and B; corresponding to the 
bending and torsional moments in the wing are obtained 
by formulas similar to Eqs. (26) and (27). In the pres- 
ent case, since the lumped masses are used, the in- 
tegral expression of Eqs. (26) and (27) are replaced by 
summations. For the determination of the accelera- 
tions at the point P in the airplane structure, the ex- 
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- — Negative Maximum — 
Per cent 


Per cent 
error Measured Computed error 
—9.5 —0.59 — 0.580 1.7 
—8 —2.00 —1.78 1] 
—18.9 Sis ; ; 
11.8 —3.17 —3.30 —4.1 
—1.8 cere 
pressions for By and B, are 
Bo = Nmax.8 
B, -_ Nmax.ZOr” gp”? (29 


where ¢p‘” is the mode coefficient corresponding to the 
point P. 

When the results are presented in terms of nop. 
dimensional quantities, the following equation results 


oS ahi ken + kw 4+ Oe 
Oref a Bo “ Bo a Bo sai r By wv 


(30 


The impact force-time curve used in the analysis 
was taken from the shock-strut load record for the 2.5. 
The curve is 
rough in the first part of the drop as a result of the 


in. drop at zero wing speed (Fig. 6). 


higher frequency oscillations in the shock strut and in 
In the 
analysis, a faired curve was used because rigid body 


the impact beam used for measuring the load. 
acceleration records indicated that this ‘“‘hash’’ was 
not transmitted through the shock strut. 

The time-histories of the bending moment at the 
wing root, the bending moment just outboard of the 
engine station, the torsional moment at the wing root, 
and the accelerations at the wing tip and at the fuselage 
were evaluated by means of the transient analyzer and 
are plotted in Figs. 11 through 15. 
the results obtained from the drop-test records are also 


For comparison, 
plotted on each of these figures. A tabular comparison 
of the measured and computed results is given in Table 
2. 

It can be seen that, except for the torsional moment, 
agreement between the test and analytical results 1s 
fairly good. The error in the computation of the tor- 
sional moment is believed to be the result of inaccuracy 
in the estimation of the basic torsional properties of the 
wing structure. 

Inspection of the curves of Figs. 11 through 15 indi- 
cates that the stress-time curve obtained from the 
analysis has the same general shape as that obtained in 
the tests. 
normal modes higher than the third is not significant. 


This is a good indication that the effect of 
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An Extension of the Theory of the 
One-Dimensional Shock-Wave Structure’ 


LEONARD MEYERHOFF?t 
Polytechnic Institute of Brooklyn 


SUMMARY 


This paper presents the method and results of an analytical 
investigation of the flow inside a shock wave in a steady one-di- 
mensional flow. Two categories are treated: The first considers 
the problem for constant gas properties, while the second is con- 
cerned with the influence of specific heat, conductivity, and vis- 
cosity, dependent on temperature. 

For the problem with constant gas properties, numerical solu- 
tions are given which show that there must be an infinite number 
of shock-wave type solutions. Only one of these solutions ter- 
minates in a region of uniform flow with the values of Rankine- 
Hugoniot. These partly by differential 
analyzer and partly by hand, are given for initial Prandtl Num- 
bers, P, = cpu/k, of 0.75 (T = 0°C.) and 0.657 (T = 500°C.) 
A comparison shows that the sbock-wave thickness for an initial 
Mach Number of 2.0 is about.three times greater for P; = 0.657 
than for P, = 0.75. 

For the problem with temperature-dependent gas properties, 
the method and results are shown in detail for the case of uniform 
exit flow. It was necessary to calculate the change in the limit 
values of Rankine-Hugoniot because of variable gas properties 
in order to perform the required numerical integrations. The 
temperature variations used for cy, u, and k correspond to experi- 


solutions, obtained 


mental data. A numerical example, given for an initial Mach 
Number of 4.0, shows that the shock-wave thickness is approxi- 
mately two and one-half times greater than that for constant gas 
properties; it is about two and one-half times greater than the 


mean free molecular path. 


INTRODUCTION 

TERM “STRUCTURE OF THE SHOCK WAVE”’ is 

used here to mean the continuous variation of the 
flow variables of velocity, pressure, temperature, den- 
sity, etc., inside that small region that is called a shock 
wave. Ashock-wave type of flow variation is character- 
ized by a rapid (and continuous) change in the flow quan- 
tities across a transition region (which is the shock wave), 
whose length, based on mathematical studies, f is known 
tobe some multiple of 10° cm. The greater the initial 
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Mach Number, the smaller the transition length. 
These analytical investigations also show that in this 
transition region the curve of any flow quantity ex- 
hibits an inflection point. In this paper, the inflec- 
tion point, which seems to be necessary analytically for 
the existence of a shock wave, is referred to as the center 
of the shock wave. It is known that the rate of change 
of the flow variables is greatest at the shock-wave 
center. For example, there the velocity and tempera- 
ture gradients are, in general, some multiples of —10” 
(cm./sec.)/em. and +10* °C./cm., respectively. ** 

Because of the mathematical difficulties involved in 
using the Navier-Stokes equations, all known studies 
of this problem have been restricted to one-dimensional 
flow and, until recently,*'° with the assumption of 
constant gas properties. With these simplifications, the 
flow variation through the shock wave is fairly well 
established by these studies. 

One question that has always been raised is whether 
it is valid to use a continuum approach to such a phe- 
nomenon that occurs in so small a region. Several 
papers have been written,'* '* which treat the gas asa 
discontinuous medium by means of molecular theories. 
While these studies have not changed the known 
phenomenologic aspect of the problem, they have 
added important knowledge to questions about sta- 
tistical equilibrium, heat capacity lag, etc., of the gas. 
The results of one of the molecular studies'* have led 
the authors of reference12 to conclude, in their paper on 
the shock-wave thickness, that the Navier-Stokes 
equations are valid for this problem within certain ob- 
vious limits—that is, as long as disassociation and 
condensation effects are not present or that the gas is not 
highly rarefied or has a high stagnation temperature. 

In the present paper, the one-dimensional Navier- 
Stokes equations are applied for shock-wave flows to 
show the following: 

(a) The influence on the shock-wave solutions due 
to the initial temperature—i.e., due to initial Prandtl 
Number, c,u/k. 

(b) The relation of the classical shock-wave solu- 
tion, which terminates in the limit values of Rankine- 
Hugoniot, to other shock-wave type solutions con- 
tained in the Navier-Stokes equations. 


** For the example given in Section (IV) of this paper for 
temperature-dependent gas properties and an initial Mach Num- 
ber of 4.0, it is found that du/dx = —10'° (cm. )/em. and 
dT/dx = +3.64 & 108 °C. 


sec. 


em. 





776 JOURNAL OF THE AERONAUTICAL SCIENCES—DECEMBER, 1950 


(c) The method and results for temperature-vari- U = velocity ratio, u/uo Firs 
able viscosity, conductivity, and specific heat, of which I = temperature ratio, 7/7) ol iced 
age : ge = viscosity coefficient 
the variations follow the existing experimental data 4 Se ee Ea. (3 
for tl , teristi k = thermal conductivity q: \ 
ce OF ar: > 4 *S - 
or these gas Cnaracteristics. Cy = specific heat at constant volume ] 
Cp = specific heat at constant pressure (# 
SYMBOLS y = ratio of specific heats, cp/c; ax 
is ; P, = Prandtl Number, cpu/k 
x positional coordinate . : 
; R = gas constant expressed in physical system of units 
u = velocity : . ; , 
P Cp — Cy After 
p = mass density E hen 2 fe iT Sy 
b = pressure 4 = intrinsic energy, J ove piacec 
T = absolute temperature ( )'== d/dx yeloct 
a@ = speed of sound Q = subscript denoting a flow value atx = — ture 
m = mass flow, up fi = subscript denoting a flow value at the shock-wave center 2a) a 
M = Mach Number, u/a 1 = subscript denoting a flow value atx = +. 
(I) THe System OF EQUATIONS OF A STEADY ONE-DIMENSIONAL VISCOUS HEAT-CONDUCTING FLow From 
The system of equations defining the special case of steady one-dimensional flow is a special case of the general a 
system of the three-dimensional equations. These one-dimensional equations are given below for (a) constant _? 
physical gas properties and (b) variable physical gas properties, where the subscript 0 denotes initial flow condi- 
tions. Eq. | 
| form 
(a) Constant Gas Properties! | (b) Variable Gas Properties? 
' l ’ 
= ' ore b(u’ 
The equation of continutty is: 
(d/dx)(pu) = 0 (] 
The integral of Eq. (1) is 
wher 
u = pollo = Mm (1. 
p Porte la the v 
where 7 is an integration constant. 
The Navier-Stokes equation of motion is: 
9 | 
du dp 4 du du dp 4d du js 
m—-=— ny) (2) || Ne * - ) 
dx dz: 3° d* dx dx 3dx\ dx Ec 
reler 
Zach of these equations has a first integral with a second integration constant mC—namely, gas | 
; ge 
; 4 du ; 1 du the 
mu — mC = = of = Mo (2a) mu — mC = —- -+- = K (Ja rath 
3° dx 3. dx 
this 
In this paper mC is derived from initial uniform flow conditions, denoted by the subscript 0; therefore way 
. e of tl 
Po + muy = mC (2b, 3b : 
| deri 
If 
stan 
The equation of state for a perfect gas is: deri’ 
= orde 
p= RT p (4 
| cons 
The equation for the thermodynamic balance of heat flow is: whe 
4 a) s a?T dT p du 0 6 4 du\? Pa d put dT du 0 6) ’ 
— Mo = 0 Cow. = ») — — me, = — = (0) Do 
3 dx dx? ” dt dx 3° \dx dx \ dx dx dx 
—- fo 
The system of equations for variable gas properties where the velocity u is the dependent variable and the I 
7 . = 0 


can be combined to form a third-order, nonlinear ordi- abscissa x is the independent variable. This equation 
nary differential equation with variable coefficients, is derived as follows: 











of units 


€ center 


eneral 
istant 
‘ondi- 


ib 





ONE-DIMENSIONAL 


First, the pressure, p, in the heat flow Eq. (6) is re- 
placed by its value in terms of the velocity given by 
Eq. (3a), so that Eq. (6) changes into 


d } dT du d [{[u* 
k — mec, — mt +m = 0 

dv ( dx dx dx dx ( 2 ) 

(6a) 
After the temperature 7 and its gradient* T’ are re- 
placed in Eq. (6a), in terms of the velocity, the desired 
yelocity equation will follow; this required tempera- 
ture relation, obtained by a combination of Eqs. 
(3a) and (4), is of the following form: 


| (C + * ( | 
u — un) u? 
R 3m 


From Eq. (3b), the temperature gradient is 


(3b) 


) 2 x!’ 
M 9) =u \7 i) 
T’ = (u?)” + - (1) | (5c) 


> 


3m — 3m 


, ‘ O\ ve 
uC — (u?)"*+ 
R | 
Eq. (Ga), with this value of 7’, assumes the following 
form in terms of the velocity 1w—namely: 


9 
(b(u?)"|’ + [(b2 — mk)(u?)’| — — wcym(u)” + 
” 


[ : 


where the coefficients of Eq. (7) are variable and have 


WwW | bo 


meu’ | (u?)’ + [/u’]’ — qu’ =0 (7) 


the values shown below: 
b = (2/3) nk 
bs = (2/3) y’k, = mkC (a 
h = m*[(c + ¢,)/2], q= m*Cc,) 


~ 


Eq. (7) is treated by a perturbation procedure in 
reference 3 by Reissner and Meyerhoff, in which the 
gas properties y, c,, and & are expressed as functions of 
the temperature 7. Because the form of Eq. (7) is 
rather complicated, it was found preferable to treat 
this problem with variable gas properties in another 
This treatment is presented in a later section 
First, however, certain results will be 


way. 
of the paper. 
derived for the case of constant gas properties. 

If all the gas characteristics are assumed to be con- 
stant, then Eq. (7) reverts to the differential equation 
derived in reference +. This latter equation is a second- 
order, nonlinear ordinary differential equation with 


constant coefficients—namely, 


(u?)"by — (u*)’fo + u’lo + u*ho — ugo = D_ (8) 
where the constant coefficients here are 


2 s Yo + l 
by = = uoko, ho = MC ms 
7) 


(b) 


> ( 
r “ e { 
fo=m (i > con) go = m*Cc,, 

, 2] 
lo = mCkp, D, integration constant 


*’=d/dx 
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It can easily be shown that Eq. (8) is identical to those 
equations developed for this same problem in different 
forms by Rayleigh® and Taylor.® 

Only one exact closed form solution of Eq. (8) is 
known which exhibits a continuous flow variation through 
the shock wave. This solution, first developed by 
Becker,’ exists if the nondimensional ratio of the gas 
characteristics, P, = c,u/k, has the special value of 
0.75. The velocity versus x curve obtained from this 
exact solution begins and ends in regions of uniform 
flow, where, at the flow exit, the Rankine-Hugoniot 
limit values appear. 

These limit values are another exact result that can 
be derived from Eq. (8) (for the trivial solution u = 
constant). If this value of u is inserted into Eq. (8), 
the following quadratic relation appears: 


u*hy — ugo = D (8a) 


The result can be interpreted in two ways: either that 
there are two constant velocities possible all along the 
flow or, if the velocity is variable along the flow, they 
can appear only asymptotically atx = +o. This is 
because, at the location of these constant solutions, all 
derivatives become zero. The latter possibility leads 


to the Rankine-Hugoniot limit values—namely, 


ui*ho — mgo = D = uo*ho — Ugo 


uo 


(8b) 


where the subscript 1 denotes the limit value at the 
exit; also, 


Co = | + (1 yoMo?) 


and 

U; = 2[v0/(ve + 1)]Co — 1 
is the exit limit ratio value expressed in terms of the 
initial flow conditions. 

The derivation of the corresponding limit values 
for variable gas properties from the differential Eq. (7) 
is more complicated. The reason is that, contrary to 
Eq. (8), in Eq. (7) only derivatives of u are present. 
In Section (III) of this paper it is shown how these 
limit values for variable gas properties can be derived 
from the individual equations leading to Eq. (7). 

Before going on to the material in Sections (III) and 
(IV) derived for variable gas properties, some new re- 
sults will be presented in Section (II) which are ob- 
tained for the case of constant gas properties. 

It was decided to examine Eq. (8) for other possible 
shock-wave solutions besides that given by Becker.’ 

+ According to the data in references 8 and 9, the value of P; = 
0.75 would seem to appear at a temperature below 0°C. How- 
ever, within ‘the experimental accuracy, it may be considered as 


the value at 0°C. 
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It was further decided to determine the influence of an of the M.I.T. Differential Analyzer* and partly by han¢ 








initial Prandtl Number different from 0.75. In Sec- computation. 
tion (II), the procedure and the subsequent results of r 
numerical investigations are presented (for an initial * The suggestion to use the analyzer machine was made by 1 
Mach Number .\/) = 2.0) as obtained partly by the use Dr. Hans Reissner. 
t 
(II) THE NUMERICAL INVESTIGATION OF EQ. (8) (CONSTANT GAS PROPERTIES) 
It seemed advisable to rearrange the differential Eq. (8) in terms of the velocity ratio u/u = U. Eq. (8) they 
assumes the form: 
(U*)"b, — (U?)'fi + U' + Uy — Un + Uih = 0 (Se 
where 
U = u/Uuo ()’ =d/dx 
bo i f2 ho m vyt1 | 
hh = = = Moo h = = Ge 
Uo Uy \3 Uo Uo 2 | 
; fo m 2 qo m* _ 
fi = —_ ko of . Crohto n = _ = CiGx j (Cc | 
Uo Uo e Uo~ Uo | 
; } 
L, =To/uo? = (m/uo) Coko D, = D/u = —hU, 
. . ° , P Yo . the Rankine-Hugoniot limit ratio 
Co = c uy = l + (1 oMo") U; = 2 CO — . = | 
Yo+ 1 value 3. 
The Boundary Conditions one of these values is assumed, in addition to LU” = 0 be ( 
There are five conditions to be satisfied by the solu- the other value (and also all successive derivatives at 
tion of the differential Eq. (Sc). x = 0) can be determined from the differential Eq. 
* Fe oa: ve “he ‘ r cP ms ‘ _ "y . : 7 7 
_jThree entrance conditions are assumed to be located (Se). hat is, for U"” = 0 (at x = 0), Eq. (Sc) shows I 
° . ° ‘ * re 2c 7 _ iy > . ecte r > 
asin previous reports, at so great a distance from the that the values of U and l are connected by the Pra 
shock wave (expressed by x ~ — ©) that the flow in quadratic relation tior 
this region can be considered as uniform with prescribed 2m - Pe 
values of velocity, temperature, and pressure. These U, = in n( c Qn ity 
. Z —_ - 0 =—!}) ? 
three prescribed values are sufficient to determine the are 
ree integration constants appearing in Eq. (Sc)—viz. a Co \° P ' : bas 
three integ pp g , m?(U.— —-m#(U.— 1)(U.-— U)| ® 
2n1 are 
(a) mM = polo , 
, ‘ slig 
p | where the subscript c denotes the center (inflection uM 
: 0 ; : Ma 
(ib) Co = +u=1+ (Sd) point of the shock wave and 
: I 
m yoM 0? 
; D) Cpéo the Prandtl The 
(c) Dy, = MU, : m=1+ r..: P, = Sy : : 
3Y0 ky Number (d I 
Two more conditions must be chosen at some point (2/3) [(yo + 1 \P cur 
; ; 9 = (2/3) )/¥ 
along the flow in order to determine the flow downstream ™ ” ‘i cho 
of the shock wave. The only clearly definable point The choice of the velocity ratio U, = u,/uo at x = 0 ee 
is the inflection point (called the center of the shock must be interpreted, as stated above, as the equivalent The 
wave) on the transition curve across the shock wave of an exit condition for the flow downstream of the isd 
between the higher and lower velocity. This point, shock wave (usually denoted by x —~ +). since — 
essential for the existence of a shock wave, is chosen at the no rigorous procedure seems to exist to determine 4 aa 
origin of the coordinate system, so that one more general relation between UL’ at x —~ + © and U, (atx = the 
boundary condition is 0), a certain special guide has been used for the selec- I 
oY ( 9 tion of a value of U. where U,” = 0. This guide 1s, Fig 
= 0, atx = 0 (se) F : . on ry . 
( eed taken from the known closed form solution of Eq. (S¢ to 
The final condition to be chosen before the numerical given by Becker? for P,, = 0.75. From Becker's solu- fact 
solution is started, can be either LU’ or U’ atx = 0. If tion, a simple formula relating LU’, and (Ui), —~ 4. cam are 
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U, VELOCITY RATIO= (4/u,) 
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“1G. 1, e-dimensional shock-wave solutions for a Pranc¢ Number of 0.75. 
Fic. 1. One-dimensional shock lut f P itl Num £0 


be derived —namely, 
UU, = VU, (10) 


Two general cases, based on two values of the 
Prandtl Number, are chosen for the numerical solu- 
tion of Eq. (8c)—i.e., P,, = 0.75, To + O°C. and 
P,, = 0.657, Ty = 500°C.* Furthermore, three veloc- 
ity ratios, which are the same for the two general cases, 
are chosen at the inflection point. One velocity is 
based on relation (10) given above; two other values 
are also selected, one slightly larger and the other 
slightly smaller than that from Eq. (10). The initial 


Mach Number for all solutions is Wy) = 2.0. 


The Influence of the Boundary Conditions 


Three velocity curves are shown on Fig. 1. These 
curves, are based on a Prandtl Number of */, and three 
chosen velocity ratios U, at the inflection point 
namely, (a) 0.70000, (b) 0.61237, and (c) 0.60000. 
The velocity ratio (b) is computed from Eq. (10), which 
is derived from the exact particular integral of Becker;’ 
it will be recalled that this integral of Becker does not 
contain any nonuniform solutions at plus infinity for 
the velocity and the other variables. 

It will be noticed that none of the three curves on 
Fig. 2 shows the characteristic uniform velocity down- 
stream of the shock wave of Becker’s solution. This 
fact is not surprising for curves (a) and (c) because they 
are derived from boundary values at the inflection 


point for which exact solutions in closed form are not 
known and probably do not exist. However, the 
velocity curve (b) should have approached uniformity 
downstream, because curve (b) is Becker’s case since 
the boundary values used at the inflection point are 
calculated directly from Becker’s exact solution. 

When this surprising feature of curve (b) first ap- 
peared in the results found by the M.I.T. Differential 
Analyzer, the Computing Group initially considered 
machine errors as a possibility of its cause. However, 
after the work was investigated by an extremely accu- 
rate hand computational method, it was found that an 
instability was inherent to the problem and could not 
be avoided by any method of computation. The use 
of the term “‘instability’’ here means that there exists 
an extreme sensitivity against minute changes of bound- 
ary values (for example, at a deviation less than 10~‘ 
percent). 

It is natural to ask what causes the velocity curve 
(b) on Fig. 1 to agree with Becker’s solution up to a 
certain positive abscissa value and then become in- 
stable? The answer, as it will now be shown, is directly 
related to the accuracy or number of significant figures 
to which the velocity ratio U and velocity gradient U,’ 
for curve (b) were computed for the numerical analysis. 

The explanation is seen clearest from the more 
straightforward derivation of Becker’s solution pre- 
sented by Morduchow and Libby.'® One sees from an 
inspection of this solution that, if P,, = 0.75, then the 
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order of the nonlinear Eq. (8c) can be reduced to the 
following first-order equation: 


3m yt 1 


UU’ — U— U,) = 
S Ho Yo ' 
3 l - 
3m yo + 1 N exp. (= x) (11) 
4 uo Yo Ro 


If the value of the integration constant N in Eq. (11) 
is taken as zero, then the resulting differential equation 
leads to the particular solution of Becker. Only from 
this particular result does one obtain a flow with a uni- 
form velocity at plus infinity. For any other real 
value of the integration constant, N, the solutions of 
Eq. (11) will show nonuniform flows downstream of the 
inflection point; such solutions, moreover, cannot be 
obtained in a closed form in terms of elementary func- 
tions. 

Now, when the boundary values U, and U,’ (calcu- 
lated from Becker’s solution for N equal to zero) to be 
used by the differential analyzer to derive curve (b) 
on Fig. 1 are inserted into Eq. (11), the left side of 
Eq. (11) should reduce to zero. This did not occur 
exactly for the numerical solution of Becker's case re- 
ported here on Fig. 1. The reason, therefore, that even 
a highly accurate numerical solution [derived from 
Eq. (8)] will not, in general, reduce this left side of Eq. 
(11) to zero is the existence of an exponential, exact 
complementary integral with an exponent of extremely 
high value. The order of magnitude of this exponen- 
tial termis 3 X 10°. 

The conclusions from this discussion and from the 
results shown on Fig. 1 are that (a) there are an in- 
finite number of shock-wave type solutions possible 
from Eq. (8) and (b) that only one shock-wave solu- 
tion, however, terminates in a state of uniform flow. 
The solution for uniform exit flow is the one that yields, 
at the exit, the limit values of Rankine-Hugoniot. 
These limit values, located theoretically at plus in- 
finity, show, as is well known, that the limit pressure 
is greater than the pressure at the flow entrance. This 
flow, while predicted theoretically and satisfying the 
conservation conditions of mass and momentum, 
could not be produced experimentally. It is to be 
anticipated that, when a two-dimensional solution is ob- 
tained for this problem, a flow will appear, as shown by 
experiments, with finitely located limit values of higher 
pressure and lower velocity, for which, however, the 
final downstream pressure is lower than the entrance 
pressure. 

The discussion above is concerned with the numerical 
solutions derived for a flow with a Prandtl Number of 
0.75. Attention will now be called to the solutions 
obtained by numerical integrations for a Prandtl 
Number of 0.657. These latter solutions are plotted 
on Fig. 2. Here, the boundary velocity ratios, U,, at 
the inflection point are the same as those on Fig. 1. 
The initial Mach Number is still 2.0. 
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When the Prandtl Number is not 0.75, there is yo 
exact relation available to calculate a velocity ratio 
U, for which uniform flow appears asx — + ©. Hoy. 
ever, by the assumption of continuity of solutions 
for the case of P,, = 0.657, shown on Fig. 2, it can be 
expected that there will be a solution that becomes 
asymptotically uniform asx — +o. (In fact such g 
solution could be derived from the general procedure 
developed in the last section of this paper.) 

An inspection of Fig. 2 shows that these numerical] 
solutions are similar in form to those on Fig. 1—ie, 
nonuniform velocity variations downstream of the 
shock wave. A significant characteristic of the curyes 
for lower Prandtl Number, shown on Fig. 2, may be 
seen in the smaller velocity gradients across the shock 
wave.* 

The dotted line curves shown on Figs. 1 and 2 were 
calculated from the perturbation solution developed 
in reference 4. 

For the purpose of comparison, the Becker solution 
is plotted on Fig. 3, together with the corresponding 
solution for P,, = 0.657. This comparison is based 
on the same initial Mach Number and mass flow and 
the same (boundary) velocity at x = 0. The curve for 
P,, = 0.657, a replot of curve (b) on Fig. 2, before it 
rises downstream of the shock wave is probably a good 
approximation to the solution with uniform flow as 
i di ae 
gradient across the shock wave is apparent for the lower 
Prandtl Number. It should be noted that, whereas 
the curves on Fig. 3 are obtained for equal initial values 
of Mach Number (but unequal initial pressures), the 
curves on Figs. 1 and 2 are calculated for equal initial 
values of Mach Number and pressure (but unequal 
Because of the interdependence of initial 


From this comparison, the smaller velocity 


mass flows). 
temperature, Mach Number, and mass flow, it is not 
possible to make the comparison for exactly the same 
initial conditions. 

The approximate solution developed by Taylor’ is 
shown on Fig. 3 for comparison. 

This completes the section on the shock-wave struc- 
ture for constant gas properties. In the sections to 
follow, an analysis is given for the case of the tem- 
perature-dependency of the gas properties yu, k, and 


Gu 
(III) THe DETERMINATION OF RANKINE-HUGONIOT 
LIMIT VALUES FOR TEMPERATURE-DEPENDENT GAS 


PROPERTIES 


In Section (IV) of this paper a method to determine 
the shock-wave structure is developed for the realistic 
case of temperature-dependent gas properties. Needed 


* The velocity gradients at the shock-wave center (x = () are 
about one-third smaller for ?;, = 0.657. This means that the 
shock-wave thickness, as defined by Prandtl,'! is therefore three 
times greater. This thickness is defined as Au/u,’, where 
Au = difference between velocities at x = +© (i.e., uy — u 
and u,’ = velocity gradient at the shock-wave center. 
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U, VELOCITY RATIO u/U,) 
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| THE ANALYSIS 
| iS BASED ON A 
PRANDTL NUMBER OF 
0.657, INITIAL MACH 
NUMBER OF 2.0 AND A 
CHOICE OF DIFFERENT 
VELOCITIES AT THE INFLECTION 
POINT OF THE SHOCKWAVE, X=0 
| (SEE TABLE I) 


COEFFICIENTS OF VISCOSITY AND 
CONDUCTIVITY ARE ASSUMED 
CONSTANT FOR THIS ANALYSIS — 
(BASED ON INITIAL TEMPERATURE) 


_INITIAL VALUES AT X-~ -co 
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Fic. 2. One-dimensional shock-wave solutions for a Prandtl Number of 0.657, 
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Fic. 3. The effect of Prandtl] Number on the velocity gradient across a one dimensional shock wave in a viscous heat-conducting 
compressible fluid 
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for that work are the Rankine-Hugoniot limit values 
for temperature-dependent gas properties. These val- 
ues will now be derived. 

The differential Eq. (7) of velocity is shown to follow 
from the system of one-dimensional flow equations in 
which yu, k, and ¢, are assumed to vary along the flow 
with the temperature. A trivial solution of Eq. (7), 
it was shown, is u = constant. This trivial solution, 
does not lead to a quadratic equation for the velocity 
(the Rankine-Hugoniot limit values) similar to Eq. 
(Sa), since Eq. (7) contains terms that are only deriva- 
tives of the velocity. Therefore, in order to derive 
the analogy to the limit values of Rankine-Hugoniot 
for the temperature-dependency of all the gas proper- 
ties, the individual equations of the system must be 
applied directly. 

The heat balance Eq. (6a) is the first to be used. 
If use is made of the ideal gas relation for the intrinsic 
energy—namely, 


caT = dE (12) 
then it is possible to integrate Eq. (6a) once. The re- 
sult of this integration is 

dT ) J mu* 
k — mE — mCu + > = constant (13) 


ax “ 


We are interested in determining the relations between 
the flow variables that will be found in the uniform 
flow regions on both sides of the shock wave; in these 
regions the term k(d7'/dx) must vanish. Eq. (13) can 
then be written as follows: 

uy" “ Uo” 
a Cuy — > = constant 


(FE, = Eo) + Cu, ~~ 
(13a) 


+ <— ©, 


where the subscript | denotes the flow at x - 


Ti 
(EF, — Eo) = f eat 
Te 


Hence, Eq. (13a) can be written in the form 


wry ae i u\" = Uo” 
c,dT + Cu, — — = Cuy — : 
To ~ - 


and is a first relation between 7}, mu, and the initial con- 
ditions denoted by the subscript 0. 

A second required relation in terms of 7; and 1 is 
given by Eq. (3b). 
two regions of uniform flow, then from Eq. (3b) this 


and 


(12a) 


(13b) 


Considering, as before, only the 


second necessary relation follows: 


T, = (m/R)(C — 1m) (14) 


Eqs. (13b) and (14) are the two equations in 7, and 1, 
from which the limit values of Rankine-Hugoniot can 
be calculated for temperature-dependent gas proper- 
ties. The value of the intrinsic energy term given by 
Eq. (12) can be simply obtained by a graphical inte«ra- 
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tion of a curve of c, vs. T as given, for example, in refer. 
ence 9. Such an integration made from the data iy 
reference 9 is shown on Fig. 7 (f;! c, dT vs. tempera. 
ture). 

It is interesting to note from Eqs. (13b) and (14) that 
the variable specific heat c, is the only gas property op 
which the new limit values depend, despite the fac; 
that all the physical characteristics of the gas ar 
taken as temperature-varying. 

It is convenient to express Eqs. (13b) and (14) jn 
their nondimensional forms before carrying out the 
solution by graphical means. These nondimensional 
forms are 


Ey cape Sie . I 
- + Cel, — f= G&- (18e 
Uo- 2 2 
and 
T1/To = yoMe?Ui(Co — Ui) (14a) 
where 
U; == 0 /Uo 
Co = C/o 
C is defined by Eq. (2b, 3b). A simple form for (, 
follows from Eq. (2b, 3b)—viz., 
= C 0 0 
Go=—= f +1l= : - I (e) 


Uo MUo polo” 


Combining Eq. (e) with the well-known sonic speed 
g | I 


relation—namely, 


Qo” = Yo(Po/ po) (f) 


the value of Cj, in terms of the initial Mach Number 


Mo, becomes 
ie = l 4 (1 0M?) (g) 


After the velocity ratio U is calculated, the corre- 
sponding pressure ratio follows from the momentum 
relation (3a), which, if written in the form appropriate 
for the present example, becomes 


pi = mC — muy (15): 


If this equation is divided by po and the initial Mach 
Number J/) is introduced, then Eq. (15) can be written 
in the following nondimensional form: 


(15a) 


Pi Po = yoMo?(C = U}) 


The density ratio p: po follows from the continuity 
Eq. (la)—namely, 


pi po = 1/U (16) 


Finally, the asymptotic Mach Number 1/; beyond 
the shock wave can be calculated from the ratio 


M,/My = (141 /Q1) /Up/ Ao 
or 
AM, == MoUi[(v0 ¥1) (T; T»)] (17) 


where 71 = C,,/C,, corresponds to the temperature 7%. 
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On Fig. +, these new asymptotic limit values,* as 
calculated from the above procedure, are plotted for 
a range of initial Mach Numbers 1.0 to 7.0. These 
jimit values are computed for an initial free-stream 
temperature of O°C. and an initial pressure of one 
atmosphere. Also plotted on Fig. 4 for comparison 
are the corresponding Rankine-Hugoniot values based 
on the assumption of constant specific heat. 

An inspection of this comparison shows that the in- 
fluence on the limit values due to variable specific heat 
begins at a Mach Number of 2.0 and that the changes 
increase with the Mach Number. It is surprising to 
note that the influence of the c, variation produces an 
increase in the pressure and density but a decrease in 
the temperature. Other conclusions based on this 
comparison are that the influence of c, variable is (a) 
the greatest on the temperature and density ratios and 
(b) the least for the pressure ratio. 

THE SHOCK-WAVE STRUCTURE FOR 
CONDUCTION, 


(IV) 
TEMPERATURE-DEPENDENT VISCOSITY, 
AND SPECIFIC HEATT 


Fairly accurate numerical solutions of the one- 
dimensional flow equations can be derived by the step- 
by-step method described below. This method is 
general insofar that it can be applied for both variable 
and constant gas properties and for both uniform and 
nonuniform exit flows. An example is given here of a 
numerical solution for an initial Mach Number Mp = 
1.0, variable u, k, and c,, and uniform exit flow. 

Four equations, given previously in Sections (I) and 
(III) are required for this analysis. These are Eqs. 
(13), (3a), (4), and (1a). 

Eqs. (13) and (8a) can, by a simple transformation, 


be written as one equation—namely, 


dT - 4u Cy + mE + mCu — (mu?/2) 


= : (18) 
du 3k mu — mC + p 


After replacing p and p, respectively, in Eq. (18) by their 
values in Eqs. (4) and (la), the following equation 


* Rankine-Hugoniot values for variable specific heat, which 
are obtained with Berthelot’s equation of state, were given in 
The limit values of reference 16 agree well with 
This agreement would seem to show 


reference 16 
those given in this paper. 
that the difference of equation of state has a small influence. 

+ The procedure for this analysis follows in some points the 
work of Becker.’ In references 
10 and 15, the value of u/k is permitted to vary with the local 


See page 32 of that reference. 
flow temperature in such a way that P, remains constant at 
0.75. In Section (IV), the local flow temperature determines the 
local P, as given by the experimental data of reference 9. 
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temperature-dependent specific heats in comparison with the 
Rankine-Hugoniot values under the assumption of constant spe- 
cific heats. 


appears: 
dT 4u Coun + uE + Cu? — (u?/2) 
du 3k s—-C+h8T (aan; 
where 
C, = Ci/m = —Eo — Cuy + (tto?/2) 
(12) 


“T 
E= | C, dt 
T. (exit) 


u, k, and c, are temperature-dependent. 
The integration of Eq. (18a) leads tof 


4 f u Cou + uk + Cu? — (u*/2) “ 
3 U = Ue (exit) k  - & + RT 


du + C3 (1Sb) 


For the problem treated here, the boundary conditions 
are prescribed as uniform at both ends of the flow. 
This procedure of integrating from exit to the entrance 
then assures that the exit condition is satisfied. There- 
fore, the integrals above use as the lower limit the values 
of T = Tx, and U = U,yi, at the end of the flow. The 
reason for this will be seen shortly. As stated at the 
beginning of this section, an example will be given only 
for the case of uniform exit flow. For this case the 
boundary values of the integral (1S8b) are then taken 
as the Rankine-Hugoniot limit ratios 7; and ™ given 
on Fig. 4. If this is done and if Eq. (18b) is rearranged 
to the more convenient nondimensional form, one ob- 
tains 


t The exit values here are characterized by the velocity and 


temperature at the exit and not by the abscissa x. 


U + 'y}} 


dU+Ih (18e) 


(II/yoMo?) — U(Cy — UV) 
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The solution of Eq. (18c) is achieved by a simple 
graphical iteration procedure. A first approximatigy 
to this solution can be obtained by drawing a curye gj 
temperature ratio vs. velocity ratio through the ty, 
sets of coordinate values (U}, Il) (1, 1) corresponding 
to the flow conditions at the exit and entrance, respec. 
tively. Here it must be remembered that the aboy 
exit values must be the Rankine-Hugoniot values fo; 
c, variable. 

As an aid in drawing the first approximation cury 
the 7 vs. u relation can be used which follows from the 
exact one-dimensional solution of Becker.’ This rely 
tion is for a parabola and has the form 


(u?/2) + Cy,T = (uo?/2) + c,To = constant 
or, in the ratio form, 


ae II ‘ae . 
= 2) 14 
2 is = 1) Mo? (m= 1) MW? 


The values of u and k required* to solve Eq. (18 
must correspond to the local temperature (which 
corresponds to the chosen velocity), while a value of | 
at this same temperature can be obtained as described 
in Section (III). Fig. 7 shows a plot of EF vs. T for 
c, variable. 

‘A plot of II vs. U is shown on Fig. 5. This curve 
gives the curve calculated for Jo = 4.0 from the 
Becker relation (19) and the corresponding curve calcu 
lated by successive approximations from Eq. (1Sc). It 
is seen from this figure that Eq. (19) is a good approxi- 
mation for developing the integration of Eq. (18c). 

After the correct temperature-velocity curve is calcu: 
lated, the x-coordinates to which these values belong 
can be obtained from a direct graphical integration oi 
Eq. (3a), which first must be changed to the following 


form: 


9 
o> ey MK 


4 ™ a : 
t= / uu [RT — u(C — u)| —' du + constant 
} 


(9) 
or, in terms of nondimensional ratios, 
Am me m II : 1 
= — — et ii 4 as 
3m ar Ue Mo vol 7 l 
constant (20a 


The lower limit U, is the value of the velocity at the 
shock-wave center (where U” = 0). If it is prescribed 


that U = U, at x = O, then the integration constant 
in Eq. (20a) is zero. The integration of Eq. (20a 18 
then performed to the right and left of U = U,. Fe 


nally, a word about determining U,. 
Eq. (3a), from which Eq. (20a) followed, can also be 


written in the form 


* The values of variables yu, k, and c,; used to derive the results 
of this paper are taken from reference 9. 
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A plot of dU’ dx vs. U is made from Eq. (21) using the 
land U values calculated from Eq. (18c). On this 
plot it will be found that the velocity gradient has a 
maximum value; this maximum value corresponds to 
U" = 0, so that the velocity at this point of the dU’/dx 


vs. Ucurve is U’,. 


Following the procedure described above and using 
the II vs. U values on Fig. 5, a plot of U vs. x was calcu- 
lated and is shown on Fig. 6. Also shown on this figure 
is the corresponding shock-wave solution for constant 
gas properties which is given by the solution of Becker. 
This comparison shows that, with variable gas prop- 
erties, the influence of the shock wave extends further 
upstream and downstream; hence, the transition region 
is greater than that with constant gas properties. 
This “‘softening’’ of the shock wave is due to the large 
increase in uw and k in the region of the inflection point. 
The thickness of this shock wave is about two and one- 
half times greater than that for constant gas character- 
istics. It is also about two and one-half times greater 
than the mean free path of the molecules. 

In closing, a word about the validity of the Navier- 
Stokes equations for this problem. The usually ac- 
cepted criterion for such validation has been based on 
a comparison of the shock-wave thickness with the 
molecular mean free path. Puckett and 
Stewart !? show that the shock-wave thickness is always 
several mean free-paths thick if the mean free path is 


However, 
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dx du l ; rhey further state, based on the work of Bethe and 


Teller,'* that a better criterion of the validity of the 
Navier-Stokes equations is the number of molecular 
collisions required for statistical equilibrium to occur 
in the shock wave. It may be concluded from this that, 
as long as dissociation, condensation, etc., effects do not 
occur or the gas is not highly rarefied, the Navier- 
Stokes equations will adequately describe the shock- 
wave phenomenon. 

It would seem that the Navier-Stokes relations are 
valid for viscous shock-wave problems within certain 
obvious limits. A difficulty in this work, not previously 
mentioned, is the formidable problem of obtaining re- 
liable experimental data for comparison with the analyt- 


ical studies. 
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The Transient Temperature Distribution ina 
Wing Flying at Supersonic Speeds’ 


JOSEPH KAYE? 
Massachusetts Institute of Technology 


SUMMARY 


The motion of a vehicle at high supersonic speeds results in a 
large nonuniform temperature rise within the vehicle. Several 
new problems, such as the magnitudes of the thermal stresses 
and of the aerodynamic deflections caused by this nonuniform 
temperature distribution, must be solved by the designer. These 
problems depend, in turn, on accurate prediction of this transient 
temperature distribution 

The calculation of the transient temperature distribution is 
illustrated for the case of.a thin wedge-shaped wing or airfoil 
subjected to constant acceleration from Mach Numbers of 1.4 
The calculation depends strongly on the state of the art 
as a result, the 


to 6 
relative to heat transfer at supersonic speeds; 
calculations may be systematically in error because of the lack 
of data in this art. 

A long and laborious numerical method of solution of the dif 
ferential equation for heat flow in two dimensions with boundary 
conditions dependent on time is presented, and the resultant 
temperature distributions in the wing are given in chart form. 
In this solution the spanwise heat flow is assumed to be negli- 
gible. The results indicated that the heat-flow problem could 
be reduced to one-dimensional heat flow in the direction normal 
to the chord without sacrifice of accuracy for a wing with a maxi 
mum thickness ratio of less than 10 per cent. 

A shorter numerical method of solution of the differential equa 
tion for heat flow in one dimension—normal to the chord—is 
presented. It is shown that this method yields temperature dis 
tributions that are in excellent agreement with those calculated 
by the method for two-dimensional heat flow. This shorter method 
was used to investigate the effects of changes in the coefficients 
of heat transfer and in the values of the acceleration. It is shown 
that this shorter numerical method of solution is applicable to the 
general case of a thin airfoil, provided the variation of the local 


coefficient of heat transfer with time is available 


SYMBOLS 


A = heat transfer area 

a = velocity of sound 

Cp = specific heat at constant pressure 

G = flow per unit area or mass velocity = pV 

¢ = acceleration given to unit mass by unit force = 32.174 
ft. per sec.? 

h = coefficient of heat transfer 

k = thermal conductivity 

l = half thickness of wing 


YW = Mach Number = V/a 


Presented at the Aerodynamics Session, Eighteenth Annual 
Meeting, I.A.S., New York, January 23-26, 1950. 
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paper is an advance release of a portion of U.S.A.F. Technical 
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Jr., and Victor C. Yeh. In addition, the assistance of R. H. 
Shoulberg and Miss Barbara M. Healy is gratefully acknowledged. 
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Nu = Nusselt Number = Ax/k 


p = static pressure 

Pr = Prandtl Number = cpu/k 

Q = heat transfer 

q = rate of heat transfer 

R = gas constant 

r = recovery factor = (taw — tm)/(t, — tm) 
Re = Reynolds Number = pVx/yu 

T = temperature, °F. absolute 

t = temperature, °F. 

V = velocity 

v = specific volume 

w = mass rate of flow 

x = position coordinate normal to chord 

y = position along chord from leading edge 
a = thermal diffusivity = k/p c, 

a’ = angle between shock and incident flow 
8 = angle of attack 

€ = turning angle in expansion 

MB = viscosity 

p = density 

@ = time 

w = vector angle of characteristic curve in hodograph plane 
w’ = wedge angle for shock 


Subscripts 


aw = adiabatic wall 

m = mean stream condition 
QO = initial condition 

§ = stagnation condition 

w = wall or surface condition 


INTRODUCTION 


| yaanatse to make air-borne devices move at ever 
increasing speeds have resulted in the discovery 
of new problems that must be solved for successful 
flight. Such is the case of the aircraft or missile that is 
required to move at high supersonic speeds——that is, 
speeds of the order of several times the velocity of sound 
in air. In this case, in addition to new aerodynamic 
problems that confront the designer, he must investi- 
gate the effects of the large temperature rise of the sur- 
face of the moving object and of the resultant tempera- 
ture distribution within the object. For example, if 
the object is required to fly steadily at an altitude of 
50,000 ft. and at a speed of 5,500 m.p.h. (or a Mach 
Number of 8), then the temperature of the surface of 
the object will be of the order of 4,000°F. 

If the surface temperature of an air-borne super- 
sonic vehicle can become so large, it is evident that the 
designer will be faced with a host of new problems that 
depend on the transient temperature distribution in 


the vehicle. Investigation of this transient tempera- 
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ture distribution within the vehicle has led to results 
that are of great value to the designer. 

Consider an object moving relative to a fluid, such 
as an aircraft or a missile moving through air. The 
velocity of the fluid relative to the object is reduced 
from a finite value at a large distance to zero at the sur- 
face of the object. The temperature of the fluid rises 
from its value at a large distance to almost the stagna- 
tion temperature (the temperature corresponding to 
zero velocity of the fluid) at the surface of the object. 
If we limit our discussion to a gas that has an equation 
of state of the form 


pv = RT (1) 


then the stagnation temperature is the same whether 
the reduction to zero velocity is accomplished at con- 
stant pressure by an irreversible mechanism, such as 
shear forces in the boundary layer, or by a reversible 
process such as diffusion to higher pressure. 

If the temperature of the object is uniform and identi- 
cal with the temperature of the air at the surface of the 
object, there will be no heat transfer between air stream 
Then the temperature of the fluid at the 
adiabatic wall tem- 


and object. 
surface of the object is called the 
perature.’’ The relation between the adiabatic wall 
temperature and the temperature of the fluid at a large 


sé 


distance—the free-stream temperature—is given simply 
in terms of the recovery factor, defined as follows: 
| ae — 1 m - i — 7 m 9 
hao no ae (+) 

Sy = Tu Vu / ile 


It is seen that the recovery factor 7 is the ratio of the ex- 
cess of the adiabatic wall temperature over the free- 
stream temperature to the excess of the stagnation tem- 
perature over the free-stream temperature. Values 
of the recovery factor are based at present on relatively 
few analytical and experimental studies. 
will be discussed more fully later in this paper. 

For subsonic speeds, the adiabatic wall temperature 
of the fluid is not appreciably greater than the free- 


These values 


stream temperature. For supersonic speeds, however, 
the adiabatic wall temperature increases rapidly with 
increasing Mach Number, until at Mach Numbers 
greater than about 5 tt becomes one of the most im- 
portant variables to be considered in the design of high- 
speed vehicles. This increase in adiabatic wall tem- 
perature with increasing Mach Number is shown in 
Figs. | and 2 for flight at sea level and at the isothermal 
atmosphere at altitudes from 35,000 to about 100,000 
ft. (labeled 50,000 ft.), respectively. The two upper- 
most curves in each of these figures represent the adi- 
abatic wall temperatures for air for a recovery factor of 
unity and, hence, also the stagnation temperatures; 
the lower curves are for recovery factors of 0.9 and 0.8. 
The uppermost curve is calculated with a constant spe- 
cific heat of air, while the other three are based on a 
variable specific heat, using values of enthalpy from 
Gas Tables.' 
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The striking feature in Figs. 1 and 2 is the extremely 
large value of the adiabatic wall temperature at the 
higher Mach Numbers. Since the adiabatic wall 
temperature is usually higher than the temperature 
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TEMPERATURE 


of the object, except in steady flight, the serious nature 
of the “aerodynamic heating’ of the object, as it is 
often called, at Mach Numbers above approximately 
4 to 5 becomes immediately apparent. 

The effects of these large adiabatic wall temperatures 
are alleviated to some extent in a supersonic vehicle by 
the finite rate of heat transfer to the object and by the 
radiation from its surface. Both of these effects will 
reduce the surface temperature of the vehicle below the 
adiabatic wall temperature. It is apparent, however, 
that flight at high Mach Numbers will raise many ques- 
tions such as the following: 

(a) What is the surface temperature of the vehicle 
for flight under accelerating conditions and under steady 
conditions? 

(b) What is the temperature distribution within the 
structure of the object? 

(c) Is it possible to estimate rapidly the temperature 
distribution for flight “under conditions of constant 
acceleration ? 

(d) What are the effects of the nonuniform tempera- 
ture distribution? Or, what are the thermal stresses? 

(e) What is the effect of temperature distribution 
on the aerodynamic characteristics of the wing, such 
as angle of attack, deflections, ete. ? 

(f) What metals will be satisfactory for use in high- 
speed vehicles? 

(g) What are the cooling requirements for safe opera- 
tion of these vehicles? 

The purpose of the present paper is to present some 
answers to questions (a), (b), and (c), not for the gen- 
eral case of any object but for the special ease of a solid, 
homogeneous, thin wedge-shaped wing. It was appar- 
ent that a detailed examination of-a single shape would 
help in discovering the basic problems of aerodynamic 
heating and would help in pointing out the basic data 
needed for future investigations. The results obtained 
in this study have justified this approach to the prob- 
lem. 


REVIEW OF PREVIOUS WORK ON TEMPERATURE 
DISTRIBUTION IN AN OBJECT MOVING AT SUPERSONIC 
SPEEDS 


A careful search of the unclassified literature did not 
disclose any previous work on the temperature distribu- 
tion in a body moving at supersonic speeds. On the 
other hand, several recent papers were found on the 
skin temperature of a body moving at high speeds. 

Wood? calculated the surface temperatures of a sym- 
metrical double-wedge airfoil for steady supersonic 
flight at Mach Numbers from 2 to 10, for emissivities 
from 0 to 1, and for altitudes from 50,000 to 200,000 
ft. He attempted to include the effects of solar radi- 
ation. He did not investigate the surface tempera- 
tures for flight at variable speed—i.e., for an accelerat- 
ing vehicle. He emphasized the large number of prob- 
lems that must be solved before his assumptions or 
extrapolations could be considered satisfactory. 
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Scherrer* presented an approximate method for 
calculating the cooling requirements of a body of revolu- 
tion in steady supersonic flight based on the assumption 
of a laminar boundary layer around the entire body. 


Stalder and Jukoff! presented ‘‘a general method, us- 
ing the methods of kinetic theory, whereby the sur- 
face temperatures of bodies can be calculated for steady 
flight at any speed in a rarefied gas.’’ The study, be- 
cause of its method, was limited mainly to extremely 
high altitudes. The effects of solar radiation were also 
included, and some effects of internal cooling were 
mentioned. 


Huston, Warfield, and Stone’ compared ‘‘the time 
history of the skin temperature measured on the nose 
of a V-2 missile with the time history of the tempera- 
ture computed by use of Eber’s experimental relation 
for heat-transfer coefficients for conical bodies under 
supersonic conditions.” They developed a_ general 
numerical method for computing the time history of the 
skin temperature of a cone, taking into account the 
effects of solar heating and radiation to the atmosphere. 

Lo® presented a numerical method for the calculation 
of the transient skin temperature of conical bodies dur- 
ing high-speed flight. He included the effect of emis- 
sion from the skin but neglected the absorption of solar 
radiation by the skin. 

Johnson and coworkers’ reviewed the literature and 
summarized the data available for recovery factors and 
heat-transfer coefficients for supersonic velocities. In 
this report they use several examples to illustrate the 
calculation of skin temperatures of an object in high- 
speed flight. 

Finally, the measurements of skin temperatures on 
the nose of a V-2 rocket by Fischer and Norris* should 
be mentioned here. Measurements made during the 
transient period of variable acceleration were analyzed 
to yield some data on heat-transfer coefficients for 
laminar and turbulent boundary layers. 
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STATEMENT OF THE PROBLEM 
Consider a solid, homogeneous, thin wedge-shaped 


wing, such as that shown in Fig. 3, which is flying 
steadily at a low supersonic velocity, such as a Mach 
Number of 1.4. The wing is suddenly subjected to a 
constant acceleration until it attains a larger supersonic 
velocity, such as a Mach Number of 6, and then it is 
flown steadily at this Mach Number. What is the 
temperature distribution in the wing at any time after 
the start of the sudden acceleration? During the period 
of constant acceleration, the temperature of the wing 
is a function of both position and time. During the 
time period beginning with steady flight at a Mach 
Number of 6, the temperature of the wing is dependent 
on both position and time, but the adiabatic wall tem- 
perature and heat-transfer coefficient are now inde- 
pendent of time. When the final steady state is at- 
tained, the temperature will be, of course, independent 
of time but not of position in the wing, since the ambient 
temperatures in the four regions of the wing will have 
different values. 

We are thus concerned with the solution of the differ- 
ential equation governing the transfer of heat for two 
distinctly different types of boundary conditions. For 
the first, corresponding to the period of constant ac- 
celeration, the boundary conditions at any point on the 
wing surface are a function of time. For the second, 
starting with the transient heating process in steady 
flight at the higher velocity, the boundary conditions 
at any point on the wing surface are independent of 
time. Many solutions are available in the literature for 
the second type of constant boundary conditions, but 
only a few solutions of an extremely specialized form are 
known for the first type of variable boundary condi- 
tions. Examples of solutions of the second type of 
constant boundary conditions are available in refer- 
ence 9 in the form of the Gurney-Lurie charts and the 
Hottel charts for a large slab. Examples of solutions 
in analytical form for the first type are available in 
reference 10, but none of these corresponds to the 
boundary conditions peculiar to the wedge-shaped wing 
undergoing acceleration. 

Excellent solutions for engineering calculations are 
available in the form of the Gurney-Lurie and the 
Hottel charts for the case of boundary conditions inde- 
pendent of time. These charts may be used as fair 
approximations to obtain the transient temperature 
distribution in the wedge-shaped wing on a basis of 
one-dimensional heat flow for the period of time be- 
ginning with steady flight. The error introduced is 
caused by the fact that the initial condition used in the 
construction of the charts was one of uniform tempera- 
ture, whereas the initial temperature in the wedge- 
shaped wing is the temperature distribution corre- 
sponding to the end of the period of constant accelera- 
For large values of the acceleration, the error 
Hence, 


tion. 
made in using these charts would be negligible. 


SCIENCES—DECEMBER, 1959 
this phase of the transient heating problem was con. 
sidered solved. The major effort in this work was qj. 
rected toward obtaining solutions of the differentia| 
equation for flow of heat for the case of boundary gop, 
ditions dependent on time. 


ASSUMPTIONS 
The following assumptions, which tend to oye. 
simplify the real problem, were made: 

(1) The wing is a solid symmetrical double wedge 
without sweepback and with a constant chord, as showy 
in Fig. 3. The chord is 5 ft. 6 in., and the maximun, 
thickness is 8 per cent of the chord. 

(2) The supersonic flow of air by the wing jis gs. 
sumed to be two-dimensional, independent of the span 
For supersonic speeds above Mach Number of |4 
Fig. 3 shows the flow pattern around the wing with q 
pair of shock waves at the leading and the trailing edges 
and a pair of expansion waves at the mid-chord. 

(3) The wing loading is fixed at 100 Ibs. per sq.ft. 

(4) The wing is made of steel with constant values of 
the thermal conductivity, specific heat, and density. 

(5) The spanwise heat transfer is negligible. 

(6) Heat transfer by radiation to or from the wing 
surface is negligible. 

(7) Flight is at a fixed altitude of 50,000 ft. 

(8) The flight of the wing begins with an almost uni- 
form temperature distribution in the wing correspond 
ing to steady flight at a Mach Number of 1.4. 

These assumptions were made primarily to reduce 
the large number of laborious computations required 
for solution of the two-dimensional heat flow in the wing. 
It is believed that they do not affect seriously the final 
results on temperature distribution in the wing during 
the transient heating period. The effects of changing 
some of these assumptions were investigated later by 
the simpler solution based on one-dimensional heat 
flow. The maximum thickness of 8 per cent was se- 
lected to keep the drag low for the supersonic velocities 
studied here. The selection of steel as the material of 
the wing was recommended because of the large tem- 
peratures that would be encountered at the higher Mach 
Numbers. The effects of radiation could have been 
taken into account, but it was decided to neglect these 
in order to retain simplicity. Although the studies 
were made for flight at an altitude of 50,000 ft., they 
apply in fact, with slight modifications, to the isother- 
mal atmosphere ranging from 35,900 to about 100,000 
ft.; for simplicity, the figure legends show only 50,000 
ft. 


SOURCES OF DATA 


Properties of Air 


The properties of air were taken directly from Gas 
Tables by Keenan and Kaye.! Here, it is assumed 
that Eq. (1) is valid for air and that the specific heat of 
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TEMPERATURE 
gir is a function of temperature. The thermal conduc- 
tivity, viscosity, and Prandtl Number of air are tabu- 
lated on page 34 of Gas Tables. 


Properties of Steel 
As stated in assumption (4), the properties of the 

steel wing were assumed to be constant; this was done 
in order to simplify the laborious calculations required 
by the method of finite differences for solution of the 
two-dimensional heat flow in the wing. The following 
values were used : 

k = 22 B.t.u. per hour ft. °F. 

c, = 0.118 B.t.u. per Ib. °F. 

p = 486 lbs. per cu.ft. 


Properties of Two-Dimensional Supersonic Flow About 

Wing 

By means of the assumptions (1), (2), (3), and (7), 
the supersonic flow conditions in the four regions around 
the wing, shown in Fig. 3, were determined. The angle 
of attack was first found at each Mach Number by 
means of the slope of the curve of lift coefficient versus 
Mach Number as given by Ivey.'! With this angle of 
attack, either Gas Tables or reference 12 was used to 
obtain the flow conditions in each of the four regions 
around the wing. From the known velocities around 
the wing, the temperatures corresponding to the ambient 
conditions in the four regions, to stagnation, and to the 
adiabatic wall were easily computed for steady flight at 
a given speed. Detailed calculations are given in Ap- 
pendix (I). 

Fig. 4 presents the results of such calculations for the 
wedge-shaped wing shown in Fig. 3. The abscissa is 
the flight Mach Number corresponding to the speed 
of the wing and is based on the sound velocity in air at 
an altitude of 50,000 ft. Separate scales are also given 
showing the wing speed in units of feet per second and 
The stagnation temperature is given 
The ambient temperatures 


of miles per hour. 
as the uppermost curve. 
of the four regions corresponding to the flow pattern 
shown in Fig. 3 are shown as the four nearly flat curves 
at the bottom of Fig. 4. Since the differences in tem- 
perature are so small for these four regions, the adi- 
abatic wall temperatures calculated from these four 
ambient temperatures and from the ambient Mach 
Number were so close together that a single line has 
been employed in place of four lines on Fig. 4 to repre- 
sent the adiabatic wall temperature. A recovery 
factor of 0.9 was used in these calculations. 

The curves in Fig. 4 start at a Mach Number of 1.4 
mainly to be fairly certain that the flow pattern given 
in Fig. 3 would be valid. For Mach Numbers less 
than 1.4 it appears that the shock is detached from the 
leading edge and that this would result in somewhat 
different temperatures for the four regions of the wing. 
Since the effects of aerodynamic heating are small at 
Mach Numbers less than 1.4, the entire range from a 
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Mach Number of 0 to 1.4 was not considered in this 
report. 

Fig. 4 shows that the adiabatic wall temperature for 
a wedge-shaped wing is about 2,200°F. for steady flight 
at a Mach Number of 6. On the one hand, if one ac- 
celerates such a wing rapidly from a Mach Number 
of 1.4 to 6 and then flies steadily at 6, the temperature 
difference available for heat transfer into the wing will 
be considerable. On the other hand, if one accelerates 
the wing slowly, the temperature difference for heat 
transfer may be small. The temperature distributions 
in the wing resulting from these widely different ac- 
celerations for the same increase of Mach Number will 
be markedly dissimilar. It is evident, therefore, that 
the acceleration of the wing must be tied together with 
the data in Fig. 4 before one can establish either the rate 
of heat transfer into the wing or the temperature dis- 
tribution at any time. This is the primary reason 
why much attention has been given to the acceleration 


of the wing in this study. 


Values of the Recovery Factor 

The recovery factor has been defined as the ratio of 
the excess of the adiabatic wall temperature over the 
free-stream temperature to the excess of the stagnation 
temperature over the free-stream temperature. It is 
given simply by Eq. (2). 

The recovery factor can be measured in a stream of 
air by several methods. One method consists of plac- 
ing a plate thermometer in the stream, waiting for 
steady-state conditions, and then noting the reading 
of the plate thermometer, 7,,,,; the free-stream tempera- 
ture, 7; and the free-stream velocity, V,,. The value 
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of the free-stream temperature is usually calculated 
from a pressure measurement. The temperature of the 
plate must be uniform—.e., no temperature gradients 

either parallel or perpendicu- 


A thin plate constructed 


should exist in the plate 
lar to the direction of flow. 
from a metal of high conductivity or a well-insulated 
plate covered with a thin metallic coating are two ap- 
proximations to this plate thermometer. 


For a plate thermometer, the temperature distribu- 
tion is given schematically in Fig. 5 as curve A for an 
insulated plate. It is seen that the plate temperature 
is the adiabatic wall temperature with zero temperature 
gradient at the surface of the plate. If the temperature 
of the plate is controlled at any desired value by ex- 
traneous devices, then curves B and C indicate sche- 
matically the temperature distributions in the boundary 
layer for.plates that are hotter and colder, respectively, 
than the adiabatic wall temperature for the same fixed 
free-stream conditions. The importance of the con- 
cept of adiabatic wall temperature will be discussed in 
greater detail in the section on coefficients of heat trans- 
rer. 


Other methods of measuring the recovery factor for 
flow over objects of different shapes are also known. 
In addition, recovery factors can be measured for flow 
of air through tubes of circular cross section; it appears 
probable that the most accurate measurements of re- 
covery factor have been obtained by this method. 
Such measurements have been published. recently by 
McAdams, Nicolai, and Keenan!* for subsonic veloci- 
ties and by Kaye, Keenan, and McAdams" for super- 
sonic flow of air in a pipe. 
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Theoretical investigations have also been Carried oy 
to determine the adiabatic wall temperature, recover 
factor, and distribution of temperature within th 
boundary layer for flow of compressible and incompre 
sible fluids. Most of the results obtained analytical, 
have been for a laminar boundary layer. Several syy, 
maries of these investigations are available; one of th, 
latest can be found in reference 7. The reader js x 
ferred to this report for details of these studies; seyery 
of the theoretical temperature distributions in th, 
boundary layer are also reproduced there. 

Recovery factors have been measured for adiabati 
flow of air over flat plates, cones, wedges, and cylinders 
and for flow of air inside tubes of circular cross section 
The available data for recovery factors for flow over 
For 


subsonic flow of air over a flat plate, this report states 


various objects are summarized in reference 7, 


the situation as follows: 

(a) The values of the recovery factor are independ- 
ent of both Mach Number and Reynolds Number for 
both laminar and turbulent boundary layers. 

(b) For a laminar boundary layer, the recovery fac. 
tor is equal to the square root of the Prandtl Number— 


5.€., 
r = (Pr) 


(c) For a turbulent boundary layer, the recover 
factor is equal to the cube root of the Prandtl Number— 
1.e., 


ry = (Pr) (4 


For supersonic flow of air over objects, the amount 
of available data and theory is so small that the situa- 
tion calls for some skillful guesswork. Johnson con- 
cludes in the above report that for supersonic flow past 
bodies one should use the values of the recovery factor 
for a flat plate with the same kind of boundary layer 
This point of view is supported, in part, by the recently 
published work on recovery factors for supersonic flow 
of air through round tubes by Kaye, Keenan, and 
McAdams. '* 


ditions where the laminar boundary layer prevailed 


The recovery factors obtained for con- 


agreed within one-half of 1 per cent with a curve repre- 
senting the square root of the Prandtl Number; such 
values are in excellent agreement with those for sub- 
On the 
other hand, the recovery factors obtained from a one- 


sonic velocities with a laminar boundary layer. 


dimensional analysis for conditions where a turbulent 
boundary layer prevailed varied from about 0.9 at a 
Mach Number of 2.4 to about 0.8 at a Mach Number 
of 1.4 and neither resembled the corresponding experi 
mental values for subsonic flow nor the rough analytical 
values for a turbulent boundary layer in supersonic 
flow. 

It is evident. fram the above discussion that much 
more experimental and theoretical work on recovery 
factors is necessary before one can select accurate values 


for supersonic-flow conditions. In fact, it should be 
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ated that data for recovery factors are practically 
aonexistent above a Mach Number of 3. It was de- 
cided, therefore, to select a recovery factor of 0.9 for 
jl calculations in the present paper for the following 


reasons : 

(1) The Prandtl Number of the working fluid, air, is 
about 0.73 at temperatures corresponding to an alti- 
ude of 50,000 ft. The square root is 0.85 and the cube 
not is 0.90. For simplicity, a value of 0.9 was chosen. 

(2) The selection of 0.9 implies that the tempera- 
ture distribution curves calculated for this report tend 
to be in error on the conservative side, since the result- 
ant surface temperatures, temperature gradients, and 
thermal stresses will tend to be somewhat higher than 
the values one would have calculated with a lower value 
of the recovery factor. 

(3) The surface temperature of the wing is always 
sonuniform in flight, so that an unknown error is 
introduced by using a recovery factor corresponding to 
uniform temperature in the body. 

(4) The effects of radiation at high surface tem- 
geratures on the value of the recovery factor ere un- 


known at present. 


THE COEFFICIENTS OF HEAT TRANSFER FOR 
SUPERSONIC VELOCITIES 


Numerous investigators have studied experimentally 
and theoretically the problem of the rate of heat trans- 
fer between a surface and a fluid that moves at sub- 
sonic velocities. The results of these studies are usu- 
ally expressed in terms of an arbitrary coefficient of heat 
transfer and in terms of some simple means of evalu- 
In turn, the co- 


ating of the properties of the fluid. 
with the 


dlicients of heat transfer are correlated 
geometry, the properties of the fluid, and the type of 
fow that prevails, such as laminar or turbulent flow. 
For moderate subsonic velocities the state of the art 
isfairly well advanced. 

In contrast to this picture, little is known about heat 
transfer for flows at high subsonic and supersonic 
speeds. Indeed, the amount of theoretical and experi- 
mental work which has been done on heat transfer at 
supersonic speeds is a small fraction of that done at 
moderate subsonic speeds during the past 60 years. 

The following reasons may account for this great 
difference in work output in the two fields: 

(1) The need for solution of practical heat-transfer 
problems involving supersonic velocities has arisen only 
recently. 

(2) The need for knowledge of the behavior of super- 
sonic flows has been emphasized only recently. 

(3) The problem of predicting coefficients of heat 
transfer for supersonic flow of a compressible fluid is a 
much more complex one than that encountered at sub- 
sonic velocities. 

(4) The experimental study of heat-transfer co- 
eflicients for supersonic flows has required a simultane- 
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ous determination of recovery factors for the same 
This difficulty will be dis- 
“effective co- 


free-stream conditions. 
cussed in greater detail in the section on 


efficient”’ of heat transfer. 


Mean Coefficient of Heat Transfer 

For low to moderate subsonic velocities it has been 
customary to use the following arbitrary definition of 
the coefficient of heat transfer between a fluid and a 
surface: 


Ps (q/A)/(tw — tm) (5) 


where /1,, is the ‘‘mean coeflicient,”’ g/A is the heat flux, 
tis the temperature of the surface, and ¢,, is the average 
or bulk temperature of the fluid at a given cross section. 

For high subsonic or for supersonic velocities of a 
compressible fluid, such as air, the temperature dis- 
tribution in the boundary layer is shown in Fig. 5 for 
the cases of adiabatic flow, of heat transfer from a hot 
plate to air, and of heat transfer from air to a cold plate. 
If we apply the above arbitrary definition of the ‘“‘mean 
coefficient’’ to this set of temperature distributions, we 
see that this mean coefficient is zero for the insulated 
plate, positive for the hot plate, and negative for the 
cold plate as long as the wall temperature is greater than 
the free-stream temperature. It is evident that the 
mean coefficient is of no value for high velocities. 


Effective Coefficient of Heat Transfer 

The situation described above can be remedied by 
using another arbitrary definition of the heat-transfer 
coefficient as follows: 
(0) 


st low ) 


h. = (q A) / (te 


where hi, is the “‘effective coefficient,’ /, is the surface 
temperature during heat transfer for given free-stream 
conditions, and ¢,, is the surface temperature for adi- 
abatic flow with the same free-stream conditions—4.e., 
the adiabatic wall temperature. 

The advantages of this definition of effective coefli- 
cient are that the heat transfer must be zero for zero 
temperature difference and that the coefficient is al- 
ways positive. In addition, some experimental data 
indicate that this coefficient is approximately inde- 
pendent of the size of the temperature difference (/, 
law). A disadvantage of this effective coefficient is that 
its use requires a knowledge of the adiabatic wall tem 
perature for the same free-stream conditions. The 
adiabatic wall temperature can be obtained from Eq. 
(2) as 


tm + (17 Vm?/2gCp) 


law 


This equation requires some knowledge of the recovery 
factor r before it can be used. In the section on recov- 
ery factors, however, we found that hardly any data 
are available on recovery factors for supersonic flows. 
The absence of such data implies, in turn, that the ex- 








794 JOURNAL OF THE 


perimental data on effective coefficients of heat trans- 
fer are probably not suited for accurate calculations. 
Herein lies one of the weakest members of the struc- 
ture used in this study. An additional question that 
introduces some uncertainty in the use of the effective 
coefficient is what temperature should be used to evalu- 
ate the properties of air in any correlation of the heat- 
transfer test results. The correct answer to this ques- 
tion is not known for large temperature differences. 
For these reasons, we have treated the effective coef- 
ficient in this report more in the nature of a parameter 
whose proper value can be assessed only through future 
experimental work. 

The effective coefficient of heat transfer has been used 
throughout this report; in all that follows, the word 


“‘effective’’ will be omitted. 


Experimental Data 


The experimental data on local heat-transfer coefti- 
cients for supersonic flow have been summarized by 
Johnson, et al.’ These data were obtained from experi- 
ments with supersonic flow of air over a cone and over a 
flat plate. In addition, a recent piece of test work by 
Fischer and Norris,* by means of analysis of skin-tem- 
perature measurements on the nose of a V-2 rocket, has 
provided more information on these coefficients. 
Finally, the recent experimental work with air flowing 
at supersonic velocities in a tube, by Kaye, Keenan, 
and McAdams," has helped to clarify the situation to a 
greater degree. 

The conclusions drawn from the available test data 
are summarized as follows: 

(1) Laminar Boundary Layer.—For a laminar bound- 
ary layer in either a subsonic or a supersonic stream, the 
test data on local heat-transfer coefficients are in fair 
agreement with the theoretical work of Pohlhausen* for 
incompressible flow over a flat isothermal plate. The 
equation used for correlation of the local coefficient of 
heat transfer for subsonic and supersonic velocities 
does not differ greatly from that used by Johnson, et 


al.’—namely, 


Nu = hx/k = 0.33Re,'*Pr“ (7) 


or, in another form, 
St = h/c,G = 0.33Re,~'“*Pr-*” (8) 


where the Nusselt and Reynolds Numbers are to be 
evaluated in terms of the distance from the leading 
edge, x. 

The choice of the temperature to be used in the evalu- 
ation of the fluid properties is arbitrary for supersonic 
flows; various experimenters have used the temperature 
corresponding to the surface, to ambient conditions 
just outside the boundary layer, and to diverse com- 


binations of these two. In this paper, the properties 


* Pohlhausen’s work is summarized in English by Goldstein," 
pp. 623-631. 
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of air were evaluated at a temperature corresponding ty 
ambient conditions just outside the boundary layer 
around the wedge-shaped wing. ] 
(2) Turbulent Boundary Layer.—For a fully ¢. 
veloped turbulent boundary layer in either a subsonj 
or a supersonic stream, the test data on local coefficien 
of heat transfer are in fair agreement with the empiric 
correlation for turbulent subsonic flow over a flat plate 
namely, 


Nu = hx/k = 0.030Re,°-8Pr' q 


or, in another form, 
St = h/c,G = 0.030Re,~-°-2Pr~* (10 


The Nusselt and Reynolds Numbers are to be eyaly. 
ated in terms of the distance from the leading edge, y. 
The experimental data for supersonic flow of air in tubes 
presented by Kaye, Keenan, and McAdams" also sup. 
port the above relation, since these data for the full 
developed turbulent boundary layer are well repre. 
sented by the empirical correlation equation for pipe 
flow of McAdams;’ this correlation equation, in tur, 
does not differ greatly from Eqs. (9) and (10) for tur. 
bulent flow of air over a flat plate. 

For this paper, the temperature at which the proper- 
ties of air were evaluated for use in Eqs. (9) and (10 
was that corresponding to ambient conditions just out- 
side the boundary layer around the wedge-shaped 
wing. This choice introduces an uncertainty in the 
calculations which can be removed only when additional 
experimental data become available. 

(3) Transition from the Laminar to the Turbulent 
Boundary Layer.—The supersonic test data mentioned 
above indicate that an appreciable length of flow is re. 
quired during the transition from a laminar to a fully de. 
veloped turbulent boundary layer. Much uncertainty 
is thereby introduced in the cboice of the proper values 
of the heat-transfer coefficient to be used for a given 
object such as the wedge-shaped wing. Indeed, the 
situation is not one inviting accuracy at all, since the 
“inception of the transition’’ process is not known with 
reasonable precision. In addition, the test data for 
transition were usually obtained under steady-state 
conditions with almost constant surface temperatures. 
For a wedge-shaped wing with a variable surface tem- 
perature in accelerating flight, it is extremely doubtiul 
that one can do more than make a judicious guess as 
to the transition region, its inception, or the length of 
flow required to attain the fully developed turbulent 
boundary layer. Some theoretical work by Lees" on 
the stability of the laminar boundary layer in a com- 
pressible fluid indicates that the point where a small dis- 
turbance in the laminar boundary layer begins to grow 
to one with large amplitudes, and then causes break- 
down or instability, is dependent mainly on the Mach 
Number and the temperature of the surface relative to 
that of the air. His results indicate that this point of 
incipient instability is displaced downstream in the 
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flow if either the Mach Number is increased to that 
of a highly supersonic stream or if the temperature of 
the surface is reduced below that of the fluid; for these 
two cases he found that the Reynolds Number was not 
the important criterion for onset of the transition from 
laminar to turbulent flow. For the case of the wedge- 
shaped wing traveling at supersonic speeds, these 
considerations of stability introduce further uncer- 
tainties into the choice of the correct heat-transfer co- 
eflicient. If the laminar layer persists for a relatively 
long distance from the leading edge of this wing, then 
the coefficients of heat transfer will be reduced in mag- 
nitude and the solution of the transient heat-flow prob- 
lem will be somewhat different from that presented in 
this paper. 

In the absence of data and for the reasons mentioned 
above, it was assumed that the “point of transition” 
from a laminar to a turbulent boundary layer is the 
same for supersonic as for subsonic flow. The value of 
the critical Reynolds Number, based on distance from 
the leading edge, at transition was therefore selected 
as 500,000. This leads to a laminar boundary layer 
only for a few inches from the leading edge of the wing 
and to a turbulent boundary layer for the remainder of 
the wing. Such a selection results in using the much 
larger surface coefficients of heat transfer over most of 
the wing and in overestimating the magnitude of the 
temperature distributions and of the thermal stresses. 
If a higher value of the critical Reynolds Number had 
been selected, the results to be given later for the one- 
dimensional heat-flow problem would not have been 
affected, but those obtained for two-dimensional heat 
flow would have been changed slightly. 


Calculated Values of the Heat-Transfer Coefficient 


Fig. 6 provides a comparison of the calculated heat- 
transfer coefficients for a laminar and for a turbulent 
boundary layer on a flat plate for several Mach Num- 
bers at the same free-stream conditions of temperature 
and density. The three marked iaminar 
boundary layer were computed from Eq. (8), with the 
assumption that it is valid over the distance shown 
The remaining three curves are 


curves 


from the leading edge. 
for a turbulent boundary layer and computed from Eq. 
(10), with the assumption that it is valid from nearly 
the leading edge of the flat plate. It is evident that the 
coefficients for the turbulent layer are 
roughly eight to ten times as large as those for the 
laminar boundary layer for the same distance from 


boundary 


the leading edge and for the same free-stream condi- 


tions. 


Inability to predict the type of boundary layer on the 
wedge-shaped wing can thus lead to serious errors in 
the prediction of the heat-transfer coefficients; the 
effects of this uncertainty on the final results will be 
discussed later in greater detail. The properties of air 
were evaluated at an altitude of 50,000 ft. for a tem- 


DISTRIBUTION 


(Btu/hr The °F’ 


COEFFICIENT OF HEAT TRANSFER 











































































IN A SUPERSONIC WING 795 
160 
140 | 
Free Stream Condition at 
50,000 Feet 
120 1 + = 
Turbulent: “ =0.030 Re?2pr 2/3 
CpG x 
Yi 
‘oolt / \ 
os y Mach Number 
ion 
~~ | —— 6.0 
60 ft — 
Se 
Ts 4.0 
40 —— 
\ -— 
20 1.4 
: h “8 ais 
&6 = 0.33 Re," Pr 6.0 
4.0 
ie) 1.4 
°O I 2 3 “ 5 


DISTANCE FROM LEADING EDGE OF PLATE (Ft) 


Comparison of Heat Transfer Coefficients for 








































Steady Flow with Laminar and Turbulent Boundary 
Layers on a Flat Plate. 
Fic. 6, 
280it-— T if T | l 
a: 
l 
h 
' | 
— 246 + + { } 
= rT 
° i] j 
x 
saad | |Surface 
£ 200 : 1 T 1 
3 Bottom 
re) 
r 160 
WwW 
Ww 
77) \ 
2 \ 
Z =a 
120 a 5 eecmmiaees' aa 
z \/ | 
= \ Mach Number 
WwW \ 6 | 
~ N 
| - 
ue 80 ; ) ae 
= | | 4 
= aol t— iBottom and top | ed 
<, : 
i on oe ed 1.4) 
Ww 1 T Fs 
S | ‘Turbulent 
1@) i las 
1 2 3 4 S$ 3&5 


1?) 
DISTANCE FROM LEADING EDGE OF WING (ft) 


Leading Edge 


Calculated Values of Coefficients of Heat 
Transfer for Flight of Wedge Shaped Wing at 


Constant Mach Number at 50 
Fic. 7. 


,OOO Feet 

















796 JOURNAL OF THE AERONAUTICAL 
200 | 
cs le Samael : a 
} Station 
| | 
c 
120 ; | me 








@ 
°o 


oN, 
aA 





> 
° 




















COEFFICIENT OF HEAT TRANSFER(Btu/hr ft? °F) 














° MACH 
a 5 6 NUMBER 
| : - — Ft/sec 
1000 2000 3000 4000 5000 6000 
} . T 7— MPH 
1000 2000 3000 4000 
SPEED 
Station a b 
WING 











lL. c d ae 
— 2.30" - 730 


Calculated Values of Coefficients of Heat 
Transfer for Four Selected Stations on Wedge 
Shaped Wing for Flight at 50,000 Feet. 
Fic. 8. 















Leading Edge 


Trailing Edge 
. y 
re —————— ——os 


Isometric Sketch of Wedge Shaped 
Wing Showing Selection of Coordinate System. 
Fic. 9. 





perature of —67°F. and a density of 0.01161 Ib. per 
cu.ft. 

Fig. 7-presents the variations of coefficient of heat 
transfer actually used in the calculation of the tempera- 
ture distribution of the wing. The lowest curve, repre- 
senting the variation along both the top and bottom 
surfaces for flight at a constant Mach Number of 1.4 
at an altitude of 50,000 ft., is actually the average of 
two curves lying close together. The other curves for 
Mach Numbers of 4 to 6 show the difference between 
the coefficients for the top and bottom surfaces of the 
wing. Each curve starts with a branch corresponding 
to the laminar boundary layer; this branch is followed 
by a dashed line that represents the transition region for 
a critical Reynolds Number of 500,000; finally, the 
almost discontinuous break at the mid-chord repre- 
The variation 
of the coefficient at the mid-chord is not known, but it 


sents the effect of the expansion wave. 
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is taken to be continuous since, in the absence of exper 
mental data, it is assumed that the boundary layer . 
continuous in the vicinity of the expansion waves, 


Fig. 8 shows the variations of the heat-transfer © 
efficient for four selected points on the wing for flight 
at several speeds at 50,000 ft. The relation betwee, 
heat-transfer coefficient and Mach Number is nearly 
linear for these four selected points and for most point 
along the surface of the wing not too close to either th, 
leading or trailing edges. In some of the calculatioy 
of this paper it has been assumed that this linear yay. 
ation is valid; it has helped in reducing both the labo, 
and time required for the solution of the problem. |; 
is evident from the discussion of the experimental dat, 
on heat-transfer coefficients for supersonic velocities 
that this assumption of linear variation of coefficient 
with Mach Number will not seriously affect the a¢. 
curacy of the calculated temperature distributions jy 
the wedge-shaped wing. 


In the computations of the coefficients of heat trans. 
fer in Figs. 7 and 8 for the wedge-shaped wing, Eqs. 
(S) and (10) were used for the laminar and turbulent 
boundary layers, respectively. For each of the four 
regions around the wedge, shown in Fig. 3, the ambient 
temperature, pressure density, and Mach Number were 
calculated from the wing loading, angle of attack, lit 
coefficient, etc. The details are given in Appendix (I), 
At any point along the surface of wing, the coefficient 
of heat transfer was then computed from the appro- 
priate equation using the ambient properties of air 
and the ambient Mach Number. Since the ambient 
conditions change downstream of the shock waves at 
the leading edge and downstream of the expansion 
waves at the mid-chord, the values of the heat-transfer 
coefficients shown in Figs. 7 and 8 do not correspond 
exactly to those shown in Fig. 6 for flow over a flat 
plate with the same free-stream conditions at the 
several Mach Numbers. 
TEMPERATURE DISTRIBUTION IN 
THE WING 


Two-DIMENSIONAL 


Differential Equation 


Employing the rectangular coordinate system shown 
in the isometric sketch of the wing in Fig. 9, the differ- 
ential equation for heat flow in the solid wing is 


eT eT oT 107 


(11 
Ox? Oy’ r 3" a 06 


Since we have assumed that the spanwise heat flow is 
zero, Eq. (11) reduces to the following one for the two- 
dimensional heat flow: 


(0?7/Ox?) + (0°T/dy?) = (1/a)(OT/20) (12) 


In addition, the adiabatic wall temperature and the co- 
efficient of heat transfer are functions of position along 
the surface of the wing and of time, so we have 
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TEMPERATURE DISTRIBUTION 
(13) 


fo(x, y, 8) (14) 


The boundary conditions required for the solution of 
Eq. (12) are so complex that an analytical solution is 
impractical. Hence, a long and laborious stepwise 
method was employed to obtain the transient tempera- 


ture distribution in the wing. 


Calculations 

The stepwise numerical solution of Eq. (12) was pat- 
terned after the method described by Emmons.'” '* 
method of 


’ 


This method is the so-called ‘relaxation’ 
Southwell’? when applied to the case of steady-state 
problems. In the case of the wedge-shaped wing, it was 
divided into a set of 113 stations, connected by rods, 
along which the heat flow was concentrated. For this 
geometry, illustrated in Fig. 10, the adiabatic wall 
temperature and heatstransfer coefficients were com- 
puted at each of the surface stations for a series of time 
intervals for a constant Mach 
Number every 30 sec. (32.4 ft. per sec.”), very nearly 


acceleration of one 
The details of these stepwise calculations are 
It required a total of 52 time 


lg. 
given in Appendix (I). 
snapshots to cover the range from an initial Mach 
Number of 1.4 to one of 6. 


Results of Calculations of Two-Dimensional Temperature 

Distribution in Wing 

The results based on the numerical method of finite 
differences for solution of the heat-flow equation in two 
dimensions are shown quantitatively in Figs. 11 to 19. 
In order to keep the number of charts within reason, 
the transient temperature distribution is shown at 
selected Mach Numbers for several planes and _ sur- 
faces and at selected planes or stations for all Mach 
Numbers. 

One important feature is immediately visible in these 
charts: The maximum temperature occurs at the 
leading and trailing edges of the wing and is equal to 
the adiabatic wall temperature at a Mach Number of 
6, or to about 2,220°F. Because of this high value of 
the temperature, the leading and trailing edges must be 
given special attention in the design of any thin sharp 
wings flying at high supersonic speeds. The transient 
temperature distribution within 1 or 2 in. from either 
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the leading edge or the trailing edge constitutes a 
special problem not treated here. If it is desired to 
yse, for the leading and trailing edges, a metal other 
than steel (such as titanium) or a good thermal insulator 
(such as a ceramic), then the transient temperature 
distribution could be obtained in terms of the proper- 
ties of these materials by numerical methods similar 
to those used here. 

These charts also show that the temperature of the 
bottom surface of the wing is always larger than the 
temperature of the top surface. This difference in 
temperature is a result of the larger coefficients of heat 
transfer on the bottom surface of the wing as compared 
to those on the top surface. In addition, both surface 
temperatures are considerably larger than the mid- 
plane temperature ; the maximum difference in these 
temperatures occurs near the mid-chord, where one 
would naturally expect to find the largest thermal lag 
of the wing. : 

These charts also illustrate the significant difference 
in chordwise and vertical temperature gradients in the 
The calculations indicate that the tempera- 


wing. 
the vertical 


ture gradient perpendicular to the chord 
direction—is several times as large as the chordwise 
temperature gradient, except perhaps in the immediate 
vicinity of the leading and trailing edges. This large 
difference in temperature gradient is important in this 
study, for it indicates that we would be amply justi- 
fied in treating the heat flow as a one-dimensional rather 
than as a two-dimensional problem. 

Figs. 11, 12, and 13 show the temperature distribution 
in the chordwise direction of the wing along the top 
and bottom surfaces and along the mid-plane or chord 
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Temperature Distribution at Midchord from 
Calculations of Two Dimensionol Heat Flow for Constant 
Acceleration of Ig from a Mach Number of 1.4 at 
50,000 Feet 
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Temperature Gradients Along Vertical Planes 
in Wing from Calculations of Two Dimensional Heat Flow 
for Constant Acceleration of |g from Mach Number 1.4 
to 6 at 50,000 Feet. 

Fic. 18. 
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Temperature Gradients Along Horizontal Planes 
in Wing from Calculations of Two Dimensional Heat Flow 
for Constant Acceleration of Ig from a Mach Number of 
L4 to 6 at 50,000 Feet 

Fic. 19. 


for selected Mach Numbers. In each chart, the dashed 
line at the bottom represents the initial uniform tem- 
perature of the wing for steady flight at a Mach Number 
of 1.4, and the dashed line at the top represents the 
maximum temperature available for heat flow at the 
given time, or the adiabatic wall temperature. Thus, 
a comparison can be quickly made of the rise in tem- 
perature at a given point of the wing at various times of 
flight and of the temperature difference that is causing 
heat flow at these Mach Numbers. 

In Fig. 11, the temperatures are shown for Mach 
Numbers of 3 and 4, after constant acceleration of lg 
from the initial Mach Number of 1.4. It is apparent 
that the surface temperatures are rising slowly, except 
in the immediate vicinity of the leading and trailing 
edges. In addition, there is a considerable lag in the 
temperature rise along the mid-plane except near the 








800 JOURNAL OF THE AERONAUTICAL SCIENCES—DECEMBER, 1950 


leading and trailing edges. Figs. 12 and 13, for Mach 
Numbers of 5 and 6, respectively, iliustrate this be- 
havior even more clearly. The temperature distribu- 
tions for the upper and lower surfaces show distinct 
dips in the vicinity of the mid-chord of the wing. 
These dips result from the sudden dips in the surface 
coefficients of heat transfer, shown in Fig. 7, which are 
caused by the change in ambient properties of the air 
in the rear half of the wing after passing through the 
expansion waves at the mid-chord. It will be noticed 
that there is no dip in the same region for the various 
curves of mid-plane temperature. The absence of 
the dips here is explained, in part, by the thermal 
lag of the mid-plane which would tend to smooth out 
such a dip and, in part, by the use of a “‘coarse’’ net- 
work of 113 points in the two-dimensional calculations 
of heat flow. 


The two-dimensional temperature distributions could 
also be drawn at selected Mach Numbers along planes 
parallel to the four surfaces of the wing. This, how- 
ever, was not attempted here because the number of 
charts would soon become unreasonably large. Hence, 
the temperature distributions were plotted for selected 
planes parallel and perpendicular to the chord. These 
planes correspond to the vertical planes labeled by 
letters A to O and to the horizontal planes labeled by 
numbers | to 15 in Fig. 10. It was again necessary to 
select only a few interesting charts from the total 
number available. 

Figs. 14, 15, and 16 present the vertical temperature 
distributions for several Mach Numbers at planes 
H, F, and D, shown in Fig. 10, and corresponding 
to the mid-chord, to about 2 ft. and 1.2 ft. from the 
leading edge, respectively. Similar curves would be 
obtained if we plotted the temperature distributions for 
vertical planes corresponding to the same distances from 
the trailing edge. The scales used in Figs. 14, 15, and 
16 are the same, but the decreasing width of the abscissa 
corresponds to the decreasing thickness of the wing as 
we move from the mid-chord to the leading edge. 
These curves illustrate again the thermal lag of the in- 
terior of the wing at all planes perpendicular to the 
chord except near the edges. Curves of temperature 
distribution along planes parallel to the chord are not 
included here, since they would be similar to those 
shown in Figs. 14, 15, and 16. 


Fig. 17 shows the temperature distribution at the 
mid-chord as a function of Mach Number and of time. 
The temperature at the mid-plane rises slightly above 
the initial value. The difference between the adiabatic 
wall temperature and either top or bottom surface tem- 
perature is much larger than the temperature difference 
between either surface and the mid-plane. This effect 
is due to the much greater resistance to heat transfer 
along the air-to-metal path, as compared to the thermal 
resistance of the metal itself, and to the fairly large 
acceleration of lg used in the calculations. 


Effect of Variable Acceleration 


Although the effect of variable acceleration was not 
studied by the numerical method used for two-dimey, 
sional heat flow, the limiting magnitudes of the acceler, 
tion are evident in Fig. 17 for finite values of the q, 
efficient of heat transfer. If we consider the limit of ar 
infinitesimal acceleration, then the temperature (jc 
tribution in the wing is that for steady flight at th 
given Mach Number; this is nearly a uniform tey, 
perature that is equal to the adiabatic wall temperatyy 
Hence, the curve so labeled in Fig. 17 represents yer, 
nearly the temperature distribution for zero acceler, 
tion with a finite heat-transfer coefficient. If we cop. 
siler the other limit of an infinite acceleration, the, 
the temperature distribution throughout the wing wij 
be a straight line representing the initial temperature 
at a Mach Number of 1.4 as long as the heat-transfe 
coefficient is less than infinite. Therefore, it is ey. 
dent that, as the acceleration decreases from infinity 
to zero, all the resultant temperature distributions must 
lie between the curve for the adiabatic wall temperature 
and the straight line representing the initial tempera. 
ture. These curves are so labeled in Fig. 17. 


Temperature Gradients in Wing 


Figs. 1S and 19 show approximate values of the tem- 
perature gradients along vertical and horizontal planes 
in the wing, respectively, at a Mach Number of 6 after 
constant acceleration of lg. The scales in Figs. 18 and 
19 are different. Inspection shows that the tempera- 
ture gradients in a vertical plane are larger by a con- 
siderable factor than those in the horizontal plane, with 
the exception of those near the leading and trailing 
edges. Since these calculations of two-dimensional 
heat flow are not valid at the edges because of the 
“coarse’’ network of points used, we feel that the edge 
problem should be treated separately. It should be 
noted that, in Fig. 18, the temperature gradients for 
the rear half of the wing are practically the same as 
those for the front half for the corresponding vertical 
planes and therefore are omitted from this figure. For 
the same reason, only the temperature gradieats for 
horizontal planes in the top half of the wing are shown in 
Fig. 19. We note from Figs. 1S and 19 that the two 
families of curves representing the temperature gradi- 
ents along vertical and horizontal planes are slightly 
asymmetrical about the mid-chord and the mid-plane, 
respectively. In addition, the positions of zero tem 
perature gradients are displaced slightly from _ the 
centers of the mid-chord and of the mid-plane; the 
small displacements are caused by the slight degree of 
asymmetry of the two-dimensional heat flow in the 
wedge-shaped wing. In conclusion, it is evident from 
Figs. 18 and 19 that the two-dimensional heat flow can 
be successfully reduced to a one-dimensional problem 
without introducing serious errors in the resultant tem- 
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ture distributions; both the large ratios of tem 


pera 
perature gradients and of areas normal to the direction 


of heat flow support this result. 


ONE-DIMENSIONAL HEAT FLOW IN WING 


The stepwise numerical method used in the calcu 
ation of the transient two-dimensional temperature 
distribution in the wing was a long and tedious process. 
it was natural, therefore, to seek shorter methods of; 
calculation for the solution of the differential equation of 
heat flow with boundary conditions dependent on time.} 
Several shorter methods were found and used; all were 
based on the assumption that the temperature distribu- 
tion along a line perpendicular to the chord of the wing 
can be obtained from the one-dimensional heat flow 
In other words, the chordwise heat 


along this line. 
The results from the previous 


fow was neglected. 
calculations of the twosdimensional heat flow showed 
that this assumption would probably introduce neg 
ligible errors at all points in the wing except near the 
leading and trailing edges. This prediction has been 
verified from the calculations based on heat flow in one 
dimension for the thin wedge-shaped wing. 

Many important advantages were obtained upon sim- 
plifying the problem to heat flow in one dimension. 
Both numerical and analytical solutions were found. 
The numerical solution required only a small fraction 
case. In addition, the number of variables which 
could be studied by this shorter numerical solution was 
practically unlimited; these variables included the 
thickness, surface coefficient, acceleration, type of ma- 
terial, etc. 


Analytical solutions of heat 


the one-dimensional 


flow were found based on special assumptions of 
the variation of coefficient of heat transfer with time. 
In order to keep the size of this paper within reason- 
able bounds, the presentation of the analytical solu- 


tions will be given at a later date. 


Numerical Method of Solution of One-Dimensional Heat 

Flow 

Consider first the special case of the wedge-shaped 
wing shown in the sketch in Fig. 9. This wing is flying 
steadily at a Mach Number of 1.4 at an altitude of 
0,000 ft.; it is suddenly subjected to a constant ac- 
celeration of 1g until it reaches a Mach Number of 6; 
then it is flown steadily at a Mach Number of 6. For 
the transient heat flow during the period of constant 
acceleration, we will show how the temperature dis- 
tribution can be estimated rapidly and quantitatively 
by the assumption of one-dimensional heat flow per- 
pendicular to the chord. In addition, the resulting 
method can be easily generalized to apply to any thin 
wing, provided the variation of the local surface 
coefficient of heat transfer with the time of flight is 


known. 
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Fic. 20. Network for one-dimensional heat flow 


The mathematical problem is greatly simplified by 
the assumption of one-dimensional heat flow. Using 
the coordinate system of Fig. 9, the differential equation 
is 

O*t/Ox? = (1/a) (Ot/08) (15) 
The heat flow is thus along a thin section normal to the 
chord, as shown in Fig. 20 where | is one-half the wing 


thickness at the given point along the chord. The 
three boundary conditions are then as follows: 
(1) At@é = 0,t = tay (16) 
(2) Atx = 0, Ot/ox = 0 (17) 
(3) Atx = 1, k (Ot/Ox) = h(taw — ts) (18) 


We assume in condition (1) that the initial temperature 
throughout the thin section is uniform and equal to the 
adiabatic wall temperature for steady flight at a Mach 
Number of 1.4. In addition, since the local coefficient 
of heat transfer and the adiabatic wall temperature are 
functions of time, we can also write 

(19) 
(20) 


h = f3(0) 
lou - fs(8) 


Since the position of this thin section from the leading 
edge is considered fixed, # and ¢,,. depend only on the 
time of flight. It is primarily for this reason that we 
can generalize this treatment to the case of any thin 
wing and solve for the temperature distribution, pro 
vided the functions in Eqs. (19) and (20) are available 
for a given point on the wing. 

The numerical method of solution of Eq. (15) with 
the boundary conditions given in Eqs. (16), (17); and 
(18S) was similar to that described by Emmons'™ ' and 
McAdams.* The thin normal to the chord 
was subdivided into a finite number of parts, and the 
temperature rise of each part was computed by a repe- 


section 


titious stepwise procedure covering the entire time 
period of constant acceleration. For large values of the 
acceleration, a small number of time snapshots was 
the entire calculation then required a few 


sufficient; 
For small values of the acceleration, many 


man-hours. 
nore time snapshots were needed, and the time of cal 
culation ran to as much as 8 man-days per tempera- 


ture distribution. 
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Comparison of Temperature Distributions 
Calculated from Heat Flow in One and Two 
Dimensions for Vertical Planes G,F,and E for 
Constant Acceleration of |g from a Mach 
Number of |.4 to 6 at 50,000 Feet 


Fic. 21. 


Results of Numerical Method of Solution of Heat Flow in 
One Dimension 


The first objective of the numericai method of solu- 
tion based on heat flow in one dimension was to discover 
the magnitude of the error that results from this sim- 
plification. This objective was accomplished by mak- 
ing a sufficiently large number of calculations with this 
shorter method so that comparisons could be made with 
the temperature distributions obtained from the time- 
consuming calculations based on two-dimensional heat 
flow for the same set of flight conditions and of vari- 
ation of local coefficient of heat transfer with time. 
The second objective was to investigate the variation 
of a number of important variables, such as local co- 
efficient of heat transfer, acceleration, thickness of the 
wing, size of the network in calculations, etc. 

Fig. 21 presents a comparison of the temperature dis- 
tributions calculated from the numerical methods of 
solution for heat flow in one and in two dimensions. 
The temperature distributions shown here are only for 
the top half of the wing for vertical planes G, F, and 
E shown in Fig. 10, so that the station numbers refer 
only to those in the top half of the wing. The solid 
curves in Fig. 21 represent the temperature based on 
two-dimensional heat flow and are identical with similar 
The 
triangles represent the temperatures based on the nu- 


lines shown previously. circles, squares, and 


merical method of solution for one-dimensional heat 


flow. The comparison of the temperatures is also made 
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in the form of the deviation plot under the main chart 
in Fig. 21. Here the deviations of the Lem peratures 
calculated from one-dimensional heat flow from thog 
based on two-dimensional heat flow are plotted againg 
station number. This deviation plot was required j; 
order to show the small differences in temperatures cal 
culated from the two different methods. 

The comparison, shown in Fig. 21, of the temperatyy, 
distribution based on the two methods of calculation 
was made for identical conditions of flow, acceleratioy 
wing geometry, coarseness of network, etc., and for th, 
same variation of coefficient of heat transfer with 
time. The largest deviation is less than 10°F. — Simila 
comparisons were made of the temperature distriby 
tions for Mach Numbers of 3, 4, and 5, but the mayj. 
mum deviation occurred in the comparison shown jj 
Fig. 21 for a Mach Number of 6. The deviations 
plotted in this figure show a small systematic trend 
that is most evident at the centerline station 8: the 
points deviate increasingly from the solid lines as we 
1.€., from 
since heat flow from the hotter lower half 


move toward the leading edge of the wing 
GtoFtoE 
of the wing was present in the two-dimensional calcu. 
lations but not in the one-dimensional points. The 
erratic scattering in Fig. 21 is believed to be due to com 
putational errors. The conclusion is that the sim- 
plification to heat flow in one dimension introduces a 
negligible error compared to the temperatures calcu- 
lated from two-dimensional heat flow with the same 
coarseness of network of stations. 
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Comparison of Temperature Distributions 
Calculated from Heat Flow in One and Two Dimensions 
for Vertical Plane H at Midchord for Constant 
Acceleration of |g from a Mach Number of |.4 to 6 
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TEMPERATURE 


Fig. 22 presents a comparison of the temperature dis 
tributions, along vertical plane H at the mid-chord, 
calculated from the numerical methods of solution 
for heat flow in one and in two dimensions. The tem 
perature deviations based on the two-dimensional 
method are also given. The comparisons were made 
jor identical conditions of flow, acceleration, wing 
geometry, coarseness of network, ete., and for the same 
variation of coefficient of heat transfer with time. 
These comparisons are indicated by the solid curve for 
the two-dimensional case and by the circles for the one 
dimensional case. In addition, another comparison 
has been added to this chart to illustrate the effect of 
replacing the triangular top of the section centered at 
plane H by a flat surface coincident with horizontal 
plane | for the method of one-dimensional heat trans 
fer: this change in geometry results in a change of 
volume of the material assigned to the surface station; 
the temperatures calculated with this geometry are 
shown by the squares. Comparisons similar to those 
shown in Fig. 22 were made for Mach Numbers of 3, 
i, and 5, but the extreme deviations were found to 
occur for a Mach Number of 6. 

For the same conditions of geometry the deviations 
shown in Fig. 22 are less than 10°F. for a Mach Number 
of 6. If the surface geometry is changed by replacing 
the triangular top of the center of the wedge-shaped 
wing with a flat surface passing through the apex, then 
the surface temperature based on one-dimensional heat 
flow will be about 40°F. lower than that obtained from 
the two-dimensional calculations; for the stations in- 
side the wing the deviation will be considerably less. 
The deviations found for this change in surface geom- 
etry are caused by the comparatively large change in 
volume of the material assigned to the surface station 
when the same coarseness of network is used. 

The comparisons in Figs. 21 and 22 establish quan- 
titatively that the complicated calculation of two- 
dimensional heat flow can be replaced by a number 
of simpler calculations based on heat flow in one di- 
mension for the same local coefficients of heat transfer. 
This simpler calculation could be applied with success 
to the temperature distribution near the leading and 


trailing edges if the local heat-transfer coefticients were 


available. 


Effect of Variable Coefficient of Heat Transfer 


The simple numerical method of solution of heat flow 
in one dimension was used to investigate the effect of 


changes in several variables. Similar studies could 


not have been made in the time available if only the 


method based on two-dimensional heat flow had been 
used. In addition to providing results quickly, the 
simpler method proved to be a valuable aid in guiding 
the analytical studies that were carried on simultane- 
ously. The first variable studied was the local coeffi 
cient of heat transfer: the comparison of surface tem 
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Comparison of Surface Temperatures 
for Several Variations of Heat Transfer 
Coefficient with Time Calculated from Heat 
Flow in One Dimension for Constant Acceleration 
of 2g at 50,000 Feet 


Fic. 23. 


peratures for several variations of coefficient with time 
or Mach Number is given in Fig. 28 for an acceleration 
of 2g. Similar charts were drawn for stations inside 
the wing, but the differences in temperature between 
the curves for the three variations of coefficient were 
less than those shown for the curves for the surface 
station. 

All temperature curves for various values of the 
curves for zero 


coefficient le between the limits of 


and for infinite coefficient; for the case of zero co- 
efficient, the surface temperature is the initial tem 
perature and, for infinite coefficient, is the same as the 
adiabatic wall temperature. The three linear rela 
tions between coefficient and Mach Number used in the 
calculations are good approximations for three different 
positions along the wedge-shaped wing to the actual 
variation of coefficient with time. It is evident that 
the temperature distributions cannot be calculated ac 
curately unless the relation of coefficient of heat transfer 


with time is also given accurately. 


Effect of Variable Acceleration 

Fig. 24 presents a comparison of surface tempera 
tures for the entire range of values of the acceleration 
calculated from the numerical method of solution of 
one-dimensional heat flow. Similar charts were con 
structed for the temperature at internal stations; on 
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these charts the curves for the various values of the 
acceleration tended to crowd together more than those 


on Fig. 24. For infinite acceleration the surface tem- 
perature remains equal to the initial temperature. 
For zero acceleration, the surface temperature equals 
the adiabatic wall temperature; hence, the tempera- 
ture curves for finite acceleration lie between these two 
limits. It is seen from Fig. 24 that the selection of an 
acceleration of lg for the calculation of the 


two-dimensional heat flow resulted in temperature 
values about halfway between the extreme values pos- 


sible. 


Effect of Size of Network 


The effect of change of size of network was studied 
by means of the numerical method for heat flow in one 
dimension. It was found that halving the length of 
the rods connecting the stations shown in Fig. 10-—or 
in other words, increasing the number of stations in 
the one-dimensional case from 6 to 12—resulted in 
a maximum difference of 17°F. in the temperaiure of 
the surface station for the extreme case of a Mach 
Number of 6. 
in the value of the coefficient of heat transfer alone 


The error introduced by the uncertainty 
would probably be larger than this figure. Hence, 
it is concluded that the numerical method of solution 
with the coarse‘ network used in the one-dimensional 
heat flow is a satisfactory solution to the problem, 
with an error probably much less than that introduced 
by the basic data. 
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CONCLUSIONS 


Calculations of heat flow in two dimensions for a Wing 
flying at supersonic speeds, with boundary condition, 
dependent on time from a Mach Number of 1.4 to one 
of 6, have shown the following: 

(1) The leading and trailing edges of the wing ay 
the most severely affected by aerodynamic heating, 

(2) The temperature distribution in the wing as ’ 
function of time of flight. 

(3) The temperature gradients at a Mach Number 
of 6 are considerably larger in the vertical direction 
than in the chordwise direction. 

(4) The flow of heat could be approximated with 
small error by one-dimensional flow. 


Calculations of heat flow in one dimension with 


boundary conditions dependent on time—namely. 


perpendicular to the chord—have shown the following 


(1) The error made in assuming heat flow in one 
rather than in two dimensions is negligible for thin wings 
for practical purposes. 

(2) The effects of variable coefficients of heat trans. 
fer, of variable acceleration, and of size of network have 
been examined. 

These calculations are strongly dependent on the state 
of the art of heat transfer at supersonic speeds. Since 
the state of the art is poor, much experimental work js 
needed to verify the basis of these calculations. 

Finally, these calculations of transient temperature 
distribution in the wing form the basis for the calcula- 
tion of thermal stresses and aerodynamic deflections of 
the wing. 
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APPENDIX (1) 
DIMENSIONS 


Supersonic Flow About Wedge-Shaped Wing 


Consider the supersonic flow pattern shown in Fig. 
25. There are four regions of the flow path wherein 
the properties of air are needed for the calculations 
of two-dimensional heat flow; these regions are indi- 
cated on this figure by subscripts 1, 2, 3, and 4. The 
symbols used, given in the beginning of this paper, are 
the same, in general, as those given on page 129 of refer- 
ence 1. In the following sample calculations, the data 
in references | and 12 were used for the properties of 


supersonic flow. 


Sample Calculations for Wing Flying at a Mach Number 
of 2 


Assumptions: (1) Mach Number is 2; (2) steady 
flight at 50,000 ft.; (3) wing loading is 100 Ibs. per 
sq.ft. 


Calculations: From Table 60 of Gas Tables, for 50,000 
ft. we get 


Ty) = —67°F. = 393°F. absolute 
po = 0.0116 Ib. per cu.ft. 
po = 1.69 Ibs. per sq.in. absolute 


From Table 2 of Gas Tables, the velocity of sound at 
393°F. absolute is 
do = 971 ft. per sec. 
The velocity is then 
Vo = Meo = 


1,942 ft. per sec. 


From the definition of lift coefficient we get 
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Sketch Showing Details of Supersonic Flow Pottern in Four 
Regions About Wedge Shaped Wing 


Fic. 25. 
: W/S 100 2 
c, = —, = = = 0.147 
pV?/2¢ (0.0116) (1,942?) /(2) (32.2) 


From Fig. 6 of reference 11, we find, fora Mach Num 
ber of 2, that the slope of the curve of lift coefficient 
versus angle of attack is 
dC,/dB = 0.041 
Hence, the angle of attack is 
B= 0.147/0.041 5.08 


The symmetrical double-wedge wing has a maximum 
thickness of 8 per cent. This results in a half angle of 
59°. The wedge angle for the shock at the upper 
leading edge of region | is then 


w,’ = 4.59 


With this value of the wedge angle and the value of 
Mach Number upstream of the shock (/p here and J/, 
in Gas Tables), we enter Table 56 of Gas Tables, and, 
for the angle between the shock and the incident flow, 
we find 

a,’ = 30.8 


From Table 57 of Gas Tables, with w’ and a’, we find 
the Mach Number downstream of the shock in region | 


From Table 58 of Gas Tables, we get the ratio of down- 
stream to upstream pressure for the shock and there 


fore the pressure in region | 


pi/po = 1.06 


pf: = (1.06) (1.69) = 1.79 Ibs: per sq.in. absolute 


From Table 59 of Gas Tables, we get the ratio of down 
stream to upstream density for the shock and therefore 
the density in region 1 


Pi, po = l O04 


pr = (1.04)(0.0116) = 0.0121 Ib. per cu.ft. 


The temperature in region | is obtained by applying 
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Fics. 26a, 26b, and 26c. Details of network for two 
dimensional heat flow. 


the perfect gas relationship for air as follows: 


T; To — (pi Po) ( px pi) 
pi = (1.06)/(1.04) = 1.02 
and so 
T, = (1.02)(393) = 399°F. absolute 


The velocity is given by 
Vi/ Vo = (My Mo) V 11/T, 
and so 


1.97 \ 299 
Vv, = 1,942 ( "= 1,920 ft. per sec. 
2.00) V303 


All ambient properties of air in region | are now known. 


The air passes through an expansion wave in flow 


ing from region 1 to 2. From Table 54 of Gas Tables 


and the Mach Number before the expansion, .1/;, we 
find the vector angle of the characteristic curve in the 


hodograph plane 


wo = 25.50 
Since the airflow changes its direction through the ex 
pansion, we get, for the vector angle of the character 


istic curve in the hodograph plane for region 2, 
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wo = wi + & 25.5" +- 9.18 34.7 


From Table 54,' the Mach Number downstream of the 
expansion in region 2, J/., is found from this valye oj 


M, = 2.32 


From Table 54,' we get 


Ave 1.97, Pr Po 0.154, Pi’ Po 0.238 
Mz = 2.32, po/po = 0.0795,  p2/po = 0.161 
Hence, we find 
Pp, = 1.06 Ibs. per sq.in. absolute 


pr 0.0822 Ib. per cu.ft. 


T»2 7 (7)(') 348°F. 
Pi pe 
| = _ 
Ve = v. (=) ‘3 
M/NT, 


The ambient properties of the air in region 2 are now 


absolute 


2,110 ft. per sec. 


) 


known. The properties in regions 3 and 4 were cal 
culated similarly. 

The properties of air computed for the four regions oj 
the wing are summarized in Table 1. 


Numerical Method of Solution for Two- Dimensional 

Heat Flow 

The preceding calculations show the method of ob 
taining the supersonic flow conditions about the wedge 
shaped wing. . The method of calculation of the two 
dimensional heat flow is presented here. Consider the 
wedge-shaped wing subdivided into a network of 113 
stations in a plane normal to the span, as shown in Fig 
10. We assume that the stations are connected by 
rods along which heat flow is concentrated between 
pairs of stations. Details of three typical stations and 
of the areas (or volume if the third dimension along the 
span is included) associated with these rods are shown 
in Fig. 26. Thus an internal station, H7, is shown in 
Fig. 2c and the area associated with it is represented 
by the shaded area centered on H7. A typical sur 
face station, F2, is shown in Fig. 26b, and, finally, the 
Fig. 26a. The two 


mid-chord station is shown in 


dimensional heat flow is then broken up into one 
dimensional heat flow between any one station and its 
closest neighbors along the net of connecting rods. 
For any pair of stations we used the differential equa 


tion for heat flow in finite difference form, given by 


Q = —kA (AT, Ax) AO (A-1 
Let us apply this to a typical internal station, to a 
typical surface station, and to the mid-chord station 
and use the symbols given in Fig. 26. 

The heat flow during 
H7 from station Hs 


(a) Typical Internal Station. 
a time interval A@ into station 


is simply 


QO = kRS2[(T ys T 7) / Si: |] 40 
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TABLE | 
Properties of Air in Regions Around Wedge-Shaped Wing Moving at Three Supersonic Speeds 


Angle 
Between 
Angle Wedge Shock 
Flight of Angle and In Turning 
Mach Attack, for cident Angle in 
Number, B Re Shock, Flow, Expansion, 
M (Deg.) gion w’ a’ € 
20 3.58 l 1.01 30.8 
2 9.18 
3 8.17 37.4 
4 9 18 
$.00 1.89 l 2.40 16.3 
2 9 18 
3 6.48 19.2 
} Y 18 
6.00 1.21 l 3.38 11.9 
2 9.18 
5 9.80 13.8 
4 9.18 


Altitude = 50,000 ft. 

Wing loading + = 100 lbs. per sq.ft 
Maximum wing thickness = 8 per cent 
Included wedge angle = 9.18 


where the spanwise length is taken as unity. The total 
heat flow in station H7 from all four surrounding sta 


tions 1s 
kA@ E T Tin 4 3 T T 
j=R ( 8 7) > ( “Aas 7) t 
¢ S, Hi H S, Ht H 7 
tr Tux) + sr Tz) | (A-2) 
37 7) 7 = - A-2 
5° 9 H7 Ss 17 H7 


The rise in temperature of the material associated 
with station H7 during the time interval 4@ is given by 


T 7’ — Tw = Q pcS,S» (A-3) 
where 7,7 and 77’ denote the temperatures at the 
beginning and the end of the time interval A@, respec- 
tively. We can eliminate 747 by proper choice of the 
time interval, such as by 


2k AO 1 +4 (=) 

i. pc," So 
By elimination of 77 from Eqs. (A-2) and (A-3), 
we get for the temperature, at the end of the time in- 


(A-4+) 


terval, 
Tur’ = (X22) [((Tue + Tus) + 7? (Voz + Tiz)] (A-5) 
where the thickness ratio, 7, is defined by 
rT = 3y/De (A-6) 
and 
| I 
A = (A-7) 


1+ (S/S)? 1+ 72 


(6) Typical Surface Station. The heat flow into a 


typical surface station in the time interval Aé is 


QO = hA(T., — T,) d0 





Ambient 
Pressure, 


Ambient 


Ambient Tempera Ambient 


Ambient Velocity, ture, p Density 
Mach v 7 Lbs. per p 

Number, (Ft. per F Sq.In Lbs. per 

M Sec Absolute Absolute Cu.Ft 
1.97 1,920 399 1.79 0.0121 
2 .ae 2,110 348 1.06 0 OO822 
1.70 1,760 $45 2.63 0.0159 
2.03 1,940 384 1.59 0.0112 
3.80 3,840 $26 2.21 0.0140 
1.54 $010 325 0). 847 0.00705 
3.53 3,750 171 3.14 0.0180 
1.19 3,900 364 1.28 0. 00950 
5.59 5,820 $52 2.75 0.0164 
6.95 6,020 308 0). 682 0.00597 
5.20 5,720 505 3.78 0.0202 
6.46 5,970 358 1.01 0.00763 
Upstream Conditions 

Temperature, 7) = 393°F. absolute 

Pressure, p = 1.69 lbs. per sq.in. absolute 

Density, pi = 0.0116 Ib. per cu ft 

Sound velocity, a@ = 971 ft. per sec 


where the bars refer to an average value of the variable 


during this time interval. The total heat flow into 


station F2, in Fig. 26b, is 


QO = hS3 (Tax Tyo) 40 + RAO X 


Dae = . a 
| - (153 - T ¥2) 7 ~ (TG2 rs) | (A-9) 
S) 2 So 


The temperature rise of station F2 is given by 


Ty’ = Ty = 20 pcS1Se (A-10 
where the factor of 2 is used because the surface station 
F2 contains only half as much metal as an internal 
station. Combining Eqs. (A-4), (A-9), and (A-10) and 


eliminating 7p, we get 


h Tr Tr , ry 
T 9" = rv ; TS3 ‘gee = T 2) 7 T 53 ie "To | T 
\ l = nN ns N)T 2 (A 11 
(c) Mid-Chord  Station..-The_ mid-chord — surface 


station, shown in Fig. 26a, is treated similarly to sur 
face station F2, except that the heat-transfer area was 


adjusted and the volume of material at station H1 


is one-quarter that of an internal station. We get 
V = h AdS3( i —_ Zn) T k AOS» S)) (Tye V1) 
(A-12 
T a’ Ta 1@) pcS)S» (A 13 
TH’ 2r[(h k) 7S3\ T,, T T Tye | 
(2r l vir \ l t 
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J A = absolute value of the path inclination with respect 
e CMMNMM1 horizontal plane 
p = air density 
¢ = absolute value of the angle of bank 
w = (V/R)cos@ = angular velocity Fl 
Discussing, for the sake of simplicity, a helicoidal unifor ner 
RIEF REPORTS of investigations in the aeronautical science * flight without sideslip and assuming that the resulting aerogy 
and discussions of papers published in the JOURNAL are jamie force lies in the plane of symmetry, the general vector straigh 
presented in this special department. The publication % ill be com condition of equilibrium of forces is 
pleted approximately 8 weeks after receipt of the material. The - 
hag ie ee eens 5 et : 4 = () straig 
Editorial Committee does not hold itself responsible for the opinions T+D+L+WrF ann 
expressed by the correspondents. Contributions should not exceed Admitting the thrust, as usual, parallel to the velocity, ¢ | 
: f : . a ‘ os ade Stralg 
800 words in length projection of all the forces on the axes of the principal trihedr ee 
of the path (tangent, principal normal, binormal) leads to th 
following scalar equations :* Vertic 
+ 
T — (1/2)CpoSV?2 + W sin 6 = 0 9 
\ertic 
. . : W V? l 
On the Speed of an Airplane in Steady Flight = cos 2.6 — - Cr pSV? sin ¢ = 0 a 
g Kk 2 Leve 
Angelo Miele (1/2)CipS V2 cos ¢ — Weos 6 = 0 di 
- , - url 
Professor of Mechanics, Escuela de Aviacion Militar, Cordoba, ’ hich tt saith — : fate) ind : ; 
Argentina, also, Assistant Professor, Aerodynamic Design, in which t ae RETO ieee, coe ficients Cz, ¢ pare nat independent - 
Escuela Superior de Aerotécnica, Cordoba, Argentina but related (incompressible flow) by the parabolic formula Straig 
September 8, 1950 Cp = Cop + (Cr?/rA.R.c) 5 
1 HAVE READ THE INTERESTING ARTICLE published by Freeman! eqs. (2), (3), (4), and (5) form a system of four equations wit! 9) 
about the velocity of an airplane in many conditions of un six variables hasic 
accelerated flight, and I have observed that the equations that ~. . e fli 
. . Cp, » Vi, 0, oN 6 coe 
give the speed can also be deduced from a more general formula pore ’ ate 
— . ° ra . . . “ef ° a 
containing, as particular cases, those obtained by Mr. Free which is determined if two parameters are given—-e.g., path in 
° ° ) 
man, clination and angle of bank 4 
j 
i ; 9). oF ee ee ee ditior 
As a matter of fact, all the cases treated by the author (straight From Eqs. (2) and (4) the lift/drag ratio is obtained, of th 
and circling level flight and climbing, descending, and vertical i ts W cos 6 . meth 
diving flight) can be studied under a general scheme, considering a ae ae pract 
that the most general steady motion of an airplane is a helicoidal ¢ : 7 ae , : = I 
motion with vertical axis, constant velocity, and constant an foams ge, (5) and (7 ), the conespentiag Hit cvelicient, ‘ 
: , ; aiid assur 
rotational speed. ? > 
I C1 Emaz gE * duce 
It is known that, for a given thrust and for an ideal ambient Ce - 2 is \ ae E , to th 
) i maz 
of constant air density, there is a triple infinity of possible heli- : ; ; Te 
coidal paths, each of which is determined if three independent Bap. (4) and (3) give the expression of the speed from 
oe se given (within certain limits, of course ) | f.. l E 2 aircr 
rhe following symbols are used: V= \ : i. L=t> v of pt 
‘ ( pppS | Ie 
A.R. = wing aspect ratio ‘ ; . ; . 
C ; 8 oa sia and Eqs. (3), (4), and (9) the radius of the cylinder to which th 
D = drag coefficien ee 
‘licoidal pz be M 
Cp, = parasite drag coefficient a ee ea whe 
CL = lift coefficient V? cos 6 7, cos 6 E *\ 
; : ; a - R= = Le aaa r= 0) 
Cre = V7A.R.eCp, = lift coefficient maximum for lift/drag g tang gCp, p Stan gl Emez 
3 ratio re a , : : ; is th 
Finally, angular velocity is deduced from Eqs. (3), (4), and 
D = drag 10 : : net 
. +> . \ ) 
e = airplane efficiency factor and 
I C./Cp = lift/drag ratio g CoS { Enc. \* 
; . : ¥ w = —tang = gtan¢ — . 
Ema V rA.R.-/4Cp, = maximum lift/drag ratio J . ra I 
F. = (W/g) (V?/r) = centrifugal force Ek \ is a 
F = V/Vwur = induced drag factor 1 + = - 1] i 
g E indi 
g = acceleration of gravity a. vee 
L = lift and lift-weight ratio from Eq. (4) 
n = L/W = lift/weight ratio r L/w a/ex ) ; 
_ z . - P == 4 = COs COs ¢ - 
r = R/cos? 6 = radius of curvature of the helicoidal path plan 
R = radius of the cylinder to which the helix belongs It is important to emphasize that Aer 
S = wing surface (1) Two possible angles of attack and, therefore, two possibi¢ : 
1 = thrust helicoidal paths exist for each couple of values 6, ¢. The greater _ 
or = bs , ~stty — - + : > rr ati ed 
7 = 7 + Wsiné@ = net thrust (+ for descending or climb velocity, the greater radius, and the smaller rotational speed Jou 
ing flight ) correspond to the smaller lift coefficient and vice versa Ma 
V = velocity (2) The exposed formulas contain explicitly the solution of the ‘ 
Var = maximum theoretical speed ideally obtainable, for a problem if thrust can be considered independent of speed; other 
given net thrust, without induced drags and with wise, they must be solved by trial s 
out compressibility effects 
‘ - : presst vility effect The sign preceding the third term of Eq. (2) must be considered posttty * 
i = weight for descending flight, negative for climbing flight i 
SUS 











respect to 


lal Unifor; 
INE aerody 
al vector; 


locity, th 
tl trihedrg 
ads to th 


ependent 


la 


ONS with 


path in 


*h the 


and 


(11 


bie 
ter 


ed 


he 
er 


READERS’ 


FORUM 80a 


TABLE 1 
Net Thrust, Lift-Drag Ratio, and Lift-Weight Ratio for Many Flight Conditions 


Flight Condition T. 
General helicoidal flight T «Wene — 
Straight level flight Tr 
Straight climbing flight T — Wsine 
Straight descending flight T + Wsine j 
Vertical diving flight T+w 
Vertical dive—no thrust W 
Level circling flight Tr 
Turning gliding flight IV’ sin @ 

, WV sin 6 


Straight gliding flight 


3) Results are closely subordinated to the validity of the 
basic hypothesis—namely, incompressible flow and flight lift 
coeflicient lower than limit lift coefficient up to which parabolic 
and experimental polars are in agreement. 

4) Results could also be extended to the study of flight con- 
ditions in compressible airflow if the dependence of the parameters 
of the parabolic polar on Mach Number is known and if iterative 


methods are applied. But this does not seem to be useful from a 


practical standpoint. 

In certain particular flight conditions, factors T,, E, and n 
assume the expressions given in Table 1. Some of them, intro- 
duced in Eq. (9), lead, through simple algebraic transformations, 
to the formulas obtained by Mr. Freeman in his paper. 

To conclude, when the assumption of independence of thrust 
from speed is correct (it is approximately the case of turbojet 
aircraft) Eq. (10) may be expressed as the product of two factors 
of physical significance* © 


V = Fi Var (13) 


where 


Ver = V2T,/CpppS (14) 
is the maximum theoretical speed ideally obtainable, for a given 
net thrust, with no induced drags and no compressibility effects, 


and where 


F; = (1/+/2)V 1 = V1 — (E/Enmaz.)? (15) 


isa reduction factor (<1) synthetically containing the effect of 


induced drags on speed. 
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An Analogy Between Open Soap Bubble and 
Compressible Flow 


V. G. Szebehely 
Associate Professor of Research, Department of Applied Mechanics, 
Virginia Polytechnic Institute, Blacksburg, Va. 


September 8, 1950 


i bees WELL-KNOWN ANALOGY between torsion and a soap bubble 
assumes infinitesimal bubble height. If the exact equation 
of an open bubble is compared with the basic equation for com- 
pressible flow, another analogy can be shown. The governing 
differential equation of an unloaded open bubble is 


O*s oz: \? O07: 072 Os Os 
-1 1 + + 
Ox? Oy 


oz \? ‘ 
; 1+ ) —2 - = () (1) 
oy? Ox Oxdy Ox OY 
and the basic equation for the steady, potential flow of a com- 


pressible, inviscid fluid 


el oi (*)'] + ela y (5*)'] + ee 
Ox? Ox oy? oy Oxdy Ox Oy 
The latter can be transformed by the 
v? = (0¢/0x)? + (0p/dy)? 
relation and using the a? = c? — 0? notation, into 
oO? O¢ \? oO? O¢ \? 0*o Od 
a [«° ’ (3) | . a [* . (3°) | aii it zi ni 
(3) 
Eq. (2) was obtained by combining the Euler equation, 
p grad (v?/2) + grad p = 0 
the continuity equation, 
div (pv) = 0 
the definition of the velocity of sound, 
c = Vdp/dp 
and introducing the velocity potential ¢. 
v, velocity vector (= grad @); and », 


Here, p is pressure; 
density. 

Eq. (3) shows complete analogy to Eq. (1) if a? = 1, which is 
the case of the “under transonic”’ flow, since the Mach Number 
for this-case is M = Vi — (1/c*). A detailed analysis-will be 
published at a later date. 
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On the Plastic Buckling of Plates According to 
the Flow Theory 


P. P. Bijlaard 
Cornell University, Ithaca, N.Y 
September 5, 1950 


- THIS JOURNAL OF July, 1950, Dr. Pearson! calculates the 
plastic buckling stress of a compressed plate that is simply 
supported at both unloaded sides, assuming plastic flow. He also 
assumes that during buckling the entire plate deforms plastically 
in analogy with the tangent modulus theory (Engesser load) for 
columns. 
Prager? had assumed that unloading occurs in part of the plate in 
analogy with the reduced modulus theory (von Karman load) 


In contrast, in a previous paper, Handelman and 


for columns. 

Handelman and Prager’s assumption led to disagreement with 
tests,* as is shown in Fig. 1, where the plastic buckling stress o-, 
is plotted against of/E (og is the elastic buckling stress and E is 

The circles refer to the experiments by 
In his paper, however, Dr. Pearson tries to 


the elastic modulus). 
Pride and Heimer]. 
show that, if the entire plate is assumed to remain plastic, a fair 
agreement with the tests is obtained, as is shown in his Fig. 10 
and also in Fig. 1 of this item. In Fig. 1 the curve according to 
the present writer’s theory* ® is also given, as calculated by Pride 
and Heimerl.* 

In the 
first place, the same case of plastic flow (Prandtl-Reuss body) 
and of the fully plastic plate was already considered in reference 
4 by the present writer. In Fig. 11 of that reference, the ; 
dashed curve from the left (Case 2, e = 0, pl.) refers to this case. 
This curve is, however, fairly close to the curve farthest to the 
left (Case 2, e = 0, el. + pl.), which was obtained by assuming 
partly elastic behavior. 
the same marked disagreement with Kollbrunner’s tests. 


Dr. Pearson’s paper is misleading in several respects. 


cond 
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For the material used in Pride and Heimerl’s tests, the cury, 
for plastic flow and fully plastic plate is given here in Fig. 1 by 
the curve denoted by (e = 0, pl.). It follows directly from Ey 
(53) of the writer’s 1938 paper,® or from Table 5 of reference 
according to which 


9 = oc/ox = 0.455 (VAD + B + 2F) 


if in reference 4, in Eqs. (21) or (21a) for A, B, D, and F, the ratiy 
8 = p»/p, between the principal compressive stresses is equated 
Furthermore, for plastic flow, in these equations ¢ ~ 


(E/E,) — 1 has to be equated to zero (E, is the secant Modulys 


zero. 


Here, too, the curve (e = 0, pl.) thus obtained for the flow theory 
is practically in the same disagreement with the Pride au 
Heimer] tests as is the Handelman-Prager curve (Fig. 1), 

Why is the curve obtained by Dr. Pearson for the same cage } 


better agreement with the tests? The reason is that he assume 
Poisson’s ratio v for the elastic deformations equal to 0.5, instea) 
of the more accurate value vy = 0.3 used by the present writer 
With »v = 0.5 in Eq. (1), the first factor to the right, representing 
0.5 (1 — v?), becomes 0.375 and the expression in parentheses 
follows from the above-mentioned equations of reference 4 with 
v = 0.5. This yields indeed the curve obtained by Dr. Pearson 

Although it would appear from Dr. Pearson’s paper that th 
better agreement with the tests is obtained by assuming the ep 
tire plate as plastic, about 85 per cent of this improvement js 
actually obtained by the arbitrary and incorrect assumption 
As stated on page 423 of Dr. Pearson's paper, he mak 
this assumption ‘‘for simplicity of algebra.”’ 


y = 0.5. 


It is obvious that, if Dr. Pearson wanted to show the improved 
agreement of his approach with tests, it is not permissible to 
introduce an arbitrary, simplifying assumption that accounts 
for 85 per cent of the apparent improvement. From a scientific 
point of view it is important to state that, actually, the assump 
tion of plastic flow and fully plastic plate leads to practically the 
same amount of disagreement with experiments as that of plastic 
flow and partly elastic plate. 

Dr. Pearson states, furthermore, that all treatments using 
deformation theories are open to objection, for which he refers toa 
1942 paper by Handelman, Lin, and Prager (actually published 
1947). 
been most rigorously treated by Handelman and Prager,? using 


He says that the problem of plastic plate buckling has 


the flow theory and that his paper is based on their results. 

In his 1938 papers, the present writer developed the general 
theory of plastic buckling of plates, including both the case ofa 
Hencky body (later referred to by other authors by the confusing 


name of ‘‘plastic deformation’’) and that of a Prandtl-Reuss body 
(plastic flow). 


as a Hencky body is not possible, since in extreme cases the ex 


He showed that, in general, a complete behavior 


ternal work done with plastic deformation would be negative 
or the work done during the total (elastic + plastic) deviatoric 
deformation would be zero.4 This objection is more convincing 
than that given by the above-mentioned authors, who state that 
Such dis 


In reference 


discontinuities in the plastic strain would occur. 
continuities also occur in other cases of plasticity. 
4, it was explained why, in the special case of concentric buckling, 
the material will nevertheless behave in accordance with the 
deformation theory. 

Furthermore, the writer showed in reference 4 that all later 
The 


treatment by Handelman and Prager is only rigorous as far as 


treatments are special cases of his earlier general ‘theory. 


their results are identical with those derived by the present writer 
in 1938 
cerned, which, for plastic flow, follow from his formulas by equat- 


i.e., as far as their excess stress-strain relations are con- 
ing his value e to zero. He considered this case of a Prandtl-Reuss 
body explicitly in his 1938 papers® but concluded that it yielded 
unsafe values as compared with the then available test results 
The newer investigations since then have confirmed thts con- 
clusion. 
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Note on Desirability of Completely Expanded 
Nozzles 


C. B. Cohen and P. J. Evans, Jr. 
Lewis Flight Propulsion Laboratory, N.A.C.A., Cleveland 


September 8, 1950 ’ 

N REFERENCE 1, the effect of velocity coefficient upon the 
| relative merit, in terms of thrust, of sonic nozzles, and com 
pletely expanded nozzles, is investigated. The conclusion is 
reached that, for overall pressure ratios less than approximately 
6, more thrust will be obtained with the sonic exit than with the 
completely expanded exit. There, for simplicity, it is arbi- 
trarily assumed that @ = @s, where @¢ is the ratio of actual exit 
velocity to the velocity corresponding to an isentropic expansion 
from the initial conditions to the actual exit pressure. The sub 
scripts s and c refer to the sonic exit and to the completely ex 
If there are no losses in the diverging 


panded exit, respectively. 
This statement follows from 


section, @- will be greater than @s. 
the physical meaning of ¢. The quantity ¢? is the ratio of actual 
to ideal kinetic energy. For the completely expanded nozzle, 
the ideal kinetic energy is greater than for the sonic nozzle 
Therefore, for the same value of @ the losses will be greater in the 
completely expanded nozzle. Thus, for a given ¢;, ¢- can be 
greater than, equal to, or less than ¢,, depending on the losses. 
In general, it appears that the losses in the diverging section will 
be excessive from the standpoint of good design if it be required 
that @ be as low as @ For example, Fig. 2, reference 1, shows 
that for a @ of 0.95 the pressure ratio for equal thrust is about 
5.75. In order for ¢ and ¢, to both equal 0.95, the total pressure 
loss in the diverging section would have to be about 12.6 per cent; 
but calculations (explained Jater in this letter) for this pressure 
ratio, for a diverging section having reasonably high losses, indi 
cate a total pressure loss of only about 5.6 per cent. Thus, the 
assumption of equal @ does not appear justified in general. It 
appears that a more correct procedure would replace Fig 2, 
reference 1, by a family of curves of P;/po versus ¢-, where ¢s 
Here P and p refer to total and static 
obtained di- 


is constant on each curve. 


pressures, respectively. These curves would be 


rectly from curves similar to Fig. 1, reference 1. However, to 
draw definite conclusions, @, must be known, as well as reason- 
able values of the @- corresponding to @ 

An analysis starting from the sonic section and considering 
only the losses downstream of this section does not require this 
information. A velocity coefficient, ¥, is defined as the ratio of 
the actual exit velocity to that obtained by expanding isentrop- 
ically from the actual sonic section (station 2) to the actual exit 
Then for each pressure ratio, P2/po, the thrusts for 
sonic and complete-expansion exits can be equated and a value 
The result of this procedure is 


pressure. 


for y explicitly determined. 
shown as the equal-thrust curve of Fig. 1 

In order to estimate on which side of the equal-thrust curve 
an actual nozzle will lie, a series of calculations were made for 


conical diverging exits having constant friction coefficient, / 


FORUM 
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P2/p) = 6, y = 13 
For the conical nozzle the effect of friction is introduced only 
through the one-dimensional equation relating areas and Mach 


Numbers: 


A Mo[1 — (yM22/2) f cot 0]% (t/T)2 |f 
A.  Ms\.1 — (yM;2/2) f cot 6} L(t/T)s 


where 
1 — (y/2) f cot 6 

eg = ‘ 3 ’ 
(y — 1), + [yv/(y — 1)] f cot 05 
(1/2) (vy + 1)/(7 — 1)] 

B= : 
1 — [y/(y — 1)|f cot @ 

t = static temperature 

T = total temperature 

y = ratio of specific heats 

@ = cone half angle 


This results in curves of P2/p) versus y for constant values of 


f cot @ as shown in Fig. 1. Since a friction coefficient of 0.005 








~ 


is higher than generally obtained in practice and a divergence 
half-angle of about 3° is as low as seems practical, the com- 
bination will yield an approximate maximum value of f cot 6 
of 0.10. 
losses of 5.6 per cent discussed in the first paragraph. 


This value was used in calculating the total pressure 
The 
0.10 curve shows that the completely expanded nozzle produces 
more thrust for values of 7?2/po greater than 2.3. The contra- 
diction between this value of 2.3 and the value of 6 arrived at in 
reference | can be attributed to the excessive losses in the diverg- 
ing section which are implicit in the assumption that @ = @s. 

There is a point that has not been considered—namely, the 
possibility of the existence of a maximum thrust at intermediate 
degrees of expansion. Therefore, calculations of thrust were 
made for conical nozzles having constant friction coefficient. 
Fig. 2, for a pressure ratio, P2/po, of 6, is a typical result. It 
shows the variation of gross thrust ratio, F;/ Fs, with e, the per 
cent expansion, defined by 


100 [(p2 — ps)/(p: — po) 


As one might expect, the maximum thrust occurs for nearly com 
From other 


e= 


plete expansion for all reasonable values of f cot @. 
similar plots, it has been observed that the optimum per cent ex- 
pansion increases as P?:/po increases. Thus, if maximum thrust 
is desired, the relative merit of sonic and complete expansion is 
more an academic question than a practical one, since the 
Maximum occurs somewhere between sonic and complete ex- 


pansion. 
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Increased Jet Thrust from Pressure Forces, Journal of the 
17, No. 7, p. 425, July, 1950. 
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Note on Base Pressure Measurements 


Dean R. Chapman 

Ames Aeronautical Laboratory, N.A.C.A. Moffett Field, Calif. 

September 5, 1950 

oe APPARENT INCONSISTENCIES in the existing base 
pressure data on bodies of revolution have not as yet been 

More- 


wherein his 


pointed out and adequately explained in these columns. 
over, in the recent communications of Gabeaud," ? 
analysis is compared with the free-flight experiments of Hill and 
Alpher,* one essential consideration seems to have been over- 
looked. The present note attempts to point out and offer ex- 
planations of these discrepancies. This note was prompted by 
the of some N.A.C.A. 
pressure‘ and by an excellent set of experiments recently reported 
by Kurzweg which did not become available early enough to per- 
N.A.C.A. 
Before discussing the 


declassification measurements of base 


mit comparison with the most recent unclassified 
No. 2137.8 


apparent inconsistencies, a few statements about some general 


measurements given in 7.N 


characteristics of base pressure, and the basis for such statements, 
are presented, since it is felt that inconsistencies have arisen in 
the past and can easily arise in the future by overlooking some of 
these characteristics. 

(1) Changes in type of boundary-layer flow can produce large 
changes in base drag; the effects uf variation in Reynolds Number 
are large in the laminar and transitional region but are small in the 
turbulent region. Measurements at low Reynolds Numbers (less 
than 2 X 10° based on body length) of the effect on base pressure 
of changes in Reynolds Number have shown that, with laminar 
flow approaching the base, the base drag easily can vary by a 
factor of 2 or greater. In this laminar region an increase in 
Reynolds Number results in an increase in base drag. The transi- 
tion to turbulent flow, in the cases measured thus far, is accom- 
panied by marked changes in the trailing shock configuration and 
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a decrease in base drag—the decrease being greater for gre 


ater 


angles of boat-tailing. In fact, for large angles of boat-t 


transition in certain cases* has resulted in a small net a 
acting on the base. For completely turbulent boundary-laye, 
flow, however, a different situation exists with regard to 4, 
variation of base pressure with Reynolds Numbers. The meas 
urements of reference 6, using different kinds of artificial] rough. 
ness, Show no appreciable effect on base pressure of variations; 
Reynolds Number between 2 X 10° and 1.6 X 107, Likewise 
the measurements of Faro® and of Hill and Alpher show no effect 
of Reynolds Number in the region between 5 X 10° and 1 x jy 

(2) Large changes in surface temperature can significantly alt 
the base pressure ‘ 
from approximately —40° to 350°F., as observed by Kurzweg 


The effect of surface temperature variatioys 


resulted in changes in base drag amounting to as much as 45 pe, 

per 
cent at M = 2.87. 
mation of base pressure shortly after combustion terminates jy , 


This effect would have a bearing on the estj. 


Changes in surface temperature an¢ 


rocket-propelled missile. nd 
in heat transfer can be looked upon as additional methods of 
changing the boundary-layer flow approaching the base. 

(3) The effect of body shape is large for length-diameter ratj 
less than about 4 but is small for length-diameter ratios greater thay 
At low and moderate supersonic Mach Numbers, the 


about 6. 
departure of flow conditions along the body surface from corre. 
sponding free-stream conditions results in a lower pressure and 
higher Mach 
much larger base drag for a short body of revolution than for g 


Number approaching the base. This results jy 


long one. A set of measurements at J = 2.46 showing large 
effects of body shape for length-diameter ratios less than unity js 
Some information is also available from 
the 


Large differences in base 


presented by Kurzweg. 


German wind-tunnel measurements’ and from Aberdeen 
firings of Charters and Turetsky.* 
drag for different body shapes, all having length-diameter ratios 
less than 4, are indicated in reference 6, where the observed dif. 
ferences amounted to several hundred per cent in a few cases 
For length-diameter ratios greater than about 6, however, no 
important effect of changes in the profile of a body of revolution 
has been found either in the N.A.C.A. measurements, the N.O.L 
measurements, the Aberdeen measurements, or the APL/JHU 


measurements. 


(4) The presence of tail fins can considerably change the base 
drag. The effect of tail fins can perhaps best be illustrated by 


reference to Fig. 1, which shows four independent sets of base 
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Figure /.- Comparison of experimental base pressure data 
for turbulent boundary -/ayer flow on bodies of revolution 
without boat tailing. 
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drag measurements presented in pressure coefficient form. All 
of these data are representative of the case of turbulent boundary- 
jayer flow approaching the base of relatively slender bodies of revo- 
istion having no boat-tailing. The APL/JHU data, indicatmg the 
greatest base drag, were taken on bodies having fins. The large 
difference is not attributable to Reynolds Number effects in view 
of]; nor, asshown in reference 6, is it attributable to interference 
in the wind-tunnel tests from supports or reflected bow waves. 
Moreover, the approximate similarity of surface temperature 
conditions in the two flight investigations, and also in the two 
wind-tunnel investigations, would eliminate temperature effects 
as the source of trouble in this case. In addition, the free- 
fight firings of Hill and Alpher at two different rates of deceler- 
ation indicate that deceleration rate is not responsible for the large 
difference. Partly on this basis of elimination of other sources 
of discrepancy, it has been concluded® that the observed differ- 
ence is due to the presence of tail fins on the APL/JHU missiles. 
These fins were 10 per cent thick, biconvex airfoil sections, of 
chord length equal to 1.5 body diameters. Simple calculations 
of the approximate average decrease in pressure over the rear 
surfaces of the tail fins show that the order of magnitude of this 
interference effect in the region of interest is ample to account for 
the difference. Subsequent wind-tunnel measurements of base 
pressure at M = 1.5 and M = 2.0 with and without these particu- 
lar fins have confirmed this. It may be noted that significant 
effects of tail fins, only.on a different configuration, are also 
evident in the measurements of reference 5. 

Turning now to the inconsistencies mentioned earlier, it is 
clear that the often-discussed discrepancy between various Ger- 
man wind-tunnel and flight measurements of base pressure” ® 
can be attributed to the differences in Reynolds Number and type 
of boundary-layer flow. A second point to be mentioned con- 
cerns an observation of Faro,’ who noted that the measured base 
pressure behind a hemisphere at J = 1.7 agreed very well with 
the free-flight data on the APL/JHU missile. This good agree- 
ment clearly is of no significance to the discrepancy between vari- 
ous base pressure measurements for the following reason: The 
hemisphere, being quite short, has considerably greater base 
drag than a slender body; also, the APL/JHU body with fins 
has considerably greater base drag than a slender body without 
fins. Base drag on the hemisphere and the missile with fins, 
therefore, might be expected to agree approximately with each 
other. 

With regard to the comparison between Gabeaud’s most re- 
cent theory? and the measurements of Hill and Alpher, it appears 
that this good agreement is fortuitous—that is, unless Gabeaud’s 
theory (the details of which are unavailable at present in this 
country) somehow accounts for the effect of four 10 per cent 
thick, biconvex tail fins. Similarly, the observed ‘‘knee’”’ in the 
APL/JHU base pressure curve and the particular values of the 
parameters f and \ introduced by Hill'® to match the APL/JHU 
data appear to be more characteristic of the particular tail 
geometry used rather than of some intrinsic behavior of 
base pressure in the range of Mach Numbers between 1.5 and 


2.0. 
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The Reflection of ‘‘Rectilinear’’ Waves in 
One-Dimensional Gas Flows* 


Lorenzo Poggi 
University of Pisa, Italy 
September 11, 1950 


iy A RECENT PAPER ON PULSE-JETS,' I described a method for 
studying one-dimensional nonsteady flows. This treatment 
was based on the use of particular solutions of the general equa- 
tions which govern such phenomena. One of these solutions 
(denoted in that investigation as the C-type solution) is of espe- 
cial interest, inasmuch as it provides, to close approximation, a 
means for representing what happens when a “‘rectilinear’’ wave 
is reflected at a wall. The present note is confined, in particular, 
to an illustration of such a solution and to its application in the 
case of a reflection process. 

As is well known, the flow of a gas in a one-dimensional duct 
under the hypothesis that the processes are adiabatic may be 


summarized by the set of equations?® 


Equation of Motion 


Ou Ou 2 Oa 
+ u =- a: 
ol Ox y¥- Il Ox 
Equation of Continuity (1) 
Oa Oa y¥-!1 Ou 
+ # = =e a: 
ol Ox 2 Ox 


wherein the following notation applies: 


x = distance along duct 

t = time 

u = local flow velocity 

a = local velocity of sound 

y = ratio of the specific heats 


This system of equations is obviously satisfied if we take as 
solution 
u = uo) 
a (A) 
a = ao\ 
that is, flow and state parameters that are uniform throughout the 


duct and constant with respect to time. 
Another well-known solution is 


2 x * aol 
“= : 
y¥7t 1 t+ to > (B) 
a =a = [(y — 1)/2]u 


where do and ¢ are constant. Solution (B) represents what may 
be called a ‘“‘rectilinear wave’’ because of the linearity of u and a 
with respect to x. 

Solutions of the two types can join up at their end points to 
produce complex wave setups such as, e.g., the one represented 
in Fig. 1, which shows an expansion wave (B) between two A-type 
regions, one of quiescent fluid on the left and one of uniform mo- 


* The author wishes to express his thanks to Dr. J. V. Foa, of Cornell 
Aeronautical Laboratory, Inc., for his interest in the translation and pres 
entation of this note and to Mr. Kramer, of the same Laboratory, for his 


careful translation 
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Fic. 1. Typical complex configuration made up of two A-type 
uniform waves and a B-type expansion wave, before reflection is 
initiated 
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Fic. 2. Approximate and tangent solutions for reflection of a 
B-type wave against the closed end of a duct. 
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Fic. 3. Configuration created in the duct during reflection at 
the closed end 


tion on the right. The expansion wave remains rectilinear but 
changes its slope with time because of the different speeds with 
which its end points are propagated toward the left. This speed, 
as well as the speed of the intermediate points, is always equal to 
u—a 

When the left-hand end point of the wave reaches the closed 
end of the duct, there begins the process of reflection which forms 
the object of consideration in this note 


OS wasn _ 
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The premise is made that, during this process of reflection, th 
state of the fluid is almost constant with respect to x (i.€., dgs 
Ox = (0) throughout a region adjacent to the wall. Substantiatjo, 
for this premise is obtained by analysis with graphical numericy 
methods, such as that of reference 3. Furthermore, the Velocity 
is obviously zero at the wall. Thus, one is led to look for a Soly 
tion of the set of Eqs. (1) for which just these properties are pa 
hibited. Such a search has led me to the third-type solution 


u = x/(t + t-) / 
a =a,[1 + (t/t-)J@ — Y/2( 


where a, and ¢, denote appropriate constants 

Although this solution is not a rigorous mathematica] fy 
simile of what takes place at reflection [even though it is an ex 
solution of the Eqs. (1)], nonetheless it furnishes a sufficient) 
approximate description of the phenomenon. — For the represent, 
tion to be rigorously exact, it would be necessary that the solutig; 
under consideration join up precisely with the arriving recj 
linear wave and that at each instant there should exist a pojn; 
of intersection between the two solutions (the a-values and , 
values of the B-type solution equaling those of the C-type) 

This state of affairs does not come about in an exact way, but 
the constants can be so selected that the C-type solution js 
“tangent” to the actual solution at the initial point of the refle 
tion process, thereby establishing the above-described conditig, 
not only when the incident wave just reaches the closed end by 
also an infinitesimal time interval later. 

More generally, it is possible to select the constants in such 
way that the above-mentioned condition is satisfied for two dij 
ferent time values 

Thus, letting ¢ = O at the initial instant of the reflection, th 
“‘tangency”’ requirement may be realized by assigning to th 
constant a, the value do of the originally quiescent fluid, and to tk 
constant /, the value 


t = [(vy + 1)/4] (L/ao) 


The symbol L denotes a characteristic length of the incident 
wave—naimely, the distance from the closed end to that point, 
the wave where the propagation velocity is zero 

The degree of approximation of such a configuration is « 
picted in Fig. 2, wherein the values of a that correspond to th 
incident wave and to the C-type configuration are represented by 
solid lines; the representation is made for those values of x 
which the respective w-values coincide. As may be seen, th 
two curves are tangent at / = 0. 

In this same figure, an analogous comparison (using dashe 
curves) is then carried out for the case in which the constant 
is assigned the value 0.5L/ay instead of 0.6L/ay [as would be o 
tained from Eq. (2) with y = 1.4]. The closeness of the curves 
is better on the average over the extensive time interval show 
when the value 0.5L /ay is used, even though near the origin they 
are less satisfactorily close. In practice, depending on 
amplitude of the range of values throughout which one desires t 
approximate the solution, as happens in many similar approx 
mation problems, one or the other (or intermediate) values maj 
be selected for the constant. 

Within such limits of approximation, immediately after th 
wave has reached the closed end, one finds, as shown in Fig 
three regions of the flow, one of type (C), one of type (B), an 
one of type (A) (constant flow velocity). This situation cot 


tinues until the intersection of the C-type and B-type solution: 
reaches the intersection of the B-type and A-type ones, which ts 


the right end of the incident wave, traveling in the opposite direc 
tion from the former. 


From this moment on (as graphical calculations carried out by 
means of the usual procedure* have confirmed ), a true and proper 
reflected wave is established. The intersection between thi 


reflected wave and the C-type wave (again calculable in an af 


proximate manner) moves toward the wall until it finally reache 
it, thereby eliminating the C-type configuration, which is thet 
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placed by a new region of zero velocity (A-type wave), com- 
ep 
oleting the reflection process. 
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Determination of Deflections in Redundant 
Structures 


HL. Cox and D. G. Sopwith 
Engineering Division, NationakPhysical Laboratory, England 
September 11, 1950 


HE POINT MADE BY B. A. Boley and R. H. Moore (Readers’ 
bi JOURNAL OF THE AERONAUTICAL SCIENCES, Vol. 
I7, No. 8, p. 526, August, 1950) seems not quite clear. For the 
determination of deflections in redundant structures, it has been 
wr normal practice to include the dummy loads P’ with the 
actual loading P and the redundancies F at the outset, to set up 
the expression for the strain energy U in terms of all three classes 
of load, and then to write 


oU/oF =0, 6 oU/oP’ 


This we believe to be the normal procedure in 


when P’ = 0. 
England. 

In computing 0OU//OP’, the redundancies F, as well as the ap- 
plied loads P, are regarded as constants, and in the derived for- 


mula for 6 the values of F are substituted from the results of 


equating OL /OF to zero. 

By this procedure the computation of the strain energy is 
carried out wholly for a simply stiff structure, the actual applied 
loads, the redundancies, and the dummy loads all being regarded 
as external loading. 

It seems that the normal procedure to which Messrs. Boley and 
Moore refer must be in some way different, but we cannot agree 
that the procedure outlined above is missing from the standard 


texts. 


On the Equations of Longitudinal Stability 


John W. Miles 


Associate Professor of Engineering, University of California at Los 
Angeles 
September 18, 1950 


en, THERE HAS BEEN some speculation as to which 

of the unsteady flow terms arising in the calculation of 
aerodynamic forces on an oscillating aircraft should be retained. 
The classical procedure! includes, in addition to the ‘“‘quasi- 
stationary”’ terms, only that term arising from the effect of wing 
downwash on the tail. 
that the calculation of pitch (or yaw) damping on a quasi- 
stationary basis may lead to serious errors,? and Statler has 
shown that unsteady flow effects may be important in other 
respects.3 


Nevertheless, it has been demonstrated 


However, an order analysis of the problem in terms of 
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its (dimensionless) characteristic parameters seems never to have 
been set forth, and it is to that end that the following is sub- 
mitted. 

Let x, y, z be a set of right-handed axes fired in direction and 
with origin at the center of gravity of the aircraft, being positive 
in the direction of undisturbed flight, to starboard and down, re- 
Uu and Uw, perturbation velocities (U’ is the flight 
velocity) along x and z and @ a rotation about y; the total 
B, the centroidal moment of inertia in pitch; and 6X, 6Z, 


spectively; 
m, 


mass; 
6M the perturbation forces due to the deviations u, w, 0. The 
equations of motion are 
m(0/odt)(Uu) = 6X (la 
m(0/odt)(Uw) = 6Z (1b) 
B(0?/ot?)0 = 6M (le) 


Then, introducing the notation of Durand (reference 1, pp. 130 
135; 
following dimensionless notation: 


however, Durand uses axes fixed in the body), we have the 


m m/(pSl) (2) 
kp = B (up.Sl*) = B/(ml?) (3) 
= (pl/U)(d0/dt) 4) 
6X = —pU*Sx (5a) 
6Z = —pU?Sz (5b) 
6M —pl*Slkam 5c) 
& = (1/2)5Cd = [Xu + (A/pu)Xu, Ju + 
[Xwo + (A/p)Xw, |w + [xe + (A u)xe |@ (6a) 
2 = (1/2)6Cr = [2uy + (A/u)ou,]u + 
[Zwo + (A/w)Swi lw + [%, + (A/u)Zq 16 (6b) 
) i 
m = — (1/2)(6Cu/ke) = [muy + (A/u)mu,|u + 
[wo + (A/p)Mw,|w + [ma + (A/pu)m, \6 (6c 
where / is a characteristic length. Eqs. (la)-(lc) then read 
hu+rx=0 (7a) 
Aw +s = 0 7b 
70 + um = 0 (7c) 


To obtain the stability quartic, we require the determinant, 


F(,, uw), of Eqs. (7a)-(7c) to vanish. Arranging this determinant 


in powers of \ and yw, we have terms of order uw" & ~» ~? and 
* % 2% 1° Now the order of magnitude considerations that dis 
tinguish dynamic stability from flutter, where such distinction is 
possible, require « to be large and \ to be of order one. (In this 
connection, we remark that, for harmonic motion, \ uk, where 
k is the reduced frequency introduced in flutter work.) Thus, 
if uw is sufficiently large to justify the original neglect of terms of 
order (A/u)? in Eqs. (6a)-(6c), which would have contributed 


terms of order »~! to the end result, we may consistently retain 


only terms of order uw! and y® in the determinant of Eqs. (7a 
(7c), whence we obtain 
F(X, w) = Ad\* + BAF + Cl? + DA +E 8 
A=1+4+ 0(z 9a 
B (Xu. + Zu, + mo,) + O(u 9b 
C = pme, + 
x 2 sm “4, Z as x ,) 
de ‘ i. 
Uy Wo Wy Ww, 8 uo 9; 
xm 
+ Ou 9c 
u, % 
xm com xsm 
D=u 4. 7" 4. 
Uy B Wo 4 Uo Wo 9 
x sm x om 
+ 9d ) 
Uo W, A Uy; Wy Oo 
xom 
E= M Qe 
Uy Wo Oo 
where the parenthetic sets of four and six variables denote Jacob 
ians; e.g., 
Su tu 
xs ™ 10 
mt lSuq Sue 20 
Uy We 0 . : 
M uo Mw MA, 
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Eqs. (9a)-(9e) simplify considerably for the case of free (power- 
off, no autopilot) flight at small angles of attack (a), in which 
case (due to our choice of axes fixed in direction) 


(0/Ow)o = (0/00)o = (0/0a) (11) 

Also, we may prove that 
(0/08); = (0/0qg) + (0/0a) (12) 
(0/Ow), = (0/0a@) (13) 


where (0/0q) implies steady flight at constant angular velocity 
(e.g., along a curved flight path, as in a pull out) and @ denotes a 
vertical acceleration, positive down. In view of the assumption 
of small angle of attack, we also may neglect drag terms relative 
to lift and moment terms, and Eqs. (9a)—(9e) reduce to 


Az 1 (14a) 
B = (Sa + mq + Ma) (14b) 
C = ume + (; a (14c) 
ad 
eal ‘i ag 4 “ad Z ”) (14d) 
uo @ Ug a@ q 
E=0 (14e) 


AERONAUTIC 
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It follows that, in this case at least, the only unsteady flow in 
of importance is mq. 

We should expect this last conclusion to be equally Valid, ins 
far as aerodynamic (as opposed to autopilot) terms alone » 
concerned, for power-on flight, but the neighborhood of stal 
flight may not be included without further investigation. Mor 
over, in interpreting Eqs. (14a)—(14e), it should be remarked th 
Zq, Za, Mg, and me, when properly calculated, all include doy 
wash effects of wing on tail.‘ 

A similar analysis of the equations of lateral stability could 
carried out, although it would be necessary to use axes fixed 
the body and to discuss the effects of yawing the wake. 
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The Transient Temperature Distribution in a Wing Flying at 
Supersonic Speeds 


(Continued from page 807) 


Final Calculation Equations for Two-Dimensional Heat 
Flow 


For the dimensions of the wedge-shaped wing with a 
thickness ratio of 8 per cent, we get, for the various 
dimensions and ratios, 


S; = 0.03142 ft. 
So = 0.3928 ft. 
S3 = 0.3941 ft. 
r = 0.08000 
h = 0.9935 


Substitution of these values results in the following 
final equations used in the actual calculations of two- 
dimensional heat flow in the wing: 


(a) Typical Internal Station: 


T'y7' = 0.4986(Ty6 + Tus) + 0.0032(T G7 + Ty 
(A-15 
(b) Typical Surface Station: 
Ty2' = 0.03132 (h/k) (Tan — Type) + 0.99357 p3 + 
0.006472 + 0.00017 ;. (A-16 
(c) \lid-Chord Station: 
Tw’ = 0.06264 (h/k) (Tax — Ty) + 
1.9871 Ty —— 0.9871 Tn (A-l7 
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ELASTIC STOP NUT CORPORATION OF AMERICA 
ELECTROL INCORPORATED 
ESSO STANDARD OIL COMPANY 
ETHYL CORPORATION 
FAIRCHILD CAMERA & INSTRUMENT CORPORATION 
FAIRCHILD AERIAL SURVEYS, INC. 
FAIRCHILD ENGINE AND AIRPLANE CORPORATION 
AL FIN DIVISION 
FAIRCHILD AIRCRAFT DIVISION 
FAIRCHILD ENGINE DIVISION 
FAIRCHILD GUIDED MISSILES DIVISION 
FAIRCHILD PERSONAL PLANES DIVISION 
NEPA DIVISION 
STRATOS DIVISION 
FEDERAL TELEPHONE AND RADIO CORPORATION 
6 ELECTRIC MANUFACTURING COMPANY 


FEDERAL TELECOMMUNICATION LABORATORIES, 
INC. (PARTIAL SUBSIDIARY) 
bag TONE TIRE & RUBBER COMPANY OF CALI- 
FLETCHER AVIATION CORPORATION 
FLIGHT SAFETY FOUNDATION 
GENERAL ELECTRIC COMPANY 
AIRCRAFT GAS TURBINE DIVISION 
AVIATION DIVISIONS 
GLOBE CORPORATION, AIRCRAFT DIVISION 
THE B. F. GOODRICH COMPANY 
THE GOODYEAR TIRE & RUBBER COMPANY 
GOODYEAR AIRCRAFT CORPORATION 
GRUMMAN AIRCRAFT ENGINEERING CORPORATION 
W. & L. E. GURLEY 
HARVEY MACHINE COMPANY, INC. 
INDUSTRIAL SOUND CONTROL 
INSURANCE COMPANY OF NORTH AMERICA COM- 


PANIES 
INSURANCE COMPANY OF NORTH AMERICA 


PHILADELPHIA FIRE AND MARINE INSURANCE — 


COMPANY 
INDEMNITY INSURANCE COMPANY 
AMERICA 
THE INTERNATIONAL NICKEL COMPANY, INC. 
IRVING AIR CHUTE COMPANY, INC. 
JACK & HEINTZ PRECISION INDUSTRIES, INC. 
JOHNS-MANVILLE SALES CORPORATION 
EARLE M, JORGENSEN COMPANY 
WALTER KIDDE & COMPANY, INC. 
a INSTRUMENT DIVISION, SQUARE D COM- 


OF NORTH 


LAVELLE AIRCRAFT CORPORATION 
LEAR INCORPORATED 

THE LIQUIDOMETER CORPORATION 
LOCKHEED AIRCRAFT CORPORATION 

AIRQUIPMENT COMPANY (SUBSIDIARY) 

LOCKHEED AIR TERMINAL (SUBSIDIARY) 

LOCKHEED AIRCRAFT OVERSEAS CORPORATION 

(SUBSIDIARY) 

LOCKHEED AIRCRAFT SERVICE, INC. (SUBSIDIARY) 
LONGINES-WITTNAUER WATCH COMPANY, INC. 
MARMAN PRODUCTS COMPANY, INC, 

THE GLENN L. MARTIN COMPANY 


McDONNELL AIRCRAFT CORPORATION 
MELETRON CORPORATION 
MINNEAPOLIS-HONEYWELL REGULATOR C0 
BROWN INSTRUMENT COMPANY (SUBSIDIARY) 
NATIONAL CREDIT OFFICE, INC. a 
NORTH AMERICAN AVIATION, INC. | 
NORTH AMERICAN AVIATION, INC., OF TEXAS» 
NORTHROP AIRCRAFT, INC. 2 
NORTHWEST AIRLINES, INC. 
PAN AMERICAN WORLD AIRWAYS, INC. 
THE PARKER APPLIANCE COMPANY ; 
PESCO PRODUCTS DIVISION, BORG-WARNER 
PORATION ; 
PHILLIPS PETROLEUM COMPANY 
PIASECKI HELICOPTER CORPORATION 
THE PURE OIL COMPANY 
RADIOPLANE COMPANY 
REPUBLIC AVIATION CORPORATION 
E. V. ROBERTS AND ASSOCIATES 
A. V. ROE CANADA LIMITED 
JOHN A. ROEBLING'S SONS COMPANY 
ROHR AIRCRAFT CORPORATION ; 
RYAN AERONAUTICAL COMPANY x 
SAVAL, INC. 
SCHRILLO AERO TOOL ENGINEERING COMPANY 
SER VOMECHANISMS, INC, ; 
SHELL OIL COMPANY, INC. 
SIMMONDS AEROCESSORIES, INC. 
SOCONY-VACUUM OIL COMPANY, INC. 
SOLAR AIRCRAFT COMPANY 
SPERRY GYROSCOPE COMPANY, DIVISION OF Til 
SPERRY CORPORATION i! 
SQUARE D COMPANY 
STANDARD OIL COMPANY OF CALIFORNIA 
STANDARD OIL COMPANY (INDIANA) Ls 
THE STEEL PRODUCTS ENGINEERING COMPANT, 
STURGESS, INC. 
SUMMERS GYROSCOPE COMPANY 
TELECOMPUTING CORPORATION 
THE H. I. THOMPSON COMPANY 
THOMPSON PRODUCTS, INC. 
TINNERMAN PRODUCTS, INC. 
TOOLKO ENGINEERING COMPANY 
TRANS WORLD AIRLINES, INC. 
UNION CARBIDE AND CARBON CORPORATION 
BAKELITE CORPORATION 
ELECTRO METALLURGICAL DIVISION 
HAYNES STELLITE DIVISION 
THE LINDE AIR PRODUCTS COMPANY 
NATIONAL CARBON DIVISION 
UNITED AIR LINES, INC. 
UNITED AIRCRAFT CORPORATION 
CHANCE VOUGHT AIRCRAFT DIVISION 
HAMILTON STANDARD DIVISION 
PRATT & WHITNEY AIRCRAFT DIVISION 
SIKORSKY AIRCRAFT DIVISION 
UNITED AIRCRAFT EXPORT CORPORATION 
UNITED AIRCRAFT SERVICE CORPORATION 
UNITED STATES AVIATION UNDERWRITERS, INC. 
VARD, INC 
VICKERS, INC. 
WESTERN GEAR WORKS 
WESTINGHOUSE ELECTRIC CORPORATION 
AVIATION GAS TURBINE DIVISION 
INDUSTRIAL ELECTRONICS DIVISION 
MARINE & AVIATION SALES DEPARTMENT 
SMALL MOTOR DIVISION, AVIATION SALES 
WESTON ELECTRICAL INSTRUMENT CORPORATION 
WYMAN-GORDON COMPANY 
YOUNG RADIATOR COMPANY 























